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A  STUDY  OF  THE  CANOPY  DESICN  FOR  THE  ADVANCED  ATTACK 


HELICOPTER  BY  USE  OF  COMPUTER  GRAPHICS 


INTRODUCTION 

The  US  Array  Human  Engineering  Laboratory  (USAHEL)  has  developed  a 
package  of  computer  programs  for  attack  helicopter  canopy  design.  The 
programs  compute  point-wise  measures  of  three  important  design  factors  for 
the  transparent  surfaces  of  an  attack  helicopter  canopy  design.  These 
factors  are  (1)  the  internal  glare,  (2)  the  external  glint,  and  (3)  the 
optical  distortions  exhibited  by  the  canopy  design.  This  work  was 
undertaken  at  the  request  of  the  Project  Manager's  Office,  USA  Aircraft 
Survivability  Equipment  (ASE),  and  is  part  of  an  effort  by  HEL  to  determine 
optimum  canopy  designs  for  the  Model  209  AH-1S  Cobra  Helicopter  and  the 
Model  YAH-64  Advanced  Attack  Helicopter  (AAH). 

A  current  problem  in  military  helicopter  design  is  the  choice  of  an 
optimum  configuration  for  the  transparent  surfaces  of  the  helicopter 
canopy.  The  design  should  satisfy  several  conflicting  requirements.  First, 
the  canopy  should  reflect  a  minimum  of  solar  glint  to  external  observers 
during  daytime  operations.  The  design  should  exhibit  a  minimum  of  internal 
reflections  of  the  internal  instrument  lights  and  external  light  sources 
during  nighttime  operations.  Finally,  the  canopy  design  should  not  distort 
the  pilot's  visibility  of  distant  scenes  during  flight.  The  need  for  these 
requirements  has  become  apparent  during  the  recent  history  of  attack 
helicopter  development. 

The  Model  209  AH-1G  COBRA,  a  predecessor  to  the  AAH,  was  equipped  with 
a  compound-curved  canopy.  The  aircraft  reflected  a  continual  solar  glint 
when  operated  in  daylight  under  clear  skies.  The  glint  was  highly  visible 
and  enabled  distant  ground  observers  to  readily  detect  and  track  the 
helicopter,  thereby  placing  it  at  a  tactical  disadvantage. 

A  flat  plate  canopy  (FPC)  design  was  developed  for  the  later  model 
COBRA,  the  AH-1S,  and  the  AAH.  The  design  reduced  the  solar  glint  to 
momentary  flashes  which  occurred  at  a  limited  number  of 
observer-aircraft-sun  angles.  A  moving  aircraft  no  longer  produced  the 
continual  solar  glint  which  was  present  on  the  earlier  compound-curved 
canopy  designs. 

However,  the  internal  surfaces  of  the  FPC  acted  as  mirrors  reflecting 
virtual  images  of  external  light  sources  during  nighttime  operations.  HEL 
has  shown  by  computer  analysis  (4,6),  that  these  reflections  can  occur  on 
most  of  the  canopy  surfaces  and  for  a  wide  range  of  external  source 
locations.  In  certain  lighting  situations,  these  virtual  images  of  ground 
level  lights  so  disorientated  the  pilot  that  he  could  not  easily 
discriminate  between  the  light  sources  on  the  ground  and  their  reflections 
from  the  canopy  surfaces.  The  use  of  the  FPC  as  a  solution  to  daytime  glint 
problems  resulted  in  a  potential  safety  hazard  during  nighttime  flight. 
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A  low  glare  canopy  (LGC)  design  was  developed  for  the  Model  YAH-64 
AAH.  The  design  has  little  solar  glint  and  lower  Internal  glare  than  the 
FPC  design.  The  LGC  consists  of  flat  transparent  panels  on  the  front,  the 
same  as  on  the  FPC,  but  the  top  and  sides  are  slightly  bowed  in  the  shape 
of  simple  cylindrical  segments.  HEL  recommended  a  similar  design  with 
several  novel  features  (6).  The  question  of  interest  is  how  effective  is 
the  low  glare  design  in  controlling  glint  and  internal  glare,  and  is  it 
possible  to  develop  an  optimal  design. 

The  computer  programs  developed  by  HEL  have  been  applied  to 
perturbations  in  the  present  LGC  design  for  the  Model  YAH-64  AAH.  The 
programs  computed  the  internal  glare,  external  glint  and  optical 
distortions  as  a  function  of  the  displacements  of  the  top  and  side 
cylindrical  panels  from  the  vertice  planes.  This  report  presents  the 
methods  and  results  of  the  study. 


METHOD 

A  set  of  computer  programs  was  written  to  compute  point-wise  measures 
of  three  important  design  factors  for  the  transparent  surfaces  of  a 
helicopter  canopy  design.  The  factors  are  (1)  the  internal  glare,  (2)  the 
external  glint,  and  (3)  the  optical  distortions.  A  description  of  each 
program  and  the  computing  method  employed  follows  along  with  a  discussion 
of  the  program's  input  data. 


Internal  Glare 

The  computer  program  (see  Appendix  C)  determines  all  points,  spaced  at 
equal  angular  increments,  where  the  pilot  can  see  reflections  on  the  canopy 
surface  of  external  light  sources.  The  program  computes  the  coordinates  and 
reflectances  of  the  reflection  points  and  the  directions  of  the 
corresponding  light  sources.  The  program  is  an  extension  of  ray  tracing 
techniques  developed  earlier  (4,5),  and  computes  secondary  and  higher  order 
reflections  as  well  as  primary  orders.  The  results  determine  the  portions 
of  the  canopy  surface  which  are  open  to  possible  reflections.  The 
reflection  points  are  plotted  on  side,  top,  and  front  views  of  the  canopy 
frame.  The  total  number  of  reflection  points  are  printed  out  for  each 
window. 


External  Glint 

The  computer  program  (Appendix  C)  computes  the  angular  direction  of 
the  specular  reflections  of  solar  glint  from  the  transparent  surfaces  of 
the  helicopter  canopy.  The  angular  directions  of  the  solar  glint  are 
computed  for  any  specified  range  of  sun-aircraft  angles. 
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The  program  first  determines  a  grid  of  equally  spaced  surface  points 
on  the  windows  of  the  canopy.  The  program  then  computes  the  direction  of 
the  specular  reflection  at  each  surface  point  for  each  of  the  specified  sun 
positions  (see  Appendix  A).  The  directions  are  next  plotted  as  points  on  a 
rectangular  plot  as  a  function  of  elevation  and  azimuth.  The  collection  of 
such  reflection  points  for  a  specified  sun-position  shows  the  range  of 
directions  from  which  the  aircraft  can  be  seen  by  distant  observers.  A 
collection  of  such  reflection  points  are  plotted  for  each  of  the  specified 
sun-positions.  The  position  of  the  sun  and  the  reflection  direction  are 
measured  as  elevation  and  azimuth  angles  in  the  coordinate  system  of  the 
aircraft  (see  Appendix  A). 


Optical  Distortions 

The  computer  program  (Appendix  C)  computes  the  light  transmittance  and 
optical  distortion  for  equally  spaced  surface  points  on  the  transparent 
surfaces  of  the  canopy.  The  program  first  determines  a  grid  of  equally 
spaced  surface  points  on  the  windows  of  the  canopy.  The  program  then 
computes  the  light  transmittance  and  optical  distortion  for  each  surface 
point  as  seen  from  the  pilot's  nominal  eye  position  (see  Appendix  B). 

The  program  computes  the  light  transmittance  for  the  straight  line  ray 
from  the  surface  ppint  to  the  pilot's  eye.  The  ray  is  traced  backwards 
through  the  canopy  surface  to  the  outside.  The  transmittance  is  the  product 
of  the  tiansmittances  at  the  inside  and  outside  surfaces,  and  that  not 
absorbed  by  the  window  material. 

The  optical  distortion  is  computed  for  each  surface  point  as  the 
square-root  error  of  the  lateral  magnification.  A  family  of  rays  is 
constructed  for  a  4-degree  cone  of  vision  centered  about  the  principal  ray 
from  the  surface  point  to  the  pilot's  eye.  The  rays  are  chosen  in  such  a 
manner  that  the  solid  angles  contained  by  adjacent  rays  are  equal.  The 
lateral  magnification  is  computed  for  each  ray  in  the  cone  by  tracing  the 
ray  backwards  through  the  surface  and  then  computing  the  angle  between  the 
incident  and  transmitted  rays.  The  square-root  error  is  computed  as  the 
square  root  of  the  sura  of  the  squares  of  the  differences  between  the 
magnification  for  the  principal  ray  and  each  ray,  divided  by  the  number  of 
rays  in  the  cone.  The  computed  measure  is  equal  to  zero  for  a  surface  which 
has  no  distortions. 

The  program  plots  the  computed  values  of  the  transmittance  and  optical 
distortion  on  side,  top,  and  front  views  of  the  canopy  frame.  The  plots 
allow  the  user  to  determine  which  areas  of  the  canopy  distort  the  pilot's 
distant  vision. 


Program  Input  Data 


The  canopy  frame  positions  and  the  shape  of  the  transparent  windows 
are  specified  as  input  to  the  programs.  The  transparent  surfaces  of  the 
canopy  design  are  specified  as  a  set  of  planar  and  cylindrical  surfaces  and 
their  corresponding  edge  vertices.  Each  planar  surface  is  specified  by  the 
coordinates  of  its  edge  vertices  and  the  consecutive  order  in  which 
adjacent  vertices  are  listed.  A  cylindrical  surface  is  specified  by 
cylindrical  parameters  and  the  consecutive  sequence  of  the  edge  vertices 
and  their  coordinates.  The  cylindrical  parameters  are  (1)  the  origin  point 
on  the  cylindrical  axis,  (2)  the  directional  cosines  of  the  axis,  and  (3) 
the  radius  of  the  cylinder.  The  edges  of  the  cylindrical  surface  are 
assumed  for  simplicity  to  be  curvilinear  lines  which  project  into  straight 
lines  when  the  cylinder  is  transformed  into  a  flat  plane. 

The  programs  are  tailored  to  the  YAH-64  AAH  low  glare  canopy  design, 
but  are  applicable  to  any  helicopter  given  the  canopy  frame  positions  and 
the  parameters  of  the  transparent  surfaces.  The  programs  compute  the 
measures  of  the  design  factors  as  functions  of  the  displacement  of  the 
cylindrical  windows.  The  displacement  is  measured  along  the  cylindrical 
radius  normal  to  the  plane  of  the  frame  vertices,  and  is  the  distance 
between  the  plane  and  the  cylindrical  side. 
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RESULTS 


The  computed  results  for  the  three  design  factors  are  shown  in  the 
figures  for  pertinent  combination  of  the  side  and  top  window  displacements 
on  the  Model  YAH-64  Advanced  Attack  Helicopter  (AAH).  The  results  are 
treated  separately  in  the  following  sections. 


Internal  Glare 

Figures  1  through  14  show  where  internal  glare  can  occur  on  the  canopy 
surfaces.  The  figures  show  the  top,  side,  and  front  views  of  the  canopy 
frame.  The  figures  are  paired  for  a  given  combination  of  the  window 
displacement  values.  The  first  figure  of  a  pair  shows  the  possible 
reflection  points  on  the  surface.  A  numeral  zero  or  one,  corresponding  to 
the  associated  reflectance  value,  is  shown  at  the  sjrface  position.  The 
second  figure  shows  the  entry  positions  of  those  rays  which  are  reflected. 
All  but  the  last  set  of  figures  correspond  to  a  2-degree  by  2-degree 
viewing  increment.  The  last  set  corresponds  to  a  0.5-degree  by  0.5-degree 
increment.  The  figures  show  primary  reflection  points.  The  secondary  and 
triplex  reflection  points,  computed  by  the  same  method,  are  discussed  in 
the  DISCUSSION  section. 

Figures  1  through  10  show  the  internal  glare  for  increasing  values  of 
the  side  window  displacement.  The  displace  icnt  is  measured  along  the 
cylindrical  radius  in  a  direction  normal  to  the  plan  of  the  vertices,  and 
is  the  distance  between  the  plane  and  the  cylindrical  sides.  The  top  window 
has  the  same  form  as  for  the  LGC  design.  The  figures  show  that  the  internal 
glare  on  the  side  windows  decreases  as  the  side  displacement  is  increased. 
The  figures  show  no  presence  of  internal  glare  on  the  side  windows  for 
displacement  larger  than  4  Inches. 
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Figure  la.  Primary  reflection  points  for  0.5  inch  side  window  displacement,  top  window 
present  design,  2°  by  2°  viewing  increment. 
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Figure  3a.  Primary  reflection  points  for  1.5  inch  side  window  displacement,  top  window 
present  design,  2°  by  2°  viewing  increment. 
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Figure  5a.  Primary  reflection  points  for  2.5  inch  side  window  displacement,  top  window 
present  design,  2°  by  2°  viewing  increment. 
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Figure  5b.  Entry  positions  for  reflected  rays  of  Figure  5a 
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Figure  7a.  Primary  reflection  points  for  3.S  inch  side  window  displacement,  top  window 
present  design,  2°  by  2°  viewing  increment. 
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Figure  7b.  Entry  positions  for  reflected  rays  of  Figure  7a. 
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Figure  8a.  Primary  reflection  points  for  4  inch  side  window  displacement,  top  window 
present  design,  2°  by  2°  viewing  increment. 
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Figure  1  .a.  Primary  reflection  points  for  0.5  inch  top  window  displacement  rotated  90°, 
side  window  present  design,  2°  by  2°  viewing  increments. 


30 


Figure  1 2a.  Primary  reflection  points  for  1 .0  inch  top  window  displacement  rotated  90°, 
side  window  present  design,  2°  by  2°  viewing  increment. 
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Figure  1 2b.  Entry  positions  for  reflected  rays  of  Figure  1  2a. 
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Figure  13a.  Primary  reflection  points  for  1 .5  inch  top  window  displacement  rotated  90°, 
side  window  present  design,  2°  by  2°  viewing  increments. 
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Figure  14a.  Primary  reflection  points  for  4.0  inch  side  window  displacement  1  5  inch  too 
window  displacement  rotated  90°  0.5°  by  0.5°  viewing  increment  P 
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The  changes  in  the  internal  glare  as  a  function  of  the  top  window 
displacement  were  computed  next.  The  computations  were  made  first  for  the 
configuration  of  the  LGC  design;  i.e.,  the  cylindrical  axis  of  the  top 
window  was  kept  parallel  to  the  pitch  axis  of  the  helicopter.  The 
computations  showed  little  reduction  in  the  number  of  reflection  points 
with  increased  displacement.  Figures  1  through  10  show  that  the  reflections 
occur  on  the  sides  of  the  top  window.  The  orientation  of  the  sides  was  not 
altered  by  the  window  displacement.  It  became  apparent  that  the  left  and 
right  sides  of  the  top  window  needed  to  be  bowed  in  to  decrease  the 
internal  glare.  This  was  accomplished  without  the  use  of  compound  shapes, 
by  rotating  the  cylindrical  axis  of  the  top  window  90-degrees  for  alignment 
with  the  longitudinal  axis  of  the  helicopter. 

Figures  11  through  13  show  the  internal  glare  as  a  function  of  the 
displacements  of  the  top  window  in  the  rotated  configuration.  The 
displacement  is  measured  along  the  cylindrical  radius  in  a  direction  normal 
to  the  plane  of  the  vertices,  and  is  the  distance  between  the  plane  and  the 
cylindrical  sides.  The  side  windows  have  the  same  form  as  for  the  LGC 
design.  The  figures  show  that  the  internal  glare  on  the  top  window 
decreases  as  the  top  displacement  is  increased.  The  glare  on  the  top  window 
is  absent  when  the  displacement  is  equal  to  or  larger  than  1.5  inches. 

Figures  14a  and  14b  show  the  internal  glare  for  one  combination  of 
side  and  top  window  displacements.  The  side  windows  are  displaced  4  inches 
and  the  top  window  is  displaced  1.5  inches  in  the  rotated  mode.  The  few 
reflection  points  which  are  possible  (for  a  0.5-degree  by  0.5-degree 
increment)  occur  on  the  upper  edges  of  the  side  windows.  Figure  14b  shows 
that  the  reflection  points  are  generated  by  light  rays  which  enter  at  the 
lower  front  corner  of  the  forward  side  windows. 
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External  Glint 


The  results  of  the  computations  for  the  external  glint  are  shown  in 
figures  15  through  21.  These  figures  were  drawn  on  a  Calcomp  plotter  which 
is  a  perpherial  unit  of  the  CDC  7600  computer  system  used  for  computations 
here  at  the  US  Army  Aberdeen  Proving  Ground  facilities.  The  figures  occur 
in  sets  and  show  the  results  for  various  combinations  of  side  and  top 
window  displacements.  The  first  figure  of  a  set  shows  the  canopy  surface 
points  used  in  the  reflection  computations  on  side,  top,  and  front  views  of 
the  canopy  frames. 


The  remaining  figures  of  the  set  show  plots  of  the  computed  reflected 
directions  for  the  side  windows,  the  top  window  and  all  windows  combined. 
The  directions  of  reflected  rays  are  plotted  as  points  according  to  their 
elevation  and  azimuth  angles.  Each  point  on  a  plot  gives  the  angles  for  a 
ray  reflected  from  one  of  the  surface  points  on  the  frame  sketch.  The  locus 
of  such  points  for  a  specific  sun-angle  forms  a  line  on  a  plot.  The  line 
shows  all  of  the  angles  at  which  solar  glint  from  the  window(s)  is  visible 
at  that  sun-angle. 
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Figure  15a.  Glint  surface  points  for  side  and  top  windows  of  present  design. 
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Figure  16a.  Glint  surface  points  for  side  window  1.0  inch  displacement,  top  window 
0.5  inch  displacement  in  rotated  mode. 
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Figure  1 7a.  Glint  surface  points  for  side  window 1 .0  inch  displacement,  top  window 
1.0  inch  displacement  in  rotated  mode. 
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Figure  1 7b.  Glint  from  top  window  and  all  windows  combined. 
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Figure  1 8a.  Glint  surface  points  for  side  window  1 .0  inch  displacement,  top  window 
1.5  inch  displacement  in  rotated  mode. 
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Figure  18b.  Glint  from  right  side  windows. 
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Figure  18d.  Glint  from  top  window  and  all  windows  combined 


ftftftftftftft 
kiftft  ft  m  ft  ft  ft" 

i «  A 

444*  4  4  M  4  *  V 

444*4*4* * 


V  1 

MW 

»HME 

4  ft  ft  ft  ft  ft  ft  ft'. 

\ 

■MU : 

fWITMT 

*\ 

HI*1!'!  it 

•wsMSM 

4  ft  ft  ft  ft  ft 

ft\ 

wwy* 

<-iaau; 

ft  ft  ft  ft  ft 

ft  *\ 

,-mmK A 

ft  ft  ft  ft 

m  ^  \ 

j  i 

J  .  <n#4»: 

*  ,  "  4:  a  *  i 

A**i»  A.-*  A  A  A  A'  A  A  A  A^sAUHl 

iiwAiui  >.  A  A  A  A  A  A  A  A  iA  *cC-W8K 

*  a  a  a  a  a  *  *  k  «-*« ir 

a  A  A  A  'A  A  *1  «  K  AA-riW 

Hm*Wkaa  a  a  a  .a  a  a  -a  a  k  <**•»»: 
HIW^A  AAA  hVa  A  A  •<  £  *»*»; 

<w-: 

*>: 

-io 


AA 


Figure  19a.  Glint  surface  points  for  side  window  2.0  inch  displacement,  top  window 
1.5  inch  displacement  in  rotated  mode. 
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Figure  19c.  Glint  from  left  side  windows. 


Figure  20a.  Glint  surface  points  for  side  window  3.0  inch  displacement,  top  window 
)  .5  inch  displacement  In  routed  mode. 
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Figure  20c.  Glint  from  left  side  windows, 
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Figure  21a.  Glint  surface  points  for  side  window  4.0  Inch  displacement,  top  window 
1.5  Inch  displacement  in  routed  mode. 
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The  elevation  and  azimuth  angles  of  the  particular  sun  positions  used 
in  this  study  are  listed  in  Table  1,  along  with  the  corresponding  plot 
symbol.  The  elevation  angle  is  measured  from  the  positive  z-coordinate  axis 
(vertical)  of  the  aircraft.  The  azimuth  angle  is  the  displacement  from  the 
positive  x-coordinate  axis  (pitch)  of  the  projection  into  the  x-y  plane  of 
the  sun- aircraft  line.  The  y- coordinate  axis  is  directed  along  the 
longitudinal  axis  of  the  aircraft. 


TABLE  1 

Elevation  and  Azimuth  Angles  (Degrees)  and 
the  Corresponding  Plot  Symbol  for  the 
Sun-Positions  Used  in  the  External  Glint  Study 


Case 

Symbol 

Elevation 

Azimuth 

1 

X 

5° 

0° 

2 

+ 

5° 

90° 

3 

V 

95° 

0° 

A 

A 

95° 

O 

O 

ON 

The  directions  of  the  reflected  rays  are  given  on  the  plots  by  their 
elevation  and  azimuth  angles.  The  elevation  is  plotted  on  the 
abscissa-coordinate,  while  the  azimuth  is  plotted  on  the  ordinate.  The 
elevation  is  measured  from  the  positive  z-coordinatt  axis  of  the  aircraft. 

A  ray  reflected  at  an  elevation  equal  to  zero  would  be  seen  by  an  observer 
above  the  aircraft.  A  ray  at  an  elevation  of  180  degrees  would  be  seen  by 
an  observer  directly  beneath  the  aircraft.  The  azimuth  is  the  angular 
displacement  from  the  positive  y-axis  of  the  projection  of  the  ray  onto  the 
x-y  plane.  A  ray  reflected  at  an  azimuth  angle  equal  to  zero  is  seen  by  an 
observer  behind  the  aircraft.  A  ray  at  an  azimuth  of  180  degrees  would  be 
seen  by  an  observer  in  front  of  the  aircraft. 
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Figures  15a  through  15d  show  the  external  glint  which  is  generated  by 
the  present  LGC  design  at  the  sun  positions  used  in  this  study.  Figure  15a 
shows  the  surface  points  used  to  compute  the  external  glint  on  side,  top, 
and  front  views  of  the  canopy  frame.  Figures  15b  and  15c  show  the  angular 
directions  of  the  reflected  glint  for  the  right  and  left  side  windows, 
while  Figure  1 5 d  shows  the  glint  for  the  top  window  and  all  windows 
combined. 

Figures  16  through  18  show  similar  plots  for  the  aircraft  with  the  top 
window  in  the  rotated  configuration.  The  figures  show  the  glint  results  for 
the  top  window  displaced  0.5,  1.0,  and  1.5  inches  from  the  vertice  plane, 
and  the  side  window  held  fixed  at  a  1.0-inch  displacement.  Figures  18 
through  21  show  plots  for  side  window  displacements  of  1.0,  2.0,  3.0,  and 
4.0  inches  from  the  vertice  plane,  and  the  top  window  held  fixed  at  a 
1.5-inches  displacement. 

The  figures  show  a  progressive  increase  in  the  angles  at  which  the 
glint  from  the  aircraft  is  visible,  as  the  displacements  of  the  side  and 
top  windows  are  increased.  However,  the  locus  of  the  glint  angles  is  a  line 
for  all  sun  positions,  and  the  glint  is  not  visible  over  a  solid  angle.  The 
line  locus  for  the  glint  angles  results  from  the  choice  of  cylindrical 
segments  for  the  LGC  design  side  and  top  windows.  The  locus  of  the  glint 
angles  for  the  remaining  flat  windows  are  single  points  for  any  one 
sun-position. 


Optical  Distortions 

The  results  of  the  computations  for  the  optical  distortions  and 
transmittance  are  shown  in  Figures  22  and  23.  The  figures  show  point-wise 
measures  of  the  optical  distortions  and  t ransm i t tances  spaced  at  equal 
intervals  on  side,  top  and  front  views  of  the  aircraft  canopy  frame.  The 
figures  for  a  particular  choice  of  side  and  top  window  displacements,  are 
presented  in  paired  sets  with  the  first  figure  showing  the  optical 
distortion  and  the  second  figure  showing  the  transmittance  values.  Figure 
22  shows  the  results  for  the  present  LGC  design,  while  Figure  23  shows  that 
for  the  design  with  a  4-inch  side  window  displacement  and  a  1.5-inch  top 
window  displacement  in  the  rotated  mode. 

A  comparison  of  Figures  22  and  23  shows  how  the  two  designs  differ  in 
optical  distortions  and  transmittance.  The  modified  design  of  Figure  23 
shows  an  increase  in  optical  distortion  and  a  decrease  in  transmittance 
over  that  portion  of  the  top  window  which  is  directly  over  the  pilot. 
However,  the  optical  distortions  and  transmi t tances  are  not  changed  over 
the  remainder  of  the  window  surfaces  from  that  for  the  present  LGC  design. 
This  is  true  to  the  degree  of  measure  exhibited  in  the  figures. 
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Figure  22a.  Optical  distortions  on  present  design. 
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Figure  23a.  Optical  distortions  for  side  displacements  4.0  inches  and  top  displacement  of  1 .5  inches 
in  rotated  mode. 
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DISCUSSION 


The  results  for  the  internal  glare  computations  are  summarized  in 
Figure  24.  The  figure  shows  the  total  possible  number  of  primary  order 
reflection  points  for  a  2-degree  by  2-degree  viewing  increment,  as  a 
function  of  the  side  and  top  window  displacements.  The  plotted  data  was 
computed  for  variations  in  the  design  of  Figure  14  with  the  top  window  in 
the  rotated  configuration.  The  figure  shows  that  the  number  of  reflection 
points  on  the  side  windows,  plotted  as  solid  lines,  decreases  as  the  side 
displacement  is  increased.  The  number  of  reflection  points  is  nearly  zero 
for  side  displacements  larger  than  four  inches.  This  is  true  for  a  2-degree 
by  2-degree  viewing  increment. 

The  figure  shows  also  the  number  of  reflection  points  on  the  Lop 
window,  plotted  as  dashed  lines.  There  is  a  slight  interaction  between  the 
top  and  side  windows,  and  the  number  of  top  reflections  is  plotted  as  a 
function  of  the  side  displacement  for  each  top  displacement  considered.  The 
number  of  top  reflection  points  decreases  as  the  top  displacement  is 
increased.  No  top  reflections  occur  for  displacements  equal  to  or  larger 
than  1.5  inches. 

The  secondary  and  triplex  order  reflection  points  were  computed  using 
the  sane  program.  Figure  25  shows  tile  total  possible  number  of  secondary 
reflection  points  (2-degree  by  2-degree  increment)  as  a  function  of  the 
side  and  top  displacements.  The  number  of  secondary  reflections  is  much 
less  than  that  for  the  primary  order.  The  side  window  reflections  vanish 
for  side  displacements  which  are  greater  than  2  incites.  The  number  of  lop 
reflections  are  a  function  of  both  the  top  and  side  displacement s .  The  top 
window  reflections  vanish  for  top  displacements  equal  to  or  greater  than 
1.5  inches.  The  number  of  possible  triplex  reflections  are  trivally  small 
and  are  not  covered  further  in  this  report. 

A  study  of  Figures  14a  and  1 4 b  shows  that  a  few  reflection  points  art- 
possible  on  the  upper,  forward  corners  of  the  rear  side  windows  and  the 
edge  of  the  upper  front  window.  Figure  14b  shows  that  the  reflections  are 
caused  by  rays  entering  at  the  lower,  front  cover  of  the  front  side 
windows.  The  few  reflection  points  which  are  possible  may  be  blocked  during 
nighttime  flight  by  the  addition  of  visors  over  the  reflection  surfaces 
area  or  at  the  entry  positions  for  the  reflected  rays. 

In  summary,  one  modification  to  the  LGC  design  which  shows  little 
internal  glare  is  that  with  a  4-inch  side  window  displacement  and  a 
1.5-inch  top  window  displacement  in  the  rotated  mode.  Figure  21  show  that 
although  the  external  glint  angles  are  increased  by  use  of  the  modified 
design,  the  locus  of  the  glint  angles  for  a  given  sun  position  forms  a  line 
due  to  the  use  of  the  cylindrical  side  and  top  windows.  Figure  21  shows 
that  the  use  of  the  modified  design  increases  the  amount  of  optical 
distortions.  However,  the  distortion  is  limited  to  that  portion  of  the  top 
window  which  is  directly  over  the  pilot's  position. 
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Figure  25.  Secondary  order  reflection  points  as  a  function  of  side  and  top  window  displacements. 


RECOMMENDATION S  EOR  FURTHER  RESEARCH 


The  following  Is  recommended  for  further  research: 

1.  Test  the  modified  LGC  design  with  a  4-inch  side  window 
displacement  and  a  1.5-inch  top  window  displacement  in  the  rotated  mode,  on 
a  mockup  of  the  Model  YAH-64  Advanced  Attack  Helicopter  available  at  USAHEE 
for  such  purposes. 

2.  Investigate  combinations  of  other  window  shapes  with  zero 
Gaussian  curvatures  (i.e.,  curvature  in  only  one  direction)  along  with  the 
cylindrical  shape.  The  locus  of  external  glint  angles  from  any  such 
combination  will  form  a  line  on  a  plot  of  elevation  and  azimuth  for  a  given 
sun  direction. 

3.  Expand  the  analysis  techniques  described  in  the  report  to  the 
computation  of  internal  glare  reflections  from  the  instrument  panel  decals 
and  lights.  This  is  easily  accomplished  by  modifying  the  attached  computer 
program,  DRWR ,  to  consider  reflected  rays  which  originate  from  the 
instrument  panel  area. 

4.  Experimentally  determine  the  correlation  between  the  measure 
of  optical  distortions  developed  for  this  report,  and  the  visual 
performance  of  the  pilot.  In  particular,  the  limiting  value  of  optical 
distortion  beyond  which  pilot-vision  is  seriously  hindered  should  be 
determined. 


5.  Develop  further  the  relationship  between  the  external  glint 
and  the  tactical  threat  by  enemy  forces.  In  particular,  the  use  of  zero 
Gaussian  curvature  surfaces  (such  as  cylindrically  shaped  windows)  results 
in  a  glint  pattern  which  has  no  solid  angle.  The  objective  would  be  a 
performance  criteria  to  weight  the  glint  pattern  that  is  produced  by  a 
canopy  design. 


6.  Investigate  further  the  use  of  mathematical  optimization 
techniques  to  design  an  optimum  configuration  by  computer  for  the 
transparent  surfaces  of  the  canopy.  The  design  would  minimize  both  glint 
and  glare  in  accordance  with  some  performance  measure,  while  maintaining 
few  optical  distortions. 


CONCLUSION 

A  computer  study  of  the  low  glare  canopy  (LGC)  design  for  the  Model 
YAH-64  Advanced  Attack  Helicopter  has  been  completed  by  the  USAHEL.  The 
study  employed  computer  programs  developed  by  USAHEL  to  compute  point -wise 
measures  of  three  important  factors  determining  canopy  performance.  The 
factors  are  (1)  the  internal  glare,  (2)  the  external  glint  and  (3)  the 
optical  distortions  exhibited  by  the  design.  The  results  of  the  study  are 
determined  that  certain  modifications  to  the  present  LGC  design  will  reduce 


the  internal  glare  with  only  slight  increases  in  external  glint  and  optical 
distortions. 

The  modifications  to  the  present  LGC  design  suggested  for  further 
study  are  (1)  increase  the  side  window  displacement  from  the  plane  of  the 
frame  vertices  to  4  inches,  (2)  rotate  the  top  window  90  degrees  and  (3) 
increase  the  top  window  displacement  to  1.5  inches.  The  modified  design 
increases  the  angular  span  over  which  external  glint  may  appear.  However, 
the  glint  does  not  appear  in  a  solid  angle  and  is  judged  not  to  be  a 
tactical  disadvantage.  The  optical  distortions  are  increased  on  the  top 
window  in  the  region  directly  above  the  pilot.  However,  the  distortions  are 
not  increased  on  the  remainder  of  the  windows  at  least  to  a  degree 
noticeable  by  the  method  of  measure  used  in  this  report. 
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The  equations  used  in  the  computations  are  discussed  in  this  section. 
The  elevation  and  azimuth  angles  of  the  sun-position  and  the  reflected 
direction  are  measured  from  the  cartesian  coordinate  system  of  the 
aircraft.  The  y-axis  is  directed  along  the  longitudinal  axis  of  the 
aircraft  from  nose  to  aft  and  the  coordinate  values  follow  the 
station-marks.  The  z-axis  is  directed  vertically  and  the  coordinate  values 
are  the  level-marks.  The  x-axls  is  orthogonal  to  the  other  two  axes  and 
directed  toward  the  left-hand  side  of  the  aircraft  (facing  forward)  for  a 
right-handed  coordinate  system.  In  the  following,  we  consider  the 
specification  of  the  sun-position  and  the  computations  of  the  specular 
reflection  from  the  canopy  surface  and  the  reflected  direction. 

a.  Specification  of  the  Sun-Position. 

The  sun-position  is  specified  by  the  angles  of  elevation  (a) 
and  azimuth  (p).  The  elevation  is  the  angle  between  the  positive  z-axis  of 
the  aircraft  and  the  straight-line  from  the  aircraft  to  the  sun.  The 
azimuth  is  the  angle  between  the  positive  y-axis  and  the  projection  of  the 
aircraft-sun  line  into  the  x-y  plane.  The  directional  cosines  (as,bs,cs) 
of  the  incident  sunlight  are: 

as  =  sin  q  slnp  (1) 

bs  =  sin  a  cosp 
cs  =  -cosa 

b.  Specular  Reflection  of  the  Sunlight  From  the  Canopy  Surface. 

The  angles  that  the  specular  reflection  and  the  incident 
sunlight  are  to  the  normal  of  the  canopy  surface  are  equal  in  value  and 
opposite  in  sign.  The  reflection  lies  in  the  plane  defined  by  the  sunlight 
and  the  surface  normal  (Appendix  A,  ref  4).  Furthermore,  the  angle  between 
the  sunlight  and  the  normal  of  the  exterior  surface  of  the  canopy  must  be 
acute  (less  than  90  degrees)  for  reflection  to  occur.  The  dot  product  (Q) 
between  the  directional  cosines  of  the  sunlight  (as,bs,cs)  anc*  surface 
normal  (an,bn,cn);  i.e., 

Q  ■  asan  +  bsbn  +  cscn» 


must  be  greater  than  zero,  Q>0.  Here,  the  surface  normal  is  directed  into 
the  cockpit.  The  directional  cosines  of  the  reflected  ray 
(an>bn,cn)  is  given  by  (Appendix  A,  ref  2). 


as~2anQ  > 


(3) 


ar  = 


br  -  bs~2bnQ, 
cr  =  cs-2cnQ. 


The  same  equations  apply  to  the  computations  for  cylindrical 
surfaces.  The  directional  cosines  of  the  surface  normal 
(an,bn,cn)  is  a  function  of  the  surface  position  (x,y,z)  where  the  particular 
ray  of  sunlight  reaches  the  surface.  It  is  a  function  also  of  the  radius 
of  the  cylinder  (rc),  and  the  origin  (x0,y0,z0)  and  directional  cosines 
(ao»bo*co)  of  the  cylinder-axis  (Appendix  A,  ref  3). 


c.  Angular  Coordinates  of  the  Reflected  Ray 

The  elevation  (ar)  and  azimuth  (Pr)  of  the  reflected  ray 
are  computed  from  the  direction  cosines  (ar,br,cr)  of  equations  (3), 
and  the  inverse  of  equations  (1);  i,e., 

ar  =  Acos  ( c r )  (4) 

Pr  =  Asin  (ar/sin  ar) 
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APPENDIX  B 


OPTICAL  DISTORTIONS 


The  computation  of  point-wise  measures  of  optical  distortions  and 
t  ransmittances  are  briefly  discussed  in  this  section.  The  program  first 
determines  a  grid  of  equally  spaced  surface  points  on  the  windows  of  the 
canopy  by  the  same  method  described  in  Appendix  A  for  external  glint.  The 
program  then  computes  the  directional  cosines  of  the  straight  line  ray  from 
the  pilot's  nominal  eye  position  to  the  surface  point  using  the  coordinates 
of  the  line  end  positions. 

The  transmittance  and  lateral  magnification  is  computed  by  tracing  the 
ray  through  the  window  to  the  outside.  Consider  a  ray  from  the  pilot's  eye 
position  to  the  window  surface  point  (xs,ys,zs).  The  ray  is  refracted  by 
its  passage  through  the  window  surface  into  the  material.  The  directional 
cosines  of  the  refracted  ray  (ar,br,cr)  are  determined  as  a  function  of 
those  of  the  incident  ray  (as,bg,cs)  and  the  surface  normal  (an,bn,cn)  by 
Snell's  Law, 


ar  =  as/n  +  an1s » 
br  =  t>s/n  bnqs, 
cr  =  cs/n  +  cn<ls- 

The  parameter,  n,  is  the  index  of  refraction  of  the  window  material.  The 
factor  qs  is  given  by: 


qs 

=>  _I 

Qs  “ 

Qr ,  where 

n 

Qs 

-  ~( 

as* a 

n  +  bs**n  +  cs*cn>> 

Qr 

=  _1 

(n2 

-  1  +  Q2) 1/2. 

n 

The  transmittance  (Tj)  of  the  refracted  ray  at  the  internal 
* jrface  is  computed  from  the  angles  of  incidence  (9y)  and  refraction 


-sing  techniques  already 

1  v«n  bv  ? 

described  (Appendix  A,  ref  4). 

The  angles  are 

9  Y  “  Acos  (9g)» 

and 

(2) 

0r  •  Acos  ( 0r)» 
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The  refracted  ray  is  next  traced  through  the  window  material  to  the 
outer  surface  using  the  surface  point  as  the  starting  point  of  the 
refracted  ray.  The  equations  for  the  intercept  of  a  known  ray  with  a  planar 
surface  are  described  elsewhere  (Appendix  A,  ref  4).  The  outer  surface  of 
the  window  is  assumed  to  have  the  same  surface  normal 
(an»bn»Cn),  but  the  edge  vertex  (xgl , yg’ , ze* )  is  displaced  outward  by  the 
window  thickness  (t) 

The  coordinates  of  the  edge  vertex  of  the  outer  window  surface  are 
given  by: 


xe1 

=  Ke^nS 

(3) 

ye1 

ye~bn  , 

ze1 

=  ze_cnt* 

in  terms  of  the  inner  surface  edge  vertex,  the  surface  nomal  and  the  window 
thickness.  Substituting  equations  (3)  into  the  appropriate  equations 
(Appendix  A,  ref  4)  gives  the  straight  line  distance  (R)  between  the 
internal  surface  point  and  the  intersection  of  refracted  ray  with  the 
external  surface, 

R  =  (an(xg-xe)  +  bn(ys-ye)  +  cn(zs-ze)  +  t)/Qr.  (4) 


The  equations  for  the  interception  point  of  a  ray  specified  by  a  known 
origin  and  directional  cosines  with  a  cylindrical  surface  are  described 
elsewhere  (Appendix  A,  ref  5).  The  origin  and  directional  cosines  of  the 
cylindrical  axis  for  the  outer  surface  is  assumed  to  be  the  same  as  that 
for  the  inner  surface.  However,  the  cylindrical  radius  (rc^)  is 
assumed  to  be  increased  by  the  window  thickness  (t);  i.e.,  rcl  =  rc  +  t.  The 
substitution  of  the  outer  surface  radius  into  the  appropriate  equations 
gives  the  straight  line  distance  along  the  refracted  ray  from  the  inner 
surface  point  to  the  outer  surface  point. 

The  intersection  point  (xo,yo,zo)  of  the  refracted  ray 
with  the  outer  window  surface  is  given  in  terms  of  the  straight  line 
distance , 


Xo  =  XS  +  arR,  (5) 

yo  =  ys  +  brR, 

Zo  =  Zs  +  crR. 

The  ray  transmittance  through  the  material  (Tm)  is  determined  by 
the  material  coefficient  of  absorption  (ara)  as: 

Tm  =  EXP(-am.R).  (b) 
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The  transmitted  ray  Is  refracted  at  the  external  window  surface.  The 
directional  cosines  (a0,bo,co)  of  the  refracted  ray  are 
determined  by  Snell's  Law, 

ao  -  nar  +  anq0,  (7) 

bo  =  nbr  +  bnq0 , 
x0  =  ncr  +  cnq0. 

The  directional  cosines  (an,bn,cn)  are  of  the  surface  normal  at 
the  point  (x0,yo,zo).  The  factor  qQ  is  given  by: 

q0  =  nQr 1  -  Qo,  where 

Qr*  *  -(ar*an  +  br'bn  +  cr«cn),  and 


Qo  =  ( 1+n2  (1-Q1r2)  ) 1/2 

The  transmittance  (T0)  of  the  refracted  ray  at  the  outer  window 
surface  is  computed  from  the  angles  of  incidence  ( 0  }  )  and  refraction  (90), 
where:  1 


9  i"  Acos  and 

®0  =  Acos  (Qo)* 


\ 

) 


However,  the  transmitted  ray  is  refracted  only  if  the  factor  Q0  is 
less  than  or  equal  to  unity;  i.e.,  Q0  <1.  In  this  case,  the  angle  of 
incidence  is  less  than  or  equal  to  the  critical  angle  and  the  transmittance 
is  greater  than  zero. 


The  material  transmittance  (T)  for  the  passage  of  the  ray  through  the 
material  is  given  by, 


T  -  Ty  •Tm»T0.  (9) 

The  lateral  magnification  (Mo)  for  the  external  ray  compared  to 
the  internal  ray  is  given  by  (ref  2,  p. 152) , 

Mo  ”  Acos  (ao*as  +  bo’bg  +  cQ*cs). 


The  optical  distortion  is  a  measure  of  the  change  of  the  lateral 
magnification  over  the  effective  field  of  view  (ref  2,  p.152).  The  optical 
distortion  may  be  computed  for  each  window  surface  point  by  using  the  above 
ray  as  a  principal  ray.  We  arbitrarily  consider  a  4-degree  cone  of  vision 
centered  at  the  pilot's  position  and  about  the  principal  ray. 
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A  family  of  subsidary  rays  is  constructed  for  the  cone  in  such  a 
manner  that  the  solid  angles  contained  by  adjacent  rays  are  equal.  The 
lateral  magnification  (Mj)  is  computed  for  each  one  of  the  subsidary 
rays  by  the  ray  tracing  method  outlined  above.  The  optical  distortion 
(mj)  is  arbitrarily  defined  as  the  square-root  of  the  sura  of  the 
squares  of  the  differences  between  the  magnification  for  the  principal  ray 
(M0)  and  each  subsidary  ray  (Mj),  divided  by  the  number  of  rays  in  the  cone. 
The  equation  is: 


md 


N 

£ 

t=i 


(Mo-Mi) 2/ N 


1/2 


(11). 


APPENDIX  C 


COMPUTER  PROGRAMS 


The  computer  programs  are  programmed  for  the  CDC  7500  computer 
batch-job  system  in  Fortran  IV  language  (1,3).  The  subroutines  which  are 
common  to  all  three  programs  are  listed  first.  The  subroutines  unique  to 
each  program  follow  along  with  the  program  listing. 
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SUBROUTINES  COMMON  TO  ALL  THREE  PROGRAMS 


The  subroutines  common  to  all  three  programs  are  those  used  to  read  in 
frame  data  and  compute  the  corresponding  surface  normals,  and  those 
controlling  the  computer  graphics  output.  The  subroutine  are: 

1.  READV  -  reads  in  canopy  frame  data  into  common  storage  from 
frame  data  file.  Called  by  main  program. 

2.  NORML  -  computes  surface  normal  for  planar  surface  segments  of 
canopy  and  plane  of  frame  vertices  for  cylindrical  surfaces.  Called  by  main 
program. 


3.  DRAW  -  subroutine  establishes  drawing  file  for  graphic  picture 
of  side,  top,  and  front  views  of  canopy  frame.  Called  by  COMPQ  in  programs 
DRWC  and  DRWE ,  and  by  DRV  in  DRWR. 

A.  EMARK  and  ETIC  -  subroutines  compute  pilot's  eye  positions  in 
drawing  file.  Called  by  DRAW. 

5.  LINS 
LINF 

LINT  -  subroutines  compute  frame  position's  in  draw-i ng- - f-Ll-e-. 
Called  by  DRAW. 

6.  BELT 

PLEM  -  subroutines  establish  drawing  file  in  form  visible  by 
Calcomp  drawing  routines.  Called  by  DRAW  and  routines  4  and  3  above. 

7.  DRW  -  subroutine  calls  on  Calcomp  routines  (see  item  9  below) 
for  drawing  canopy  side,  top,  and  front  views  using  established  drawing 
file.  Called  by  COMPQ  in  programs  DRWG  and  DRWE  and  by  DRV  in  DRWR. 

8.  DRWP  -  subroutine  calls  on  Calcomp  routines  (item  9)  for 
drawing  surface  points  on  canopy  frame  drawing.  Called  by  COMPQ  in  programs 
DRWG  and  DRWE,  and  by  DRV  in  DRWR. 

9.  PLTBEG 
PLTSCA 
PLTSYM 
PLTDTS 

PLTPGE  -  subroutines  for  processing  Calcomp  drawing  (ref  3). 


INTERNAL  GLARE 

The  subroutines  peculiar  to  this  program  are  listed  below.  The  list 
does  not  include  those  routines  described  elsewhere  (4,5).  A  copy  of  the 
program  is  attached. 
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1.  DRWR  -  main  program  controlling  read  in  of  canopy  frame  data, 
computation  of  surface  normal  and  Indexing  process  for  computing  Internal 
reflections. 


2.  DRV  -  subroutine  reads  in  parameters  for  incrementing  radius 
of  cylindrical  windows,  calls  for  computation  of  internal  glare  and 
controls  computer  graphics  output.  Called  by  DRWR. 


EXTERNAL  GLINT 

The  subroutines  are  listed  below  followed  by  the  program. 

1.  DRWG  -  main  program  controlling  read  in  of  canopy  frame  data, 
computation  of  surface  normals  and  indexing  processes  for  reflection 
computations. 


2.  DRV  -  subroutine  reads  in  parameters  for  incrementing  radius 
of  cylindrical  windows,  specifies  sun  angles  for  which  reflections  are  to 
be  computed,  and  calls  for  computation  of  reflected  directions  as  a 
function  of  sun  position  for  each  cylinder  radius  considered.  Called  by 
DRWG. 

3.  COMPQ  -  subroutine  calls  for  computa.t  ion....flf  xei  l.e-c_t.e^_. 
d  TrectTons  ~  a  s  a  Tu  nctl  ono  f  sun  angle  and  drawings  of  canopy  frame  with 
reflection  points  and  plots  of  the  reflection  directions  for  each  cylinder 
radius  considered.  Called  by  DRV. 

4.  DRVO  -  subroutine  established  grid  of  equally  spaced  surface 
points  on  canopy  surface  and  computes  reflection  angles  for  the  specified 
sun  position.  Called  by  COMPQ. 

5.  INTCY  -  subroutine  computes  intersection  point  of  straight 
line  with  cylindrical  surface.  Called  by  DRVO  to  establish  grid  of  surface 
points . 


6.  TRSCY  -  subroutine  converts  coordinates  of  a  cylindrical 
coordinate  system  into  a  rectangular  system.  Called  by  INTCY  when  testing 
an  intersection  point  of  a  line  with  a  cylindrical  surface  against  the 
frame  edges  of  the  surface. 


OPTICAL  DISTORTIONS 

The  subroutines  are  listed  below  followed  by  the  program. 

1.  DRWE  -  main  program  controlling  the  read  in  of  canopy  frame 
data,  computation  of  surface  normals  and  indexing  processes  for 
computations. 
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2.  DRV  -  subroutine  reads  in  parameters  for  incrementing  radius 
of  cylindrical  windows  and  calls  for  computations  for  each  cylindrical 
radius  considered.  Called  by  DRWE. 

3.  COMPQ  -  subroutine  calls  for  the  computation  of  the  optical 
measure  and  the  drawings  showing  the  distribution  of  the  optical  measures 
on  the  canopy  frame  for  the  specified  cylindrical  radius.  Called  by  DRV. 

4.  CANL  -  subroutine  establishes  grid  of  equally  spaced  surface 
points  on  canopy  surface  and  calls  for  computation  of  optical  measures  for 
the  angular  direction  to  each  surface  point  in  turn.  Called  by  COMPQ. 

5.  CALC  -  subroutine  determines  the  transmittance  and  optical 
distortions  for  the  specified  angular  direction.  Called  by  CANL. 

6.  INTEC  -  subroutine  determines  the  intersection  point  of  a 
straight  line  ray,  from  the  pilot's  normal  eye  position,  with  a  planar 
surface  and  whether  or  not  the  point  lies  within  the  planar  segment.  Called 
by  CALC. 

7.  IN1CY  -  subroutine  computes  intersection  point  of  straight 
line  with  cylindrical  surface  and  determines  whether  point  is  enclosed 
within  the  cylindrical  segment.  Called  by  CALC. 


8.  TRSCY  -  subroutine  converts  coordinates  of  a  cylindrical 
coordinate  system  into  a  rectangular  system.  Called  by  INTCY  when  testing 
an  intersection  point  of  a  line  with  a  cylindrical  surface  against  the 
frame  edges  of  the  surface. 

9.  WKTEC  -  subroutine  computes  the  exit  position  and  directional 
cosines  of  a  ray  passing  through  a  planar  surface,  along  with  the  optical 
transmittance  for  the  ray.  Called  by  CALC. 

10.  WKTCY  -  subroutine  computes  the  exit  position  and  directional 
cosines  of  a  ray  passing  through  a  cylindrical  surface,  along  with  the 
optical  transmittance.  Called  by  CALC. 

11.  TRANS  -  computes  the  transmittance  for  a  ray  passing  through 
a  dielectric  surface  from  Fresnel's  Law.  Called  by  WKTEC  and  WKTCY. 
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SMYTH, STMfZ*T1000. 

ACCOUNT »  HE*** , 

BEGIN, ATTACH,PLOTLIB. 

REQUEST, TAPE13,*PF. 

FTN<SL,P ) . 

MAP, ON. 

I  GO  • 

BEGIN, PLOT, CALC OMP,TAPE13. 

EXIT. 

BEGIN, PLOT, CALC ONP.TAPE13. 

EXIT. 

PROGRAM  DR WR( INPUT, OUTPUT, TAPEA,TAPE13,TAPE2-INPUT1 
C  MAINLINE,  RAY  TRACING  FOR  PILOT 

C  3  DIMENSIONAL  FLAT/CYLINDRICAL  CANOPY  KITH  INSTRUCTIONS 
COMMCN/LB/STRG(<t) 

DATA  STRG/"INTtRNAL", "GLARE", "REFLtCTl ON", »>"/ 

CALL  READV 
CALL  NORML 
CALL  ORV 
STOP 
ENO 

SUBROUTINE  DRV 

COMM ON /CAN/NT, NB,NC,NP,NA,ND,NV(10C),PXV(100, 8), PYV(  100, b),PZV(  IOC 
C  ,  8  ) 

COMMON/C YL/NC Y, N3C (10),NSPI10,10),XC(10),YC(10),ZC(10),AE(10),Bc(l 

00)*CE(10),RC(10>  _  _ _ _ 

COM  M  ON  y.C.QJNV-.T.R^JSLj  A  X  1-00-r  1-S) - - 

COMMON/ LR/8,RCC,B1,RCC1 
DIMENSION  A W ( 9 ) , A R ( A ) 

DATA  AWt9),AR(4)/2*lH>/ 

nc  o*nc ♦  i 

R£AD(2,d0C) 

READ(2,801)XH»XP,YP,ZP»AN,BN,CN,YA,7A 
READ(2,801)XH1,XP1,YP1,ZP1,AN1,BN1,CNI,YAI»ZAI 
READI2, 800) 

READI2, 802)80, BM,DS,NB 
RE  AD ( 2,  800) 

READ(2,802)B01,BM1,0B1,NB1 

REA0(2,800) 

R  E  A  D (  2,  80  3  )  KR 

800  FORMAT ( ) 

801  F0PMAT(2X,A(r-I0.A,2X)/fct2X,F10.4)> 

802  F0RMAT(?X,3(F10.A,2X),I2) 

803  FORMAT ( 2X, 12  ) 

OJ  ?C  IC«1»KR 
KRR«IC+1 

I  R  «  1C 
B»  ttO-DB 
DO  10  I  *  1 , NB 
B  *  B ♦ DB 

RCC«.5*(Xh**2+B»*2) /B 
DO  3  IJ-1,2 
RC  t  1J)»RCC 

XCUJ)»-(XP-(RCC-B)*AN)*((-1)**IJ) 

YC ( IJ)*YP— IRCC-B) *  B  N  ♦  Y  A 
ZC ( IJ I -ZP-<  RCC-B ) *C N  +  Z  A 
3  CONTINUE 
B 1 ■ BO 1-D8  1 
DO  1C  K *  1 , N Bl 
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Bl-eitDBl 

RCCl-.5*(XHl**2tBl**2)/Bl 
R C ( 3  I -RCC 1 

XC ( 3  >-XPl-(RCCl-Bl » *  AN  1 
YCt  3  )-rPl-(RCCl -81  )*8N1*1TA1 
zcm-zPi-<RCci-ei>*cNi+ZAi 
NVLL-0 

00  5  IK-1,KRR 
DO  5  I P-NCO, ND 

5  NSR ( IP,  IK  )-0 
CALL  ORAM 
CALL  DR W ( 13) 

REWINO  <. 

CALL  CANUNVLL) 

CALL  ORWP(NVLL) 

ENC0DE(73»1000»AW)B»RCC»8I,RCC1 

1000  FQRMAT(2X,"SIDE  WINDOW", 2X,2(F10. 4, 2X), "TCP  W 1 N 00-", 2 ( 2 X , F  1 0  .  A  )  ) 
CALL  PLTSYM( .12b, Aw.0.,600., 5.  ) 

ENCODE (23, 1001, AR)KR 

1001  F0RMAU2X, "TOTAL  R  E  FL  E  C  T  I  ON  S  »»  2  X ,  I  2  ) 

CALL  PLTS YM( . 125, AS ,  0.,600. ,-7.  ) 

ENC0DE(79,1002,AW)(IP,IP-NC0,ND) 

1002  F0RMAT(2X,«SURFACE",10(4X,I3») 

CALL  dLTSYM(.125,AW,0.,600.,-19.) 

YCK-O. 

DO  6  IK-1, KRR 

ENCODE ( 79, 1003, AW ) (NSR ( IP, IK )» IP-NCO,ND ) 

1003  FGPMAT(2X,  "POINTS"/  IX,  10(1XU6>)  _  .  . . . 

L  AXIL  PT  T  S  YK (  .12 57 A W » 6  .  »  6 00  .7  -  31.-Y C  'K  ) . 

YCK-YCK+12. 

6  CONTINUE 

ENC0DE(22,100A, AR)NVLL 
10  0  A  F0RPAK2X, "TOTAL  PU  I N  TS"  ,  2  X  ,  16  ) 

CALL  PLTSVM( .125, AR,0.»b00.,-31.-YCK ) 

CALL  PLTPGE 

10  CONTINUE 
20  CONTINUE 

RETURN 

END 

SUBROUTINE  CANUNVLL) 

C  CONTROLS  CALCULATION  OF  REFLECTION  POINTS 

C0MhCN/CAN/NT,N6,NC,NP,NA,ND,NV( 100) »PXV(100»  b) >PYV(100»fc),PZV(100 

0,6  ) 

DATA  All, AI2/2.,2./ 

C  INDEX  PILOT  VIEWING  DIRECTION 

C  A 1  ELEVATION,  ANGLE  FROM  Z  AXIS  1CWARU  X  AXIS  IN  X  Z  PLANE 
C  A2  AZIMUTH,  ANGLE  FROM  X  Z  PLANE  TOWARD  Y  AXIS 

C  Z  AXIS  TOWARD  UPWAKO*  X  AXIS  TOWARD  LEFT  FACING  FRONT  OF  CANOPY  AND  Y 
C  TOWARD  BACK  OF  CANOPY  ALONG  LONGITUDINAL  AXIS  OF  AIRCRAFT 
A  I  -0 . 

1C  A2-AI2 

11  A2-A2-AI2 


I F (A2.GE.-9C.  ) 

GO 

TO 

15 

Al-Al+AIl 

IFU1.GT.18C.  ) 

GO 

TO 

299 

GO  TO  10 
15  CONTINUE 
I  S  -  0 

CALL  CALC ( Al, A2, IS.NVLL) 
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IF  (  IS.EQ.NAICALL  CALC(A1,A2,NA,NVLL) 

GO  TO  11 
299  CONTINUE 
RETURN 
ENO 

SUBROUTINE  CALCU1,A2,NK,NVLL) 

C  PRIMARY  REFLECTION  ONLY,  SURFACE  NK  TRANSMITTER  ONLY 

COM  NON/C AN/NT, NB»NC,NP,N A, ND»NV( 100) ,  PX V< 100*d),PYV(100»F),PZV<K 
3,8  ) 

CJMMON/PILOT/XO,YO, ZO 
CC1MMCN/LINE/AS,BS,CS,XS,  YS.ZS,  AC,BC,CC 
CUM MON/ CON /KM, N5R I  100,10  1 
COMMON/ FACT /SX, AX0,BXC» A  XI, BX1 

DIMENSION  IPP<10>»ASP(10>»BSP(10>,C$P(10),XP(10),YP<10)»ZP<10),AI 
C(10),8IP(10),CIP(10),ARP(101»BRP(10),CRP(10),RP(10) 

DATA  PI/3.14159/ 

A1M«A1*PI/180. 

A2M«A2*PI /180. 

AS«CCS ( A2M) *s  IN( A  1M  1 
BS-SINU2M) 

CS«C03(A2M)*CCS ( AIM  1 

R«l. 

X  S  «  X  C 
Y  S  ■  YO 
ZS«Z( 

I P  U  N  ■  0 

'IF  CONTINUE 

DO  2C  I S* 1 »  ND 

IF ( IS .EG. INK  JGOTU  20 

I  SK-O 

I F (  IS.LE.NA1CALL  INTEC(1SK,IS,X*,YR,ZR) 

IF ( 1S.GT .NA1CALL  INTCY(1SK,IS,XR,YR,ZP) 

IF ( ISK ,GT .0 1G0TG  25 
2 C  CONTINUE 
RETURN 
25  CONTINUE 

IF (  IS.LE.NC1RETURN 
RX«AS*AC+BS*8C*CS*CC 
IF (Rx.GT.C. IGCTU  20 
CALL  COMP ( ANG,R T, TT1 
T*  TT*« 

R • R  T*R 

IF l IS.EQ.NK1G0T0  40 
I P  UN  *  I  PUN  +  1 
A  I  *  —  A  S 
a  I  .-PS 
Cl «-C5 

AP«-AS+2. *RX* AC 

BR»-BS+2.*RX+bC 

C  R  ■— C  S+-2.*RX*CC 

IF  (  IPUN. EG.  (KM4-1  IIGOTC  30 

IF (R.LT. .  00001 1  RETURN 

1PP (  I  PUN  1 ■ TS 

AS  P (  I  PUN  I  * AS 

BSP ( IPUN  )  *BS 

CSP  ( IPL'N  1  »CS 

XP(  I  PUN  1  *  XP 

YP<  IPUN  1  *  YR 
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ZP ( IPUN ) ■ ZR 
A  I P  <  IPUN I • 4  I 
BIP< IPUNI-BI 
C I P ( IPUN ) «C  I 
ARP( IPUN)»AR 
BRP( IPUN ) *BR 
CR  P ( IPUN ) «CR 
RP<  IPUN)  .R 
AS  «-AR 
BS«-BR 
CS--CR 
X  S  “  X  R 
YS-YR 
ZS«7R 
GO  T C  10 
3G  CONTINUE 

IF  U .LT . .COOOl ) RETURN 
NVLL  «NVLL  +  1 
OG  3  5  I  K ■  1  #  K H 

wRITE(4,1000)XP(IK),YP(Ih),ZP(I*) 
1C0C  FORNAT « 4< F10.4)  ) 

IPS-IPP(IK) 

NS«(IPS.IK)»NSRUPS*IK)  +  1 
3 5  CONTINUE 

—  ;»?  I  TL  •  4#1  000  R  *  ZRf  T  _ 

NiR(lS»KM*U"NSRll$.KN  +  l)  +  l 
NK-IPPI1) 

RETUPN 
40  CONTINUE 

1 N  K ■ NK 
GO  TC  18 
END 

SUBROUTINE  DRWPINVLL) 

CUMNON/CON/KM#NSM 100# 1G) 

C0NN0N/FACT/SX>A0*B0»A1,B1 

REWIND  4 

DO  10  I ■ 1  *  N VL  L 

DO  5  I K ■ 1 # K M 

RE  AD( 4* 1CCCIXP, YP, ZP 

XX-(YP-5.5)*SX*10C. 

YY-(ZP-111.52)*SX+100. 

CALL  °LTDTS<3/0»XX.YY,1*0) 

5  CONTINUE 

REA0L4. 1C0C)XR,YR»ZR»T 
IC— AL0G10(T)*1. 
XX.(YR-S.5)*SX»1C0. 
YY«(ZR-111.52)*SXfl00. 

CALL  PLT0TS(3»IC»XX,YY.1»0) 

10  CONTINUE 
1000  F0RNAT(4(F10.4)  ) 

REWIND  4 
DO  20  I«1#NVLL 
DO  IS  IK-l.KM 
READL4, 1CC0)XP,YP,  ZP 
XX-IYP-5. 5)*SX*1CC. 

YY«(XP*104. )*SX  +  433  .44 
CALL  PLTDTS(3»0»XX»YY#1#0) 

15  CONTINUE 

RE  AD( 4» 10C0)XR» Yf .ZR#T 


IC--AL0G10LTH1. 

XX-<YR-S.5)*SX+100. 

YY-LXR+104.  HSX  +  433.44 
CALL  P  L  TD  T  S ( 3»IC» XX.YY.1.0) 

20  CONTINUE 
REWINO  4 
DO  30  I  •  l »  NVL  L 
00  2 t  I  K  »  1 »  K  P 
READ  1 4.  1 000 ) X  P  »  YP  *  Z  P 
XX-(ZP-111.52)*SXF690.42 
YY«tXP+104.|*SX+433.44 
CALL  PLTDTSL3#0,XX,YY»1.0> 

2S  CONTINUE 

R  E  A  C  L  4  » 10  00)XR»YP»Z.R»T 
IC--AL0G10L T  > ♦ 1 . 

XX»(ZR-lll.i2)*SX+690.A2 
YY-LXR+10 A.  HSX  +  433.44 
CALL  »LT0TS13#IC.XX,YY,1,0> 

3C  COfTINUE 
RETURN 
FND 

SUtiRTUTINE  READV 
C  R  t  A  P  IN  SURFACE  VERTICES 

CUPP  I'N/CANANT.NE.NC.NP.Ni.NQ. NVL ICO) »  P  X  V  (  1 0  0  #  c  )  »  P  Y  V  (  10  0. f 1 •PZYtlOC 

0*81  .  _ .  -  -  - 

'  CQPPON/ VEPT/XVL 200 )* YVL200) . ZV(200) ,NVR ( 100, 8 ) 

COPPLN/CYL/NCY.NSCL 10)»NSP( 10* 10)*  XC (lOJfYC (10) «ZC( 10)*AE( 1 0  )  »  8  E (  1 
C  0  >  *  C  E ( 10) *  R  C I  10) 

CUPPCN/PILCT/XU.YO.ZO 
RE  ADL  2. 1000) 

1000  FOPPATO 

RE  A0(2.  1000  ) 

R£ACi{2»10Cl)NI»N8»NC»NP»NA*ND 

1001  FQRPATLM2X.I3)) 

RE  AE ( 2,  1000) 

KE  ACL  2*  1002  ) LnV ( 1 ) ,  I*1»ND ) 

1002  PORPATL  10  (2X,  13  )  ) 

8  E  A  C ( 2*  1000  ) 

00  1C  J  *  1 . 8 

READ12»10Q2) ( NVk ( I »  J ) » I » 1 , NO ) 

10  CONTINUE 

RE  A0( 2» 1000 ) 

RE A0L2*  1003  I ( XV<  I  ) *  I  -  1  *  N  T  ) 

1003  F  OR  P A  T ( 8 ( 2X.F7.4 )  ) 

R  E  A  D ( 2*  1000) 

RF  A0(  2,  1003  > l YV (1) .  1-1  .NT  ) 

RE A0L2. 1000) 

R  E  A  0 ( 2 .  1003  M  Z V (I) *  1- 1»N  T  ) 

RE  AO ( 2.  1000  ) 

Rt  AT  (  2»  1 0  02  )  N C  Y 
READL2.1002HNSCm.I-l.NCY) 

00  20  I-l.NCY 
k P -NSC L  I ) 

20  REAPL2.10C2)  LNSPL  I.K),K-1.kP) 

READ(2»10  0AMXCII)»YCH),ZC(I). ALLI), 3E<I).CEtI>»RCll), 1*1,  NCY) 

1004  FORMATL  7<  2X.F7. 4 ) ) 

READL2. 1000) 

REAf L2.1C03)xn. YO.ZO 
R  E  T  U  E  N 
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END 

SUBROUTINE  NOP  ML 

C  ESTABLISH  SURFfcCE  NORMAL  FOR  EACH  PLATE  SURFACE 
C  SURFACE  NORMALS  DIRECTED  TOnARO  COCKPIT  INTERIOR 

COMMON/ CAN/ NT# NB# NC»NP#NA#ND»NV< IOO) »  PXV( 100# S >#FYV(lC0.fc)#PZv(10Q 
0*8) 

COMMON/ VE  RT/XVI 200).YV(200)»ZV(200)»NVR(1CO#9) 

COMMCN/NORM/AXN ( 100)#  AYN{ 100  )» AZN( 100) 

CQMM0N/CYL/NCY»NSCU0)#NSPI10»10)#XCI10)#YC(10)»ZC(10>»AE«10).BE11 

C0l»Cc(10)#RC(10) 

COMMON/PILOT/XO#  YO#  ZO 
DO  1C  I  ■  1 »  ND 
NK»NV( I ) 

DO  S  K  »  1  #  NK 
KV-NVRI  I  »  K  ) 

PXVd.K  )  -  X  V  (  K  V  ) 

P  Y V ( I.K).YV(KV) 

PZ V( I .K ) »  ZV(KV) 

5  CONTINUE 
K  ■  c 

7  A  1  *  P  X  V ( I»K)-PXV(1»  1  ) 

A2«PXV<  I»K  +  1)-PXV(  J#1  ) 

Bl-KY V< I.K ) -P  Y  V  l  I# 1 ) 

P2»PY V ( I*  K*1 )-PYV ( I . 1 ) 

C1*PZV(I#K)-PZV(I#1) 

C2«PZV(  I»K+1)-PZV(  1 . 1  ) 

Pl«SCRT(Al**2+bl**2+Cl**2) 

P  2  »  S  CR  T (A2**2+BE*A2+C2**2) 

A«(Al*A2  +  Bl*e2«-Cl*C2)/(Pl*P2) 

IFIABS(A).LT.l.)  GO  TO  9 
■  K  »  K  ♦  1 

IF  (K.EO.NK)  GO  TO  10 
GO  10  7 
9  A  N  »  A  C  0  S (A  ) 

R«  1  .  /  ( P  i*P2*S  1  N ( AN  )  ) 

A XN  (  1  )  ( B1*C2-C1*B2 ) *R 

A  Y  N ( I  )•-( A1*C2-C1*A2)*P 
AZN (I) ■♦ ( A1*B2-A2*B  1 )  *P 
1C  CONTINUE 
PETORN 
END 

SUBROUTINE  IMEC(IiK,IS#*R#YR»ZR) 

C  DETERMINES  IF  RAY  STRIKES  CONVEX  SURFACE 

COMMON /C AN /NT#  NB#nC»NP#N A, ND#NV ( 10C)#PXV(1CO#H).BYV(100»B)»PZVI100 
0#  8  ) 

C0MMCN/NJRM/AXN(10C)#AYN( 100)>AZN( 100) 
CQMMON/LINE/Aj#BS#CS#Xi,YS»ZS»AC»BC#CC 
CK»AXN(IS(*AS+AYN(IS)*BSaAZN(IS)*CS 
IFICK.GE.O.)  RETURN 

C  RAY  STRIKES  SURFACE  IN  OUTUAkO  UlRtCTICN 
I«1 

IF(XS.E3.PXV(IS»I).ANO.YS.EO.PYV(IS»I).ANC.ZS.E/.PZV<IS»I))I*I+1 
S»(AXN(IS)*iPXV(IS#I)-XS)*AYMlS)*(FYV(IS#I)-YS)  +  AZN(IS)*(PZV(lS#I 
C  >-ZS> )/CK 
>R»AS*b+*S 
YR«8S*S+YS 
ZR-CS*S*ZS 
A1«PXV(IS.1)-XS 
B1  - P VV (  IS# 1  )-YS 
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HUMAN  EN6INEERIN6  LAB  ABERDEEN  PROVING  6R0UND  MD  F/G  1/3 

A  STUDY  OF  THE  CANOPY  DESIGN  FOR  THE  ADVANCED  ATTACK  HELItOPTER— ETC(U) 
AUG  79  C  C  SMYTH 

HEL-TM-13-79  NL 


Cl-PZVI IS#1  )-ZS 
Pl-XR-XS 
P2-YP-YS 
P3-2P-ZS 
IN-NVI IS) 

OQ  1C  I  ■  1  »  I  N 
1C-I+1 

IFd.50.lN)  IC-l 
A2«PXVdS»IC)-XS 
8 2-FYVI  I S  » I C ) - YS 

C2 •  p z v ( is#io-zs 

Q»P1*(B1*C2-B2*C1)-P2*IA1*C2-C1*A2)+P3'Ma1*6Z-B1*A2) 

IP (O.GT .0. ) RE  TURN 

C  RAY  STRIKES  SURFACE  UN  ENCLOSED  S10E  Of  SURFACt  cDGt 
A 1 «  A  ? 

61  »B? 

C1«C? 

10  CONT iNUc 

C  RAY  'TRIKE*-  ENCLOSED  SURFACE 
I  S  K  ■  1 

A  C  *  A  XN (  IS  ) 

FC- AYN<  IS  ) 

CC-AZNl  IS  ) 

RE  TURN 
ENO 

SUP  ROUTINE  INTCYIISK.Is.Xk.YR.ZR) 

C  COMPUTES  INTERSECTION  POINT  OF  LINE  WITH  CYL1NDE* 

COMMCN/CAN/NT»NB>NC,NP#NA#ND#NVdOC) »  PXV<1  0C#6)  ,»YV<  ICO#  8  )  »Ri  VdCC 

0,8) 

COMMON/ LINE /AS, BS#CS»XS#  YS#  ZS» AC. 6C»CC 

C0MFLN/C»L/NCY,NSC(10),NSP(1U,10),XC(10)#YC(10)»ZCI10)#AE(10),BEC1 
OC) >CE ( 10) #RC ( 10) 

CPMMGN/PILOT/XP.YP.ZE 

C  CtTEFMINE  CYLINCERICAL  SURFACE  WHICH  CONVEX  SUkFACE  IS  A  HART  GP 
C'O  2  I  ■  1  #  NC  Y 
KP-NSC  (  I  ) 

DC  2  K  *  1 #  K P 

IF  (  N  S  P  ( l  *  K  )  ,  El> .  1  S  )  GOTO  5 
2  CONTINUE 
RETURN 
5  CONTINUE 

C  OF  Tr RHINE  INTERSECTIUN  POINT 
XO-XC ( I ) 

YO • YC (  I  ) 
iO«ZC<  I  ) 

RO-RC (  I  ) 

AC  * Af  (  I  ) 

BO -BE  (  I  I 
CO-Ct  1 1  ) 

RC3*SORT((XS-XO)**2+(YS-YO)**2+(ZS-ZG)**2) 

ACO* ( XS-XC ) /PCS 
POO- ( YS-YO)  /ROS 
CCO*  tZS-ZO) /POS 
A-AO*AS+BC*BS+CD*CS 
B-KCS*UO*AOO  +  eO*BOO  +  CO*COO) 

Al« as-a*ag 
A2-FS-A+B0 
A3-CS-A+C0 
B1»ACO+ROS-B*aO 
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B2*BCO*ROS-B*BO 
B3»C00*R0S-B*C0 
A.Al**2tA2**2tA3*»2 
B- Al*BltA2*B2»A3*B3 
C»Bl^*2tB2**2+B3**2 
AB  —  B/ A 

AX«B**2-A*(C-RO**2> 

IF(AX.LT.0.>AX-0. 

BB-SORTUXJ/A 

BS-AB-BB 

IFUS.LE.O.  »RS«Ab*BB 
IF  (PS.LE.O. IRETURN 
XR-XS+AS+RS 
YR  » Y  5  »BS*  RS 
2R«Z5+CS*RS 

R«RS*(AO*AS*bO*bS+CO*CS)+ROS*<AO*AOO+BO*bOO+CO*COO) 

Xl»XC+AO*R 

Yl-YO+BO*R 

ZWC  +  CO*R 

AC* IX1-XR)/R0 

bC"( Yl-YR  )  / R  0 

CC* ( 21-ZR  I  /RO 

C  DETERMINES  WHETHER  RAY  STRIKES  ENCLOSED  SURFACE 
XSS-XP 
YSS-YP 
ZSS-ZP 

CALL  TRSCYI I  ,  X  S  S  ,  Y  $  S  ,  Z  S  S  ) 
xRR«XR 
YR  R ■ YR 
Z RR.ZR 

CALL  TRSCY( I , XRR, YRR, Z RR  ) 

Pl-XPR-ASS 

P2-YPR-VS5 

P3*ZRR-ZSS 

X1«PXV(1S»1> 

Y 1 • P  Y  V ( IS»1 ) 

Z1«PZV< 15,1) 

CALL  T  R  SC  Y { I, XI, Y1,Z1 ) 

A  1 ■ X 1 “X  S  S 
91» Yl-YSS 
C1-Z3-ZSS 
I  N  ■  N  V  ( IS  ) 

03  10  1 1 ■  1 ,  I N 
IC-IIfl 

IF ( 1 1.E3. INJIC-I 
X2"PXV(IS,IC) 

Y2 ■ P  YV ( IS,  ICI 
Z2*  PZV ( I S  »  IC  ) 

CALL  TRSCY(I,X2,Y2,Z2) 

A2-X2-XSS 
B2-Y2-YSS 
C2-Z2-Z SS 

0-Pl*(81»C2-b2ACU-P2*(Ai*C2-Cl*A2)+P3*Ul*ii2-Pl+t2) 
IF(Q.GT.O. (RETURN 

C  RAY  STRIKES  SURFACE  ON  ENCLOSED  SIDE  OF  SURFACE  EDGE 
A1-A2 
P 1  ■  B2 
C  1  »C2 

1C  CONTINUE 
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C  KAY  STB  IK  t  S  ENCLOSED  SUBFACE 
ISK-1 
RETURN 
END 

SUBROUTINE  TRSCY<I»XR,YR,ZR> 

C  CONVERTS  COaRDlNATES  IN  CLINDRICAL  COORDINATES  INTO  RuCT  ANGULAR  iPACt 
COPKON/CYL/NCY,NSCnO)*NSP(  10,  10),XC  <l0)*YCll0)*ZC(10>,At<l0),BE<i 
QC) .CE(IO) »RC( 10) 

CN-SCRT ( l.-CE (I  )  **2  ) 

AN  * -A  e ( I ) *CE (I)  /CN 
BN«-BE(II*CE(n/CN 
AP-BE  <  I ) *CN-BN*CE (  I  ) 

3  P  ■  -  l  A  E (I )*CN-AN*CE ( I )  ) 

CP-AE  < I  )  *BN-AN*BE (  I  > 

BO»AE(I)*(XR-XC(Il)+BEm*(YR-YCm)+CE(I)*(ZR-ZC(I)) 

XX-XCU  )+AE(I)*RO 
YY«YC(I  ) ♦ BE  ( I  )  *R  0 
ZZ-ZC ( I )+CE (I)*RO 

R«SaPTl(XP-XX)**2MYR-YY)**2  +  tZR-ZZ)**2) 

AA* (XR-XX  )/R 
BB« ( YR-YY  )  / R 
CC«!ZR-ZZ)/R 
A«AN*AA+BN*BB+CN*CC 
IFU.GT.l. ) A  ■  1 . 

If  (  A  .LT.-l  .  I  A  —  l. 

ANG*ACOS( A  ) 

G»AP*AA+BP+B8+CP*CC 
IF (C.LT .0. I ANG»-ANG 
XR*ANG*RC  (  I  ) 

YR«PC 

ZR«F-RC ( I ) 

RETURN 

END 

SUBROUTINE  C OMP ( A NG , R T , T T  ) 

C  COMPUTES  INCIDENCE  ANGLE*  f>  t  F  LE  C  T  A  N'C  E  *  AND  T  R  An  UN  I  T  T  A  NC  E 
C  NATURAL  LIGHT*  ADDITION  OF  POLARIZATION  COMPONENTS  IGNORED 
C  XN,  INDEX  OF  REFRACTION,  TX,  INTEkNAL  TRANSMTTANCE 
COMPUN/LINE/AS,BS,CS»XS»rS»ZS,AC,BC.CC 
DATA  XN,TXX/1.5» .05/ 

TT«0. 

RT«1. 

A»-AS* AC-BS*BC-CS*CC 
IFIA.LE.O.IRETURN 
IF ( A.GT. 1 .  IA-1. 

ANG*ACOS( A) 

AN  C  P  *  AS  1  N(  SIN  (  ANG  ) /XN  ) 

CA-COS( ANG) 

SA«S IN( ANG) 

SI *SCRT ( XN**2-S A**2  ) 

RO»mCA-Sl)/(CA  +  Sl))**2+((CA*(XN*  +  2)-il)/(CA*(xN**2)  +  Sl))**2Jx2. 
T0«(1.-R0)*CA/C0S(ANGP) 

CA-CrS( ANGP) 

SA-SINI ANGP) 

TX«EXP(-TXX/CA) 

S 1  *  SQR  T  <XN**2-SA**2) 

RI«(((CA-S1)/(CA  +  S1))**2-*-((CA+(XN**2)-S1)/(CA*(xn**2)  +  S1))**2)/2. 
TI-(1.-RI)*CA/CQS(ANG) 

TT»T0*TI*TX/(1.-(RI*TX>**2) 

RT«RO+RI*TT 
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RETURN 

END 

SUBROUTINE  DRAW 

C  CRAWS  GRAPHIC  PICTURE  OF  CANOPY  IN  THREE  FOLD  LAYOUT 

COMMON/ DR A/MX(10C5»NS<100)*X$(  1000 )  #  Y  S ( 1000 )  .  NS Y  ( 1 00 ) # NE  »  NX 
COMHON/CAN/NT»NB»NC*NP*NA,NO*NV<100)»PXV<100*8>»PYV< 100# 6 ) # P Z V <1 00 
0*8) 

C0MM0N/N0RH/AXN(100)»AYN< 100 ) » AZN ( 100 ) 

COMMON/PILOT/XO. YO*  ZC 
COMMON/ F ACT /$X»  AC* BO# Al* B1 
NX  ■  0 

SX-2.122 
NS  S -NC ♦ 1 
AO  —  5.5 
B0-1C0. 

Al  —  111.52 
Bl-100. 

XX«(Y0+A0)*SX+B0 
Y  Y- (ZOtAl )*SX+B1 
CALL  BELT(I*XX» YY, 1,0) 

CALL  EM  ARM  XX»YY) 

NE  *  1 

DO  2  I-NSS*ND 
QS— AXNtl  ) 

IfUI.LE.NC.AND.CS.LT.O.I.OR.U.GT.NC.AND.QS.GT.O.))  GOT  U  2 
C  SURFACE  FACES  VIEWER  FROM  SIDE  VIEW 
XX-(PYV(I»l)+AG)*SXiBO 
YY-(PZV(I»1)4A1)*SX+B1 
NE -NE*1 
MD-1 

IF ( I.LE.NC.OR.I.EQ.NAJMD-4 
CALL  BELT (NE»XX, YY»MD*0) 

CALL  L INS ( I  I 

2  CONTINUE 
A  1*104. 

Bl-433.44 
XX»(YO+AO)*SX+BC 
YY- (X0+A1  )  *  S  X  +  B 1 
NE -NE+1 

CALL  BELT (NE»XX»YY, 1,0) 

CALI  EMARM  XX*  YY  ) 

DO  3  I-NSS.ND 
QT-4AZNI I  ) 

IF((I.LE.NC.AND.OT.LT.O.).OR.(I.GT.NC.AND.QT.GT.O.)IGOrO  3 
C  SURFACE  FACES  VIEWER  FROM  TU*  VIEW 
XX-(PYV(I#1)+A0)*SX4B0 
YY- (PXV( I  *  1 ) + A 1 )*SX+B1 
NE -NE+1 
MO-1 

IF( I<LE<NC<OR .1 .E0.NA)MD-4 
CALL  BELT (NE»XX#YY,MD*C) 

CALL  LINT  1 1 ) 

3  CONTINUE 
A0  —  1U.B2 
BO ■ 600 .42 
XX-(Z0+AC)*SX+80 
YY- (XO  +  Al  )*SX  +  B1 
NE -NE ♦ 1 

CALL  BFLT(NE*XX,YY,1,0» 
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CALL  ETIC (XX.  YY  ) 

DU  10  I-NSS.ND 
OF  «-A  YN ( I  ) 

IF((1.LE.NC.AND.0F.LT.0.).0R.(I.GT.NC.AND.QF.GT.0.))GCTU  10 
L  SURFACE  FACES  VIEWER  FROM  FPONT  VIEW 
XX-(PZV(I.1J*A0I*SX+B0 
YY* (PXV( I  .  1 )  +  A 1  ) *S  X ♦ B  1 
NE  -NE  ♦  1 
HD-1 

IF ( I.LE.NC.OR.I  .EQ.NA1H0-4 
CALL  BELT (NE.XX, YY.HD.O) 

CALL  LINF(I) 

10  continue 

NS  <  NE  +  l ) «  N  X ♦ 1 

RETURN 

END 

SUBROUTINE  BELT ( IE*  X.Y.MOD.NSS  ) 

CONHUN/DRA/MO(100)*NS(lCO)»XS(luCC)»YS(iOCJ)»NSY(100)»NE#EP 

NP-NP+1 

MS  1 T f  > -NP 

NSY (  IE  )  *NSS 

HD ( IE ) -MOD 

XS ( NP ) -X 

YS (NF  )  * Y 

RETURN 

END 

SUBROUTINE  P  L  E  M  (  X  »  Y  ) 

CGHHGN/DR  A/MD (100).NS(100)#XS(1000)»Y5(1000)»N5Y(100)»NE.NP 

NP-NP+ 1 

XS (NF  )-X 

Y  S ( NP ) =  Y 

RETURN 

END 

SUhRCUTINt  EMARK(X.Y) 

C  HARK  EYE  POSITION  IN  CANOPY 
X-X-10. 

Y-Y-5. 

CALL  PLEM(X.Y) 

X-X+10. 

Y - Y  +  5 . 

CALL  PLEH(X.Y) 

X-X-10. 

Y-Y+5. 

CALL  PLEH(X.Y) 

X  *  X  ♦  3  • 

Y-Y-1.5 

CALL  PLEM(X.Y) 

Y-Y+3. 

CALL  ? L  F M ( X  *  Y  ) 

Y- Y-14. 

CALL  PLE  H ( X  *  Y  ) 

RETURN 

END 

SUBROUTINE  ETIC(X.Y) 

Y-Y-5. 

CALL  PLEM(X.Y) 

Y-Y+10. 

CALL  PLEH(X.Y) 

Y-Y-5. 
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CALL  PLEM(X»Y) 

X-X-5. 

CALL  PLEM(X.Y) 

X-X+10. 

CALL  PLENUM) 

RETURN 

END 

SUBROUTINE  LINSII) 

COMhON/CAN/NT#NB*NC#NP»NA»ND»NV{ 1001 #PXV  < 100#  b ) » P V V  < 100#  M # KZ V < 1C C 

0.8) 

CGMMON/FACT/SX.AO»BO*A1#B1 
NK-NVl I ) 

DO  10  K* 1 »  NK 
K 1 »K  +1 

IF (K.E0.NKIK1-1 
XX«(P'YV<I,Kl)tAO)*SX*BC 
YY-(PZV(1*K1)+A1)*SX+81 
CALL  PLEM(XX#YY) 

10  CONTINUE 
RETURN 
F  N  C 

SUBROUTINE  L I NF (  I  > 

CGf'MON/CAN/NT#Nb*NC  #NP,NA.ND»NV(  ICG  >»PXV  (10^*0), <>YV<100»i.)*PZy(1U0 
Q  *  P  ) 

COMMON/FACT/SX, A0#B0»A1, il 
NK  »NV( I  ) 

DO  10  K • 1 »  NK 
K  1  *  K  ♦  1 

IF (K.EO.NK)Kl-l 
XX»(PZV(I#K1>+A0)*SX+PD 
YY«  l  PXV<  I  #K  1)  +A1  )  ♦S*  +  t)l 
CALL  PLEM(XX,YY) 

10  CONTINUE 
RETURN 
END 

SUBROUTINE  L I  NT ( I  ) 

C0NMCN/CAN/NT»NB#NC#NP,KA#ND»N\M100 )*PXV( 100»8) #PYV< 10C»8l #PZV( IOC 
0*8) 

COHMON/f ACT/SX* A*#BC#Ai*Bl 
NK»NV( I ) 

DO  10  K ■ 1  *  NK 
K  1  ■  K  ♦  1 

IF(K.E0.NK)K1«1 
XX«<PYV<I#Kl)*AC)*iX+BO 
YY»<PXV( 1#K1)+A1 ) *  S  X*  B 1 
CALL  PLEM(XX.YY) 

10  CONTINUE 
RETURN 
ENO 

SUBROUTINE  DRW(IL'NIT) 

CONMON/DRA/MD(10C)»N$nOO)*XS(  1000)»YS(  1000>#NSY(10C)»N£*NP 
COMrON/LB/STPGU) 

C0MM0N/DG/XP(1CC)#YP(100) 

DIMENSION  L  ABEL ( 4  ) 

DATA  LABEL /"SMYTH"#  "6520".  "X365A"#  "HE  379"/ 

C  PRCGRAr  COMPUTES  LAYOUT  ON  102A  UMTS  FITTED  TJ  25  TNChES 
CBY  25  INCHES  DISPLAY  SHEET 

CALL  PLTBEG(25.#25.#l.#lUNlT#LA8EL) 

CALL  PLTSCA(5.#5.#2  00.*0.*5C.<)6,A0.Ofe) 
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CUL  PLTSYM(.25#STRG#0.*201.»1.» 

DO  10  I E ■ 1»  N£ 
f'DD  »MD  ( I  E  I 
NSS-NSY(  IE  1 
N1-NS(IE> 

N2«NS( IE*1»-N1 
DO  5  I ■ 1 »  N2 
XP( I ) ■ X  S ( N1 *1-1 ) 

Y  P I  1 ) ■ Y  S ( N1 ♦ 1-1 ) 

5  CONTINUE 

CALL  PLTDTS(MD0.NSS»XP>YP»N2»0) 

10  CONTINUE 
RETURN 
END 

C  A  AH  CANCPY  DATA 

C  NUM3ER  CF  VERTICES.  5 Uk F AC E S — C ONS T R A  I N T .  ELAT.  ANC  CYLINDRICAL 


166 

lb 

56 

56 

59 

6  A 

NO. 

VERTIC 

c  S  1 

PER  SURFACE 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

£ 

5 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

7 

7 

7 

7 

A 

A 

6 

A 

A 

A 

A 

A 

A 

A 

5 

A 

A 

5 

7 

7 

7 

7 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

6 

ft 

ft 

6 

6 

A 

VFRTICE  ASSIGNED  TO 

EACH 

SURFACE 

089 

09  C 

066 

121 

125 

129 

133 

1  35 

139 

1 A  3 

14  6 

1  A9 

153 

157 

137 

1  A5 

0C1 

002 

OC  3 

CCA 

007 

001 

005 

012 

Oil 

021 

015 

028 

022 

032 

029 

033 

110 

035 

036 

111 

038 

039 

03b 

109 

0  A  0 

0  A  1 

112 

053 

CA7 

060 

0  5  A 

061 

061 

089 

C  9  A 

093 

100 

099 

C  9  8 

087 

077 

081 

101 

068 

071 

09A 

ICO 

035 

117 

1  1  A 

115 

122 

126 

130 

13A 

136 

1  A  0 

1  AA 

1 A  ft 

150 

1  5  A 

158 

1  36 

1A8 

002 

003 

COA 

005 

0  1  A 

007 

OOA 

01  A 

012 

020 

016 

027 

02  3 

031 

030 

0  3  A 

109 

03ft 

037 

112 

039 

038 

1  1  1 

111 

0  A  1 

11C 

110 

052 

0A6 

059 

055 

06  A 

Oft2 

090 

089 

09A 

0  °5 

100 

099 

098 

078 

082 

102 

069 

0  72 

093 

099 

086 

1  1A 

115 

119 

123 

127 

131 

130 

137 

l  A  1 

1  AG 

1  A7 

151 

166 

159 

162 

163 

009 

010 

Oil 

012 

012 

005 

003 

008 

009 

019 

017 

026 

02  A 

C  3  0 

031 

109 

035 

OA  3 

OA  A 

044 

C  A6 

034 

109 

037 

045 

112 

0  A  2 

051 

0A9 

05B 

0  5  6 

063 

063 

066 

0  65 

C7C 

071 

076 

075 

07A 

079 

083 

103 

070 

073 

092 

096 

087 

090 

066 

065 

1 2  A 

12b 

132 

129 

138 

1  A  2 

13° 

1  Ae 

152 

155 

160 

161 

1  6  A 

006 

OC  9 

CIO 

Oil 

005 

002 

002 

009 

010 

018 

Olfi 

025 

025 

G2  9 

032 

11C 

0A2 

0  A  2 

0  A  3 

037 

0  A  5 

033 

0  3  A 

036 

0  Aft 

0  A  5 

0A3 

050 

C50 

057 

057 

062 

06  A 

065 

0  70 

069 

076 

075 

07A 

086 

080 

084 

1 0  A 

065 

0  7  A 

091 

097 

068 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

165 

156 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

017 

019 

02A 

026 

0 

0 

0  A  1 

0 

0 

0 

0 

0 

0 

0 

035 

0 

0 

0  A  A 

0A9 

051 

056 

05  e 

0 

G 

0 

99 


0 

0 

0  0 

0  0 

0 

0 

105 

C66 

075  C9C 

096  0 

0 

0 

0  0 

0  0 

0 

0 

0 

0 

0 

0 

0  0 

0  0 

0 

0 

0 

0 

0 

0 

0  0 

0  016 

020 

023  027 

0 

C 

40 

0  0 

0  0 

0 

0 

0 

0 

0 

0 

0  048 

052  055 

059 

0 

c 

C 

0 

0 

0  0 

0  0 

0 

0 

106 

067 

076  089 

095  0 

0 

0 

0  0 

0  0 

0 

0 

0 

0 

0 

0 

0  0 

C  0 

0 

0 

0 

C 

0 

0 

0  0 

0  015 

021 

022  028 

0 

0 

0 

0  0 

0  0 

0 

0 

0 

0 

C 

0 

0  047 

C  53  054 

060 

0 

0 

0 

0 

0 

0  0 

0  0 

0 

0 

107 

0 

0 

0 

0  0 

0 

0 

0  0 

C  0 

0 

0 

0 

0 

0 

0 

0  0 

c  c 

0 

0 

0 

0 

0 

0 

0  0 

0  0 

0 

0 

0 

0 

0 

c 

0  0 

0  0 

0 

0 

0 

0 

0 

0 

3  0 

0  0 

0 

0 

0 

0 

0 

3 

0  C 

0  0 

0 

0 

108 

c 

0 

0 

0  0 

vti  net 

X-8TS IT  ION 

2.00 

2.00 

2.00 

2  •  OC 

2. CO 

2 .  G  0 

2.00 

-2  .OC 

-2.0  0 

-2  .00 

-2.00 

-2.00 

-2.00 

-2.00 

19.06 

18.70 

14.25 

12.90 

11.2C 

14.76 

17.42 

-le.oo 

-17.60 

-13.84 

-13.30 

-13.24 

-15.96 

-17.96 

8.50 

-11.50 

- 12 • v4 

9.ce 

19.63 

19.63 

11.46 

11.46 

11.46 

19.63 

19.63 

-19.63 

-19.63 

-11.46 

-11.46 

-11.46 

-19.63 

-19.63 

19.16 

18.84 

14.22 

13.10 

11.24 

14.76 

17.50 

-17.90 

-17.60 

-14.00 

-13.30 

-13.28 

-15.96 

-17.96 

6.56 

-11. 4« 

-12.14 

c.  lb 

23.08 

20.92 

18.08 

13.72 

22.72 

23.23 

23. 2C 

13.60 

13.00 

15.16 

16.36 

22 . 84 

16.72 

-16.72 

-14.28 

14.28. 

13.60 

-13.60 

-9.  20 

9.20 

10.56 

-1C. 56 

-11.56 

11.56 

-23.08 

-20.92 

-16.08 

-13.72 

-22.72 

-23.28 

-23.20 

-13.60 

-13.00 

-15.16 

-16.36 

-22.84 

19.63 

13.26 

11.46 

-11.46 

-13.26 

-19.63 

-12.00 

12.00 

13.26 

-11.46 

13.26 

-11.46 

-23.75 

-4. 

4. 

23.75 

-4. 

-23.75 

4. 

2  3.75 

24. 

36. 

36. 

24. 

-36. 

-24. 

-24. 

-36. 

26. 

104  . 

104. 

26. 

26. 

104. 

67.6 

61. 

61. 

67.6 

-104. 

-2b. 

-26. 

-104. 

-104. 

-26. 

-61  . 

-67.6 

-67.6 

-61. 

-46. 

-24. 

-22. 

-56. 

24. 

46. 

56. 

22. 

-10. 

10. 

10. 

-1C. 

61. 

61 . 

-61. 

-61. 

-56. 

56. 

VEtTICE 

Y-P0S1T1UN 

57.50 

57.50 

66.37 

71.57 

53.98 

58.98 

60.30 

57.50 

57.50 

68.37 

71.57 

58.98 

58.98 

60. 30 

84.66 

82.92 

83.46 

86.98 

98.18 

98 . 66 

96.24 

84.66 

t2.9  2 

83.46 

86.98 

98.18 

98,66 

96.24 

87.48 

87. 4r 

97.56 

97.56 

98.45 

103.49 

100.92 

110.92 

115.61 

121.40 

121.50 

98.45 

103.49 

100.92 

110.92 

115.61 

121.40 

121.50 

144.26 

142.52 

142.54 

146.00 

157.52 

158.40 

156.32 

144.26 

142.52 

142.54 

146.00 

157.52 

158.40 

156.32 

148.00 

148.00 

156.66 

156.66 

59.20 

59.20 

7C.  14 

108.42 

110.90 

69.78 

116.49 

113.25 

139.69 

156.97 

156.97 

140.09 

58.30 

58 .3C 

67.16 

67.16 

69.57 

69.57 

108.01 

108.01 

112.13 

112.13 

145.44 

145.44 

59.20 

59.20 

70.14 

108.42 

110.90 

69.78 

116.49 

113.25 

139.69 

156.97 

156.97 

140.09 

115.52 

115.52 

114.82 

114.82 

100 


115*52 

115.52 

113.01 

113.01 

103.69 

103.69 

121.60 

121.60 

57.5 

5.5 

5.5 

57.5 

5.5 

57.5 

5.5 

5  7.5 

59.2 

65.2 

132.61 

158.6 

65 . 2 

59.2 

156.6 

i32 .61 

186.61 

191.61 

226. 61 

236.61 

186.61 

191.61 

169.61 

189.61 

180.61 

180.61 

191.61 

186.61 

236.61 

226.6  1 

191.61 

186.61 

1 1  9 , 6  1 

189.61 

180.61 

180.61 

216.61 

216.91 

216.61 

216.61 

215.61 

216.61 

216.61 

216.61 

198.61 

198.91 

198.61 

19e.61 

198.61 

226.61 

196.6 

226.61 

216.61 

216.01 

c 

VEFTICF  l- 

POSITION 

1  3  1 . 8  1 

138.66 

166.27 

161.77 

133.97 

133.97 

132.83 

131  .el 

136.66 

166.27 

161. 77 

133.97 

133.97 

132.8  3 

126.60 

1  3C . 36 

165.26 

168.56 

168.08 

136.56 

127.96 

126.60 

130.36 

1  6  5 . 2  o 

168.56 

168.08 

136.56 

127.96 

118.61 

118.61 

156.00 

15  6. 0  C 

129.20 

168.00 

157.58 

160.51 

160. se 

139.20 

129. 20 

129. 2C 

16«. OC 

157.58 

160.51 

160.88 

139.20 

129.20 

167.56 

169.60 

1  0  6  ■  2  6 

167.66 

167.66 

156.20 

167.56 

167.5b 

169. 6C 

166,26 

167.66 

lb7 .68 

156.20 

167.56 

137. ffc 

137.-6 

175. oF 

1  7  5 . 6  u 

129.05 

161.61 

15*..  PC 

172.60 

15C.6b 

130.65 

166.36 

176.^9 

1  79.26 

179.26 

176.60 

150.80 

163.52 

163.5? 

166.60 

156.62 

156.72 

156.72 

175.15 

175.15 

177.93 

177.93 

182.12 

1 62  .  12 

129. C5 

161.61 

156.80 

172.60 

150.  tt 

1  3C  .  H  5 

1  6  0 . 3  c 

17  6.5-. 

175.26 

179.26 

176.00 

150 . oO 

152.56 

152.56 

159.03 

1 5  9 . 0  ? 

152.56 

152.56 

175.79 

175.79 

168. OC 

168.00 

1  39. 2C 

139.20 

166.2 

129.  1 

129.1 

166.2 

122.1 

127.25 

122.1 

1  2  7  .  C  5 

126.12 

126.12 

126.  12 

126.12 

126  .  12 

126.12 

126.12 

126.  1^ 

162.52 

160.52 

160.52 

162.52 

132.1c 

132.92 

132.52 

13?.5? 

111.52 

111.52 

160.52 

162.52 

162.52 

160.52 

132.92 

132.1? 

132.52 

132.52 

111. 52 

111.52 

167.92 

167.92 

176.52 

178,5? 

167.92 

167.92 

178.52 

178.52 

161.52 

181.52 

221. 52 

221.52 

111.52 

132.52 

111.52 

132.52 

157.92 

157.92 

c 

cylindrical  data 

o 

c  2 

60  61 

62  63 

66 

1 

-31 .9956 

135.89 

167.6713 

•  V* 

.  9bC5 

.196,0 

55.991o 

21.9956 

135.89 

167.6713 

.  0 

.  9805 

.i960 

55.9916 

0. 

133.6172 

139. 3729 

.  0 

-.9837 

-.1796 

62.9375 

c 

PILOT  EYE 

POS  IT  ION 

-1.5 

162.33 

171.2 

c 

CANOPY  EDOES  FITTtD 

TO  CYLINDRICAL 

SIDES 

2C.  763 

17.875  130 

.  91 

161.175 

.9592 

.0 

• 

2828 

6.98 

•  ' 

9995 

11.25 

0. 

126 

.175 

180.135 

0. 

-.1796 

9637 

0. 

0. 

c 

CANOPY  displacement 

Of  CYLlNOKlCAL 

ilCtS  FROM 

PLAT  w 

INOCW 

1 . 

5. 

1 . 

1 

C  TOP  WINDOW 

.5  2.5  .5  1 

C  ORDER  OF  REFLECTIONS 
1 
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smyth, stmfz»tiooo. 

ACCUUNT#HE***. 

d6GIN,ATTACH*PL0TLItt. 

REGUEaT, TAPE13#*PF. 

FTN(SL.R). 

NAP#  ON • 

LGO. 

BEGIN#  PLaT#CALCOHP# TAPE13. 
exit. 

bcGlN#  PlOT»CALCOMP»TAPE13. 

EXIT. 

PROGRAM  DR*G< INPUT# OUT  PUT, 1  APEO»TAPE13#TAPE2* INPUT) 

C  MAINLINE,  GLINT  ACCEPTANCE  ANGLES 

C  THREc  0  I  ME  NS  I  UAL  ELAT /CYLINDRICAL  CANOPY  INCREMENTAL 
CJMMUN/LB/STRGm 

DATA  ST RG /"EXTERN AL"#"GLINT"#"REFLECTI ON", ">"/ 

CALL  RE  AD V 
CALL  NURML 
CALL  o RV 
STOP 
tNU 

SuokuUT I Ne  DRV 

C0MM0N/CAN/NT»NB#NC»NP,NA,N0,NyI  100  )#PXV(100»8),PYVI  100,0), 

UP  ZV ( i0u#8  ) 

V.UMM0N/  SUN /K1NC#  ORK  I,  AU,  8u,  CU 
CuMNUN/LR/AW(9)#AR(  IS) 

An  (  9  )  •  j.H> 

AR  I  I  a ) ■ 1M> 
kcAJ<2*800> 

RcA0(2#8j1)XH»XP»YP»ZP»AN*BN»CN,YA*ZA 
RtAu(2#8v.l)XHl,XPi,YPl,ZPl,ANl,BNl#CNi,YAl,ZAL 
Rt AD(2#800 ) 

RE  AD ( 2  1 802  )  BO 
RtAOU#  800) 

RE  AD ( 2#  002 ) B01 
RE  AD ( 2, 800  ) 

800  FORMAT ( ) 

801  FORMAT  (2X»<*(F10.‘*,2X)/6(2X»F1u,4)  ) 

822  FORMAT  12X,3(F10.N,2X), 12) 

C  AN 3  *  8N  S  SUN'S  DIRECTIUNAL  ANGLES 

C  ANS  cLt  I/AT  ION  FROM  Z-  AXIS*  BNS  AZIMUTH  FROM  X-AXIS  TOw  ARO  Y-AXIS 
C  START  SUN  ANGLES#  ANS — 5  DEGREES#  BNS  — 3  DEGREES 
K1NC-2 
KFULL'180 
ORM-KFULL/KINC 
B*dO 

RCC*.a*(Xrt*A2»B**2)/B 

CALL  C0MPR(1,2«RCC*3#XP»YP*ZP#AN#BN»CN#YA,ZA) 

B 1 ■ BO  1 

RCCi* •  a  *  ( XH1**2*B1**2) /B1 

CALL  C0MPR(3#3#RCC1#B1#XP1> YP1#ZP1» AN1#BN1#CN1,YA1,ZA1) 
tNCOOtl 75,995,A*)B,RCC,B1,RCC1 

995  PORMAT ( 2  X# "S  IDc  h  INOOm"* 2 ( 2 X* F 10 . A ) *  2 X # " T OP  W I NDUk" , 2 ( 2 X , F 10 . A  )  I 
CALL  C  OMP  0 ( 1 #  3 ) 

RE  TURN 
tNO 

SUBROUTINE  COMPK(lJl*IJ2*RCC#B*XP#YP*ZP*AN,BNfCN,YA#2A) 
CUMM0N/CYl/NCY,NSC(1G),NSP(10»10),XC( 10 )» YC (10  )»ZC(10)#AE(1C)»8E(1 
UC ) #  CE ( 1 0 ) #  RC ( 1 0 ) 
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DU  1C  lJ«lJi,IJ2 
KCl  I J ) • RCC 

XC(1J)»UP-(RCC-B)*AN>*<<-1I**IJ) 

YU  UJ)-YP-<RCC-B)*BN4YA 
1C I1JI *ZP-(RCC-BI*CN,ZA 
10  CONTINUE 
RETURN 
cND 

SUBROUTINE  C  OPIPQ  (  IJ1,  I  J2) 

CONNON/C  AN/NT,  NB,NC,NP,NA,ND»NmOC),PXV<  100,  6 )  ,  P Y V ( 1 Ou » 8 )  ,PZV(lCu 

0  ,  o  > 

iJMMUN/CYL/NCY,NSC(10),NSP(10,10),XC(10),YC(lG),2C(10),AE(lU)»8E(l 

Jj),CEU0),RC(10) 

COMMUN/ SUN/KlNC,DRKI,AU*BO»CU 
C0Hi10N/SUL/NKP,NRT  (10,10,20) 
luMNUN/FACT/SX,  AO,BO,Al,t»i 
CuMNUN/LR/AW(5),AR(i:>) 

UIHcNilON  L  A  9c  L  (  4  ) 

J  A  I  A  LABEL  /"SMYTH", "B520»,"X3b54",  "HE  3  7V"/ 

OATA  Pl/3. 14139/ 

ICC-L 

R  t  m  1  NO  4 

OU  9  1RU»1,KINC 

XlU* 1 R J- 1 

AN3-XIU*uRKl*:>, 

Ali»ANj*Pl/180. 

00  4  KU ■ 1, K INC 
XRU-KU-1 
BNi*  XRUAORM 
A22»3NS*PI/ibO. 

AU--SlN(All)*CUS(A22) 

3G*-SIN(A11)*3IN(A22) 

CU--COS ( All) 

1CC-1CC+1 

DO  7  1J-1J1,  1J2 

1 C  «  NSC (  lJ  ) 

Du  1  1 K ■ 1 » 1C 

LU'NSPI IJ,  IK) 

CALL  DRVO ( I J, I Q ) 

NR  T ( IJ, IO-NC, ICC ) -NRP 
3  CONTINUE 
1  CONTINUE 
9  CONTINUE 
CALL  URA* 

CALL  DRW i 1 3 ) 

A0--3 • 3 
B  v *  1 00  • 

A1--111.32 

31-100. 

CALL  DRWPU,IJ1,IJ2) 

Al-104. 

31-433. 44 

CALL  DRwPi 2,  I J 1, 1 J2) 

AD--1U.32 
80*690. 42 

CALL  DRWPt  3, I Jl, I J2) 

I C T  »K I NC *K  INC 

CALL  PLTSYNI .123, AW, 0., 600.,  5  •  ) 

ENCODE!  132,1020, AR) ( ICC, ICC-1*  ICT) 
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1020  FuKMAT ( 2X» "SUN  AN&LES"»2Q< 2X»  14  )) 

CALL  PLTSYMt .125, AK, 0 . , 600 . » -7 . ) 

yck-u. 

OCi  lo  I J  *  I  J 1 ,  l  J  2 
IC-NSC(IJ) 

00  10  IK-1, 1C 
IQ-NSP(  IJ,  IK) 

yck-yCk*i2. 

cNCuDc (  132,1021, A*) IQ, (NRT< IJ, IQ-NC, ICCJ,1CC-1,ICT) 

1021  FORMAT  (6X,21(2X,  I  4  )  ) 

CALL  PlTSYMC .125,AR,0.*600.»-7.-YCK) 

10  CONTINUE 
CALL  PLTPGb 
00  25  I  J •  I  J  1 »  I  J  2 
IC  *  NSC (  I  J ) 

00  25  IK- 1, IC 
IO-NSPI IJ, IK) 
tNcUOt ( 40, 1C22, AW ) 10 
iG2  2  FORMAT (2X,NWlN0UW",2X, 12 , 2  b  X ) 

CALL  PLTdEG(5.»5., 1., 13, LABEL) 

CALL  pLTiCA(.25,.25»0.»0.,40.,40.) 

CALL  PLTAXS<  10.,1G.,0.,180.,C'.,180.,8) 

RcWiNO  4 
ICc  «c 

00  21  1R0- 1, KiNC 
00  21  KU*l»KlNC 
ICC-ICC+1 
00  21  I J J» IJi, I J2 
ICE- NSC (  IJJ  ) 

00  21  IKK-l,ICfc 
IOO-NSP ( I J J,  IKK) 

NRP-NRTI I JJ,  IOQ-NC, ICC ) 

IF(nRP.EQ.O)GOTO  21 
00  2w  I-l, NRP 

REA0(4,1w24)N,AVN,AWN,AN,  IS, XR, YR, ZR, BAT  A 
10  24  FORMAT (1 1C, 3F 10. 4, 14, 4F 10. 4) 

I F ( IOQ.NE .  IOJGUTU  20 

XX- A VN 

YY-AWN 

CALL  PLT0TS(3,ICC,XX,YY,1,0) 

20  CONTINUE 

21  CONTINUE 

CALL  PLTjYMI . 125, AW, 0., 40., -5. ) 

CALL  PLTPUE 
25  Continue 

IF(lJl.E0.IJ2.AN0.IC.t0.1)RETURN 
ENC 002(40,1025, AW) 

1023  F0RHAT(2X,"ALL  SURF  AC ES", 26X ) 

CALL  PLTBEG(5.,5.,1.,13,LABEL> 

CALL  PLTSCA(.25,.25,0.,0.,40.,40.) 

CALL  PLTAXSU0.,10.,0.,180.,w.,180.,8) 

KtWlNO  4 
1CC-0 

00  27  I  RU- 1 , K l NC 
00  2 1  KU-1,K1NC 
ICC-ICC-l 
00  2  7  IJ-  IJI,  IJ2 
1C-NSC ( I J ) 

00  27  IK- 1,  IC 
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1Q-NSP  « 1 J» 1K» 

NRP-NRT  ( I  J,  1  Q-NC  *  ICC  ) 

IF lNRP.Ew.O)GOTO  27 
DO  2b  1 ■ 1  *  NR P 

REa01A»102A)N»AVN/AWN,AN»IS,XR»YR>ZR>BA1A 
AX» AVN 
Y Y»  AWN 

CALL  PLTDTSU,  ICC/ XX,  YY,x,G) 

2b  CONTINUE 
27  CONTINUE 

CALL  PLTiYMT  .12b,AW»C.»A0.,-;>.) 

CALL  PL  I  PC  t 
RcTOKN 
C  NO 

SodRQUTINE  UR«P l I C  S » I J 1 » 1 J 2  ) 

C  GMMUN/C  AN/NT  , NB,NC,NP,NA,ND,Ntf(l00),PXv(10C,B),PYV(lUC,6),PZV(lC0 

U,8) 

COMMON /C YL/NCY# NSC (lo »»NSP 1 10,10, Ad 10)/ YCT1C),ZC  C 1C), At  1 10),Bt< 1 
00 ) , t  E( 10 ) , RC ( 10  I 
CONMON/SUN/lUNC»ORKI,  au>  80,  CO 
COMMON/SUL/NRP.NRT (10# 10, 20  I 
COMMON/ F ACT/SX, AO*  80,  Al, 31 
REWIND  A 
ICt«0 

00  20  IkU«l,KINC 
DO  2D  KU»1,KINC 
ICC«iCC*l 
00  IS  1  J  *  I  J  1  /  I  J  2 
I C  *  N  Sc ( IJ) 

00  ic  IK* 1, 1C 
IO-NSP  (  I  J,  IK  ) 

NRP-NRTT  IJ,  IO-NC#  ICC  ) 

IFINKP.EQ. 0)6010  12 
00  1 w  1*1, NRP 

REAim,1002)N,AV,AW,AN,IS,XR,YR,ZK,3ATA 
10o2  FURMAT(I10,3Fi0.A,lA,4F10.«) 

XI*  Yk 

I F ( ICS.Eo.3)Xl*ZR 
Y1*ZR 

IFCICS.Gt.2)Yl*XR 

XX*  I  X  1  A 0  )  * S X *  80 

YY»( Y1*A1  )*SX  +  B1 

CALL  PLT0TS(3,ICC,XX,YY,1,C) 

10  CONTINUE 
12  CONTINUE 
lb  CONlINOt 
20  CONTINUE 
RETURN 
END 

Sod ROU I INE  kt AO V 
C  REA 0  In  SURFACE  VERTICES 

COMMON /CAN /NT »N8,NC,NP,NA,N0,NV(10o) »PXV< 100,8 ) , P YV ( 100,8)  ,PZV (IOC 

0,8) 

LOMMON/VcRT/XV(20C)»YV(230)«2V(230),NVR(100,b) 

C0hM0N/CYL/NCY,NSC(i0j,NiP(lC,lQ>,XC<J.0)»YC(10)»ZC<lC),AEtlC)»BEti 

wO),CE(10),RC(10) 

COMMON/ PILOT/XO, YO,ZO 
REA0I2, 1000) 

1000  FURMATO 
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ReAUtZ,  1000) 

REAQ(2,1001INTfNB#NC»NP,NA,ND 
lCti.  i-URMAT<6(2X,i3l) 

RfcAQU, 10001 

Re  AO (2# IOC  2  I (NV (11*1-1* NO) 

XD02  F0RMAT(10(2X,I3)  ) 

Re  AO ( 2* 1G00 ) 

DO  ID  J  •  1 »  B 

READ(2»lO02) (NVR ( I  *  J  )  * 1  *  1  *  ND ) 

1-  CONTINUE 

RL  AD ( 2* 1000 ) 

RE  AD ( 2* 1003) (XVII )  * l *  1 / NT ) 

1003  t-  jRMaT  <B(2X,F  7.  A)  ) 

Rt  AD l 2  » 1000 ) 

RtAOU»  1003HYV(l)»I-i,NI> 

Re AD( 2# 1000) 

ReA0(2»10vj3)(ZV(l)»I*l»NT) 

R  £  AO l 2  » 1000 ) 

R  r.  A  b  (  2  *  100  2  )  NC  Y 

KtAOU*  1002)  (NSC  (  I),  1-1  ,NCY) 

OJ  20  I  *  1  *  NC  Y 
KP*NSC (  I  ) 

2  0  Rt.  AO  12#1002><NSP(I*K)#K-1,KP) 

KLAD<2,lCb4)(XC(i>,YC<I)»ZC<l)»A£<I),BE(I)»CE<i),RC(I)»l-i,NCY) 

1004  F  OR  MAT ( 7 l 2  X»  F  7 . 4  )  ) 

Re  AO ( 2/ 1000 ) 

R£AD(2*1003)X0»YU»ZQ 
Re (URN 

c  ND 

SUBROUl'lNe  NGRML 

C  e  S  T  ABL I  id  SURFACE  NORMAL  FOR  EACH  SoRFACb  NORMAL 
C  SURFACE  NORMALS  OiRcCTED  TOWARD  COCKPIT  INTERIOR 

COM  MON  /CAN/NT,  NB , NC > NP , NA, NO* N V ( ICC ) , PX V ( 1 00,  B),PYV<1  CO  »B)*PZV(100 
D»  0  ) 

CQMMUN/VeKT/XY(200)*YV(20U)»ZV(20C)#NVR(iOC*&) 

COMMON /NORN/ A XN (100  I # ATM ( 100  I »  AZN( 100  ) 

COMMON/CYL/NCY/NSC  « 10)»NSP ( 10*10)# AC ( Id), YC (Id  )  #ZC( ID  )»AE (  1C  )  ,  BET  1 
00 ) »  C E ( 10 ) *  RC ( 10 ) 

CGMNuN/PlLOT/Xj* Y3#ZC 
uj  10  I • 1 #  NO 
NK • N V  T  I ) 

DO  t>  K»1#NK 
KV-NVR (  I*K  ) 

PXV( I#R)«XV(KV) 

PYY(i»K)-YV(KY) 

PZV(1»K)-2V(KV) 

5  CONTINUE 
K  ■  2 

7  A  1 ■ P  X  V ( I#K)-PXV< I#i> 

A2-PXV(I»Ka1)-PXV(I#1) 
d 1 • P  Y  V ( I#K)-PYV(I»1) 

B  2  •  P  Y  V  (  I#  K  +  l  1  -  P  Y  V  (  I*  1  I 
C 1 • P  Z  V ( I  #  K  )  -P  Z V ( l  »  1 ) 

C2-PZVT I#K*1)-PZV(I»1) 

Pi-SJRT ( A1**2+B1**2»C1**2 I 
P 2* SORT (A2**2+B2P*2*C2**2) 

A-( Ai*A2tBi*B2ACl*C2)/(Pl*P2) 

I F< ABS ( A ) . L T .1 . I  00  TO  9 
K  »K  ♦  * 
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IRK.Ed.NM  Uu  TO  10 
GO  (0  7 
V  AN  *  A  C  OS (A) 

K-1./(P1«P2*SIN(AN>I 
AXNm»*<81*C2-CU82l*R 
AYN(II»-Ul*C2-Ci*A2)*R 
AZNU)»RA1*B2-A2*B1)*R 
iO  CONTINUE 
kc I  URN 
t  NO 

iUcsRQU  TINE  IRSC Y l 1/ Xk, Yk/ Zk > 

C  CONVERTS  COORDINATES  IN  CLINDRICAL  C00RD1NAIES  INTO  RECTANGULAR  SPACE 
C0HM0N/CYL/NCY,NSC(10)/NSP<10/10)/XC<1G)/ YC<10)/ZC  <1G),AE(1C)/BE<  . 
00)/CET10)/RC<10) 

CN*Sa«r< 1 . -C  £ ( I ) **2  ) 

AN--AE  (I)*Ct  (D/CN 
dN*-6£(l)*Ct(il/CN 
AP-Bt( I ) *CN-BN*CE I 1 ) 
dP»-UE<I>*CN-AN*CETI)) 

C  P • AE ( 1 ) ♦BN-AN*BE  < I ) 

Ru*  AE  ( 1 ) * ( XR-XC ( I) )  +  BE (I)*(YR-YC( 1 1  I ♦ C  t  C I )  *  { ZR-ZC (I  )  ) 

XX*XC(I ) ♦ AE  ( I ) *RQ 
YY»YC(I)*8E(  I )  *R0 

zz«zcm+c£m*Ro 

k»SukT(UR-XX)**2RYR-YY)**2*(ZR-ZZ)**2) 

AA« ( Xk-XX  )  /R 
Bb« ( YK-YY)/k 
CC»(ZR-2ZJ/R 
A«AMAA*dN*8B*CN*CC 
IK  A.GT.i.  >  A  *  1  . 

IFIA.LT. -l.)A«-l. 

ANu« ACJS ( A ) 

U*  AP*AA  +  BP*BB  +  CP*CC 
IRQ.LT.Q.  )ANG  —  ANG 

xr-ang*rci  n 

YR-RO 
ZR-R-kC ( I ) 

RETURN 

cND 

subroutine  ORvu(iCzii) 

CQMNGN/CAN/NTzNd/NCzNP/NAzNDzNV(10G)zPXV(100/8)/PYV(lG0zd)/PZv(10Q 

C/8) 

C0MMQN/CYL/NCY,NSC{10)/NSP(10zi0)/XC(10)/ YC<10)zZC(1C)zAE<10)z8E(1 
CO ) / C  E ( 10 ) / RC ( 10 ) 

COMMON/LINc/ AS/ BS/CS/XS/ YS / ZS / AC / 8C / CC 
CuMMuN/PILOT/XPz YP/ZP 
COMMON/SUN/KlNCzDRKl/ ASS/BSS/CSS 
COMMON/ SOL /NRP/NRT (10/ 10/20) 

DAIA  PI/3.1Ali»V/ 

NkP  *0 
Yi« YP*20. 

10  CONTINUE 
XS*0. 

ZS«100. 

YS»Y  S“5  * 

IF( YS.LT .0. )RETURN 
DU  20  IA*l/180/2 
An» I A-91 
AN“ AN*P I / 1 80 • 
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AS-SIN(AN) 

8S*0. 

C  S-CUS ( AN  I 
ISR-0 

CALL  INTCYIIC»ISK»  IS#XR# YR»ZR) 

IF l ISR.EQ .0) GOTO  20 

BATA*ASS*AC*BSS*BC>CSS*CC 

IFIBATA.LE.0.1G0T0  20 

AR*ASS-2.*AC*BATA 

8R*8S$-2.*BC*8ATA 

CR*CSS_2«*CC*BATA 

IFICR.GT.1.  ICR<1. 

IFICR.LT.-l. JCK--1. 

AVW*ACQS ( C R  ) 

ASN-SINtAVV) 

IF  I ASN.EQ.O. )G0T0  18 
ASN-BR/ ASN 
IF (ASN.GT.l. JASN-l. 

IFUSN.LT.-1.)A$N*-1. 

AHH  »  AC  OS  (ASM 
oO  TO  19 
la  CONTINUE 
AHH*  0. 

19  CONTINUE 
NkP«NRP*l 
AV*AVV*l0U./PI 
AH*AHH*180./PI 

rfRlTE<A»1002)NRP#AV,AH,AN,IS,XR,YR,ZR#BATA 
IOC  2  FORMAT  I I10#3F10.A»  14.4E10.A  ) 

2C  CONTINUE 
GO  TO  10 
END 

subroutine  i ntcyi i# isk# is# xk, yr#zr  > 

C  CuMPUftS  INTERSECTION  POINT  OF  LINE  WITH  CYLINDER 

COMMON/C  AN/NT#NB»NC»  NP»NA,  ND#NV  <  100  )#PXV(10G»  8  ),PYVUOO#  6  )»PZV(  100 
C  #  8 ) 

CUMN0N/CYL/NCY,NSC(10)#NSP(10»10)#ACU0>,YC(10»#ZC(ig  » # AE 1 10 ) » bE ( 1 
C0)#CE(10)#RC(10) 

CQMMUN/LINE/AS#8S#CS#  XS» YS» ZS» AC#  BC»CC 
COMMON/PILOT/XP#  YP»ZP 
C  DETERMINE  INTERSECTION  POINT 
xo  *XC ( 1 ) 

YO*  YC ( I ) 

zo*zcm 

RO-RCII) 

AO-AcI  I  ) 

BO*  BE (  I ) 

CC-CE(I) 

k OS* SORT! C  XS-XO  »  * * 2 ♦ « YS-YO ) **2> ( ZS-ZO » **2  ) 

AOo*(XS-XO)/ROS 
BOO-I YS-YOI/ROS 
COO- ( ZS-ZO I / R  0  S 
A»AO*ASm*BStCG*CS 
B«RCS*(AO*AOO*BO*BOO*CO*COC ) 

A1*AS-A*A0 

A2*BS-A*80 

A3*CS-A*C0 

ttl*A00*R0S-8*A0 

B2«B0o*R0S-B*B0 
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**53*'.!* 


B3«C00»R0S-B*C0 

A-Al**2*A2*«2*A3**2 

B«A1*81*A2*B2*A3PB3 

C-B1**2+B2**2»83**2 

AB--B/A 

AX*B**2-AMC-RQ**2I 
1  FI AX.LT.O. » AK«0. 

B8-SQRT ( AX)/A 
RS-AB-88 

IF  US.LI.Q.  )RS«AB*BB 

xr-xs»as*rs 

YR» YS*BS*RS 
ZR-ZSKS*RS 

R«RS*(AO*AS^BO*BS*CO*CS)tRCS*(AO*AOO+BO*BOO^CO*COO) 

X 1 *  X  Q+  AO  *  R 
Yl»YO+BO*R 
Z  WO  +  COAR 
«C* ( Xl-Xk  )  /RO 
3C«(Yl-YRI/RQ 
CC»(Z1-ZR)/R0 

C  Ob  T  E  RHINE  nHblHER  RAY  STRIKES  ENCLOSED  SURFACE 
XSS«  xP 
YSS*  YP 
ZSS-ZP 

CALL  TRSCY(I*XSS*YSS*ZSSI 

XRR  «  XR 

YRR-YR 

ZRR-ZR 

CALL  TRSCY(i,XRR,YRR,ZRR J 

PI-XkR-XSS 

P2» YRR-YSS 

P3-ZRR-ZSS 

X1»PXV(1S*1» 

Yl*PYV< IS, 1» 

Zi-PZV ( IS* 1 ) 

CALL  TRSCY(I*X1*Y1*Z1) 

Al-Xl-XSS 
bl» Yl-YSS 
C  Wl-ZSS 
IN-NV(IS) 

00  10  1 1  ■  1  *  I  N 
IC  -  1 1 ♦  1 

IF( II.tO.IN)  IC - 1 
X2*PXV l IS* 1C  ) 

Y2-PYWI IS* IC  ) 

Z2«PZW(  IS,  IC  > 

CALL  TRSCY(I*X2*Y2*Z2) 

A2-X2-XSS 

B2-Y2-YSS 

C2-Z2-ZSS 

Q«P1A(81*C2-B2*C1  )-P2*(Al*C2-Cl*A2  )*P3MA1*B2-B1AA2) 
IMO.br.O.  (RETURN 

c  ray  strikes  surface  on  enclosed  sioe  of  surface  eoce 

A 1*  A2 
31  *  B2 
Ci«C2 

10  CONTINUE 

C  RAY  STRIKES  ENCLOSED  SURFACE 
I  S  K  •  1 
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RETURN 

tNO 

SUBROUTINE  DRAM 

C  DRAMS  graphic  picture  op  canopy  in  three  polo  layout 

COHHuN/ORA/NXUOOI»NS(100)»XS(1000)>YS<1000)#N5Y<100)#NE#NX 
C0HN0N/CAN/NT,NB,NC,NP,NA,ND,NV<100)»PXV(i00# 8)#PYV( 100# 8 ) » PZ V ( 1GG 
Q#B) 

CONNUN/NQRH/AXN(lO0l*AYN(lO0l»AINUOOI 
CON NON/PILOT /X0«YO»20 
COHMUN/FACT/SX»AO»BO» A1#B1 
NX*  0 

SX-2.122 

NSS-NC*1 

A0--5.5 

60*100. 

Al  — 111.52 
b I* 100. 

XX- I  YO*AO)*SX  +  eO 
YY  -  ( Z0»A1 ) *SX»B1 
CALL  BeLT(1#XX#YY#1#0) 

CALL  cHARKC  XX# YY I 
NE-1 

DO  2  1-NSS.ND 
QS--AXN<I ) 

IF  (  (I.LE.NC.AND.OS.LT.O.)  .OR.d.GT.NC.AND.US.GT.O.  )  )GOTU  2 
C  Surface  FACES  WIEMEk  PROM  SIDE  VIEW 
XX-(PYV(1#1)+A0)*SX+B0 
YY- (P2V(I»1)+A1)*SX+B1 
N  E  *  N  E  ♦  1 
NO*  1 

IF ( 1.LE.NC.0R.I.E0.NAIHD-4 
call  belt (ne#xx#yy#md#o> 
call  L INS (  II 

2  CONTINUE 
A1-1C4. 

61-433.44 

XX«(Y0*AC)*SX+B0 

YY-(X0*A1)*SX+B1 

NE-Ne*1 

CALL  BELT (NE#XX#YY#1#0) 

CALL  ENARK ( XX#  YY ) 

DO  3  l-NSS.NO 
OT ■♦AZN ( I  ) 

IF U I. LE. NC.AND.QT.LT. 0. > .OR. I  I. GT.NC. AND. 0T.G1 .0.) ) GOTO  3 
C  SURFACE  FACES  VIEWER  FRON  TOP  VIEW 
XX-  IPYVII#1)*A0)*SX+B0 
YY-(PXV(I#l)*Al)*SX+Bl 
NE-NE*i 
HD-1 

IF  < I.LE.NC.0R.I.E0.NAIND-4 
CALL  BELT INE»XX»YY»NU»GI 
CALL  LINT(I) 

3  CONTINUE 
AO— 111.52 
Bu-690 .42 
XX-(20>A0)*SX+BO 
YY-C  XOtAl  )*SX*B1 
NE-NE*1 

CALL  BELT (NE#XX#YY#1#0) 


no 


CALL  E  TIC ( XX* Y Y  ) 

00  10  l  -N  S  $  *  NO 
QF--AYN< I » 

IF  ((  I.  L£.  NC.ANO.QF.LT  .0.  »  .OR.  (  I .  GT  .  NC  .  AND  .  OF  .  (i  T  .  0  .  t  (GOTO  1C 
C  j URF  AC  E  FACES  VIEWER  FROM  FRONT  VIEW 
XX  ■ ( P  Z  V ( I*1)*A0(*SX»B0 
YY-(PXV(I»l>*Ai)*SX*Bl 
NE ■ NE ♦ 1 
MO-1 

IF( I.LE.NC.0R.I.E0.NA1MD-4 
CALL  B  t  LT (NE*XX#YY»MD»0) 

CALL  LINF(I) 

10  CONTINUE 

NS(NE*i)»NX*l 
RETURN 
t  NO 

SUBROUTINE  BELT(1E*X*T*MDD*NSS) 

COMMON/ OR A/ MOT I00)*NS( 100),XS(1000) »YS(1000> #NSY C 100)*NE»NP 

NP-NP+1 

NS  ( IE ) - NP 

NS Y ( IE ) -NSS 

M  J  (  I E  )  - MOD 

XSINPT-X 

YS(NP»-Y 

Rt TURN 

C  N  j 

SUBROUTINE  PLEM(X,Y) 

COMMUN/ORA/MDl 100 )»NSl 100 )*XS( 1000 )*YS (1000 )»NSY( 100)* NE#NP 

NP • NP ♦ i 

XS1NPI-X 

YS(NP»-Y 

RETURN 

END 

SUBROUTINE  eMARK ( X*  Y ) 

C  MARK  EYE  POSITION  IN  CANOPY 
X-X-10. 

Y-Y-5. 

CALL  PLEMi X* Y ) 

X-XtlO. 

Y-YO. 

CALL  PLEM(X*Y» 

X-X-10. 

Y- Y ♦ 5 • 

CALL  P L E M ( X*  Y I 
X-X»3. 

Y-Y-1.5 

CALL  PLEM( X*  Y ) 

Y-Y+3. 

CALL  PLEM( X*  Y ) 

Y-Y-i*. 

CALL  PL  EM ( X*  Y ) 

RETURN 

ENO 

SUBROUTINE  ETIC(X*Y) 

C  MARK  EYE  POSITION  IN  CANOPY 
Y-Y-5  # 

CALL  PL  EM ( X*  Y ) 

Y • Y ♦ 10 • 

CALL  PLEM(X*Y) 
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y«y-5. 

CALL  PLEMXfY) 

>•*-5. 

CALL  PLEMUfV) 

X«X*10. 

CALL  PLEMIX#  Y ) 

RETURN 
E  NO 

SUbRuuTINE  LlNS(I) 

C0«hUN/CAN/NT,Nd#NC*NP»NA#ND»NV«i00)#PXV(lut,8),PYV(  100*  (3 )  ,  Pi  V  ( 1 00 

QfB) 

CuNMUN/FACT/SX, AOfBOf Ai# Si 
NK • N V ( I  I 
00  10  R-i#NK 
Ki-KU 

IF(K.EQ.NK)K1-1 
XX*(PYV(1fK1)+A0)*SX+BU 
YY-(P2V(I#K1)«A1)*SX*B1 
CALL  PLEN(XXfYY) 

10  CONTINUE 
K  t  r URN 
c  NO 

SUBROUTINE  LlNF(l) 

CaH«0N/CAN/NT#NB»NC,NP,NA,ND#NY(100)»PXUi0C»8)#PYV(lC0Fb),PZV(lGC 

a»e > 

connon/fact/sxfao#bo#ai#bi 

NK,» NV (  I  I 
00  In  K« 1 f  NK 
Kl-KU 

IF  IK.cO.NMKi-l 
XX> (PZV( I#K1)+A0I*SX*BL 
YY-IPXVC XfK1J*A1)ASX»B1 
CALL  PLENIXXfYY) 

10  CONTINUE 
RETURN 
cNO 

Subroutine  linhii 

C0«n0N/CAN/NTFNBFNCFNPFNA,N0FNV(100)FPXVli0CFS)FPYV« IOOf d ) f P Z V  1 100 

0,0) 

CUNNQN/F ACT/ SXf AOfBOf AIfBI 
NK'NVIl) 

00  lu  K« If  NR 
K 1 ■ K  *1 

IF(K.cO.NK)Kl-l 
XX«(PYV(IfK1)*A0>*SX*B0 
YY«(PXW(I»K1)*A1)*SX+Bl 
CALL  PLEH(XXfYY) 
lu  CONTINUE 
Rt  I  URN 
END 

SUBRUUTINE  DRWUUNIT) 

COMMuN/ORA/HO( 100 )fNS ( ICO )fXS( 1000 )fYS( 1000 )fNSY( 100 JfNEfNP 
CUMMUN/L  B/ STRGI 9 ) 

C0NM0N/D6/XPU00I#  YPdCOI 
OIHtNSIQN  LABEL  ( 9  I 

DATA  LABEL /"SHYTH"f"B520"f"X3659"f "HE  379"/ 

C  PROGRAM  CQMPUTtS  LAYOUT  ON  1029  UNITS#  WHICH  IS  FITIEU  TO 
C  25  INCHES  BY  25  INCHtS  DISPLAY  SHEET 

CALL  PLTBEG(25.#25.#1.# IUNITf LABEL ) 
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WALL  PLT  SC  A15.*5.»200.»0.»40.96#  40.96) 

CALL  PLTSYHI  .25»STRS,0.*20l.#l.  » 

00  10  IE«l*Nt 
MOU-flD (  It  ) 

NSS-NSY (  IE  ) 

Ni«  NS  1 1 E  ) 

N2-NST lt*l»-Nl 
OQ  5  1*1,  M2 
XPl I )»XS(N1*I-1) 

ypu  i-ystnim-i) 

3  CONTINUE 

Call  PLTDTS(MDD»NSS»XP»YP»N2»0) 
lu  CuNTINUt 
RETURN 
ENO 

C  A AH  CANOPY  DATA 

C  NUMtftR  Of  VERTICES#  SURFACES --CONSTRAINT*  FLAT*  AND  CYLINDRICAL 


1  bo 

lb 

56 

5b 

39 

64 

NO. 

VERTICES  I 

PER  SURFACE 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

c 

3 

t 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

7 

7 

7 

7 

4 

4 

b 

4 

4 

4 

4 

4 

4 

4 

5 

4 

** 

* 

7 

7 

7 

7 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

b 

6 

6 

b 

6 

4 

Vc  R  T I C  c 

ASSlbNcO  TO 

EACH 

SURFACE 

oa9 

090 

066 

121 

125 

129 

133 

135 

139 

143 

14b 

149 

15  3 

157 

137 

145 

001 

002 

003 

004 

00  7 

Owl 

005 

012 

Oil 

021 

015 

028 

022 

032 

029 

033 

110 

035 

036 

111 

038 

039 

038 

109 

040 

v>41 

112 

Ob  3 

047 

060 

054 

Ool 

061 

089 

-94 

093 

100 

099 

098 

067 

077 

081 

101 

066 

0  7i 

0  94 

100 

w85 

117 

1x4 

115 

122 

12b 

130 

134 

136 

140 

144 

14b 

130 

154 

158 

1 3o 

146 

002 

003 

004 

005 

014 

O'o  7 

004 

014 

012 

020 

Olt 

027 

023 

031 

010 

034 

109 

036 

037 

112 

039 

038 

111 

111 

041 

11U 

110 

052 

048 

059 

055 

064 

u62 

090 

D&9 

C94 

>*95 

IOC 

u99 

098 

07b 

082 

102 

069 

072 

093 

099 

Obb 

114 

1x5 

119 

123 

127 

131 

130 

137 

141 

140 

147 

151 

166 

159 

162 

163 

009 

010 

311 

012 

012 

Oc5 

003 

006 

009 

019 

017 

026 

024 

030 

031 

109 

033 

043 

044 

044 

046 

034 

109 

u  3  7 

04  5 

112 

042 

051 

049 

058 

05b 

063 

063 

Obb 

365 

070 

071 

0  76 

u  75 

074 

079 

C83 

103 

070 

>*73 

092 

C  98 

087 

w90 

w  66 

i*65 

124 

128 

132 

129 

138 

142 

139 

148 

152 

155 

160 

161 

164 

006 

009 

010 

Oil 

00  5 

0-2 

002 

009 

010 

016 

Clb 

025 

025 

029 

032 

ilC 

042 

042 

043 

037 

045 

C33 

034 

036 

04b 

045 

043 

050 

030 

057 

057 

Co  2 

064 

065 

070 

069 

u7b 

075 

074 

088 

080 

C  84 

104 

065 

074 

C91 

097 

086 

0 

w 

0 

0 

0 

0 

0 

0 

0 

0 

0 

it>5 

15b 

0 

0 

0 

0 

0 

C 

0 

u 

0 

0 

0 

0 

017 

019 

024 

026 

0 

0 

G4i 

c 

0 

0 

0 

0 

C 

V* 

035 

113 


j  w 

C44  o  49 

051  056 

C5b 

0 

0 

0 

0  w 

C  0 

J  0 

0 

0 

105 

066 

w75  090 

u96  0 

J  w 

o  C 

0  0 

0 

0 

0 

0 

&  0 

0  0 

0  0 

0 

0 

0 

0 

C  w 

0  0 

0  016 

020 

023  027 

0 

0  40 

0  0 

0  0 

0 

0 

0 

0 

0  0 

0  04b 

052  055 

059 

0 

0 

0 

0  w 

&  C 

0  0 

C 

0 

106 

067 

u76  wov 

995  0 

»  u 

0  c 

0  0 

9 

0 

0 

0 

w  0 

0  0 

0  0 

0 

0 

0 

0 

V  0 

0  c 

0  013 

021 

0 22  028 

0 

j  0 

C  0 

0  0 

0 

0 

0 

0 

w  V 

0  047 

053  054 

060 

0 

0 

0 

0  w 

0  C 

0  0 

C 

0 

107 

0 

V  0 

w  0 

j  0 

0  0 

0  0 

0 

0 

0 

0 

0  Irf 

0  0 

0  0 

0 

0 

0 

0 

r 

4  V 

0  0 

C  0 

0 

0 

0 

0 

0  ^ 

9  G 

C  0 

0 

c 

0 

0 

w  V 

9  9 

w  3 

9 

c 

0 

0 

C  0 

0  0 

0  0 

0 

0 

10b 

c 

0  0 

C  0 

VERTICE  X 

-POSITION 

2.0c 

2.0C 

2.00 

2.00 

2.00 

2 .00 

2. 00 

-2.00 

-2.GC 

-2  .Gw 

-2.00 

-2.00 

-2 .00 

-2.00 

19.06 

18.70 

14.25 

12.90 

11.20 

14.76 

17.42 

-lb. 00 

-17.60 

-13.84 

-13.30 

-13.24 

-15.96 

-17.96 

8.50 

-11.30 

-12.04 

9.08 

19.63 

19.63 

11.46 

11.46 

11.46 

19.63 

19.63 

-19.63 

-i.9.63 

-11.46 

-11.46 

-11.46 

-19.63 

-19.63 

19.16 

18.84 

.4.22 

13.10 

11.24 

14.73 

17.50 

-17.90 

-17.60 

-14.00 

-13.30 

-13.2b 

-15.9b 

-17.96 

8.56 

-11.48 

-12.14 

9.16 

23.0b 

20.92 

18.08 

13.72 

22.72 

23.28 

23.20 

13.60 

13.30 

13.16 

16.36 

22.84 

16.72 

-16.72 

-14.2b 

14.2b 

13.60 

-13.60 

-9.20 

9.20 

10.56 

-10.56 

-11.56 

11.56 

-23. cb 

-20.92 

-18. Ob 

-13.72 

-22.72 

-23.2b 

-23.20 

-13.60 

-13.30 

-15.16 

-16.36 

-22.84 

19.63 

13.26 

11.46 

-11.46 

-13.26 

-19.63 

-12.00 

12.00 

13.26 

-11.46 

13.26 

-11.46 

-23.75 

-4. 

4. 

23.75 

-4. 

-23.75 

4. 

23.75 

24. 

36. 

36. 

24. 

-36. 

-24. 

-24. 

-36. 

26. 

104. 

104. 

26. 

26. 

104. 

67.6 

61. 

61. 

67.6 

-104. 

-26. 

-26. 

-104. 

-104  . 

-26. 

-61. 

-67.6 

-67.6 

-61. 

-46 . 

-24. 

-22. 

-5b. 

24. 

46  . 

56. 

22. 

-10. 

10. 

10. 

-10. 

61  • 

61. 

-61. 

-61. 

-56. 

56. 

VtRT ICE  t 

-POSITION 

37.50 

57.50 

68.37 

71.57 

58.98 

58.98 

60.30 

57.50 

57.50 

68.37 

71.57 

58.98 

58.98 

60.30 

64.66 

82.92 

63.46 

86.98 

98.18 

98.66 

96.24 

84.66 

82.92 

63.46 

66.98 

98.18 

98.66 

96.24 

87.48 

87.48 

97.56 

97.56 

98.45 

103.49 

100.92 

110.92 

115.61 

121.40 

121.50 

98.45 

133.49 

100.92 

110.92 

115.61 

121.40 

121.30 

144.26 

142.52 

142.54 

146.00 

157.52 

158.40 

156.32 

144.26 

142.52 

142.54 

146.00 

157.52 

138.40 

156.32 

148.00 

148.00 

136.66 

150.66 

59.20 

59.20 

70.14 

108.42 

110.90 

69.7b 

116.49 

1x3.25 

139.69 

156.97 

156.97 

140.09 

38.30 

58.30 

67.18 

67.18 

69.57 

69.57 

198.01 

108.01 

112.13 

112.13 

145.44 

145.44 

59.20 

59.20 

70.14 

10b. 42 

110.90 

69.78 

116.49 

113.25 

114 


139.09 

150.97 

150.97 

1*0.09 

115.52 

115.52 

' 

11*. 82 

.  1 

11*. 62 

115.52 

115.52 

113.01 

113.01 

103. *9 

103. *9 

121. *0 

121. *0 

57.5 

5.5 

5.5 

57.5 

5.5 

57.5 

5.5 

57.5 

59.2 

05.2 

132.01 

158.0 

05.2 

59.2 

138.0 

132.01 

1 60 • 01 

191.01 

226.01 

230.61 

160.01 

191.01 

189.61 

169.01 

loo  .01 

180.61 

191.61 

186.61 

230.01 

220.01 

191.01 

166.01 

x  d9 . 6 1 

189.01 

180.01 

180.61 

21*. 01 

21*. 01 

21*. 01 

21*. Ol 

219.01 

21*. 61 

21*  .01 

21*. 01 

198.01 

198. 61 

198.61 

196.61 

198.01 

22*. 01 

198.0 

22*. 01 

21*. 61 

21*. 01 

c 

VtRTICE  2- 

POSITION 

131.81 

i 38 . 06 

1*6.27 

1*1.77 

133.97 

133.97 

132.83 

131.61 

13o . Oo 

i*6. 27 

1*1.77 

133.97 

133.97 

132.63 

128. *0 

130.3* 

1*5.20 

1*8.5* 

1*8.08 

130.5* 

127.96 

128. *0 

130.3* 

1*5.20 

1*8.5* 

1*8.08 

130.5* 

127.96 

118.01 

116.01 

150.00 

15o.0L 

129.20 

1*0.00 

157.56 

100.51 

100.68 

139.20 

129. 2C 

129.20 

1*8.00 

157.56 

100.51 

IOC  .68 

139.20 

129.20 

1*7.56 

1*9. *C 

10*. 20 

107.6* 

107.06 

156.20 

1*7. 5* 

1*7.50 

1*9. *0 

16*. 2b 

lo7 . 0* 

107.68 

156.20 

1*7.5* 

137.80 

137.60 

175.68 

175.68 

129.03 

1*1 .*1 

15*. 60 

172. *0 

150.66 

130.65 

1*0.36 

17*. 39 

179.28 

179.2* 

170. *0 

150.60 

1*3.52 

1*3.32 

15*. *0 

15*. *0 

150.72 

150. 72 

175.15 

175.13 

177.93 

177.93 

x  6  2 . 1 2 

182. 12 

129.03 

1*1. *1 

15*. 60 

172. *0 

150.08 

130.65 

1*0.30 

17*. 59 

179.28 

179.2* 

176.00 

150.80 

152.5* 

152.3* 

159.03 

159.03 

152.5* 

152.5* 

175.79 

176.79 

1*8. Oo 

1*8.00 

139.20 

1 39.2C 

1**.2 

129.1 

129.1 

1  **.2 

122. 1 

127.05 

122.1 

127.05 

12*. 12 

12*. 12 

12*. 12 

12*. 12 

12*. 12 

12*. ±2 

12*. 12 

12*. 12 

1*2.52 

1*0.52 

1*0.52 

1*2.52 

132.12 

132 .92 

132.52 

132.52 

111.52 

111.52 

1*0.52 

1*2.52 

1*2. 52 

1*0.52 

132.92 

132.12 

132.52 

132.52 

111.52 

111.52 

1*7.92 

1*7.92 

176.52 

176.52 

1*7.92 

1*7.92 

170.52 

178.52 

181.52 

161.52 

221.52 

221.32 

111.52 

132.52 

111.32 

132.52 

157.92 

157.92 

c 

CYLINDRICAL  DATA 

2  2 

1 

00  01 

62  63 

0* 

-31.995* 

135.89 

1*7. *713 

.0 

.9805 

.1966 

55.991b 

31.995* 

135.89 

1*7. *713 

.0 

.9805 

.1968 

35.9916 

0. 

133.0172 

139.3729 

0. 

-.9d37 

-.1790 

*2.9375 

C  PILOT  EYE  POSITION 


-1.5  1*2 • 33  171.2 


C  CANOPY  EDGES  FITTED  TO  CYLINDRICAL  SIDES 


20.783 

17.875 

130.91 

101.173 

.9592 

.0 

.2828 

*.98 

.9995 

11.25 

c  • 

120.175 

180.133 

0  . 

- • i 796 

.9837 

0. 

0  # 

C  CANQPY  OISPLACEMtNT  OF 

CYLINDRICAL 

SIDES  FROM 

FLAT  WINDOW 

*. 


C  (OP  wlNuUW 
1.5 
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SMYTH,  STMP2, 11000. 

ACCOUNT, HE***. 

BtGlN, ATTACH, PL0TLI8. 

RtUUES  T»  TAPei3, *PP  . 

ftn(sl*r»  . 
nap, qn . 

LGO. 

8EGlN,PL0T,CALCUMP,TAPE13. 

EXIT. 

PEG  IN,  PLOT,CALCOMP,TAPtl3. 

exit. 

FkGuKAM  ORmEI INPOT, OUTPUT,! AFc4,TAPE13, I APt2« INPUT, UPt 3-UOTPO  I  ) 

C  MAINlINc,  RAT  TRACING  P  3K  CANOPY  OPTICAL  cVAloAIIGN 
C  j  DIMENSIONAL  FlaI/CYlIndRICAL  CAnOPY  »ITH  lBSTRUCTIuNS 
CUHN0N/LB/STRGI4) 

J AT  A  STRto/"GPTICAL", "DISTORTION","  ",">»/ 

CALL  RtAuv 
CALL  NuRML 
CALL  DRV 
,1  JP 
ENO 

SudRGUTlNt  RcAuV 
C  READ  IN  jCkE ACE  VERTICES 

COMNUN/C  AN/NT,  N8,Nu  ,NP  ,NA,  ND,NV  (iDL),PXWlL.»'!l,PYV(i.Lw,t.),r2V<i: 
J  ,  tf  ) 

COMMuN/VERT/XV(2GOl, YVl20o)»2V<20. >»NVk( Uu,t) 

l JMMUN/C  Tl/nCY,NjC  (LJ»,NSP»lC,10»,XCIi0»,YCa-»,xC(l>.),At<.0),bc( 
0.),C'E(10),RC(I0) 

COMMON/ P1LO I /Xu, Yu, 20 
REAOU,  lOGw) 

1000  E  uk  M A  T I I 

Rt  AO ( 2, 10C0 ) 

RcAD(2,10«.l)Nl  ,  Nb  ,  NC  »  NP  ,  NA  »  NO 

1001  PuRMATI  6UX,  13  )  ) 

READ ( 2, 1C0C ) 

Rc AO  I  2,1002  I ( MV ( I  J,i-1,ND> 

1002  P  JR  MA  T  < 1 0 ( 2X , 1 3  )  ) 

RcAJ  I  2, 1000) 

OU  1  w  J  ■  1 ,  t) 

RcA0(2,1002) ( N/R t I, J ) , I « 1 , NO ) 

10  CON  TINGE 

Rc  AO ( 2 , 1 000 ) 

RtA0(2,iCC3) ( X V ( II, I ■ 1, NT ) 

1003  P  JRMAT(8CX,F7.4  )  ) 

Rt  Ao ( 2 , luOO ) 

READ <2,1003)1 YV(I),I«1,NT) 

Rt.  adi  2, 1000 ) 

RcAOl  2, 1003  )  (2Vin,I>l,NT) 

K c  AO ( 2 , 1 000 ) 

RtA0(2, 10C2INC Y 
REA0(2,luw2) (NSC ( 1), I *1 , Nl  Y ) 

Ou  2-  1  ■  j. ,  NC  Y 
KP-N^C I  I ) 

20  RcAD(2,1002) INSPI 1,K),K«1,KP) 

RcA0l2,20^A)UCm,YC<I),2cm,Acm,8E(l),CL(i),RCm,i*l,NCY) 

1004  FORMA  I ( 71 2 X,F7  .h )  I 
RE  AO (2, 1000) 

RcA0l2, 10031X0, YU, 2J 
Rc I  URN 
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t  N  J 

SUBROUTINE  NJKML 

C  ESTABLISH  SURFACE  NORMAL  fUK  eACH  PLATE  SURFACE 
C  iURFACE  NURMALS  OlPtCTEJ  TOWARD  COCKPIT  1NT  t  k i OP 

COMMON /CAN/NT, Nb#NC»NP»NA»ND»NVllOQ),PxV  (10C»8J»FYV( i  GO,  o  ) »  PZ  V  (  1  C . 

0  i  D  I 

COMMUN/ VtPT/ XV  T200>,YVi20  0),il/(200  ),NVR(*lL,fc) 

COMMON/ NOP M/AXNtl 00), AYN<  1Co),AZN<1Gu) 

COMMON/  C YL/NCY  /NSC  <lC),NiP(lG#iG),XC<lL),YCUu),ZC(i1.>,AE<lv)»Et.<] 
0C),LE(1G),PC(1G) 

CDPMUN/P1LUT/Xu,  YG,ZG 
o 0  10  1 • 1 #  NO 
NP  *  N  V  (  1  ) 

DO  3  K * 1 » M K 
KV*NVP(I»K) 

PxV<  I»KI«  XVtKV ) 

P  Y  *  <1,K)*YV(KV) 

PZv ( 1 »  K  )  *  z  V ( K  V ) 

3  LONTINOt 
K*2 

1  Al»PA»TI»K)-PXV(l»l) 

A2«P A  <  I, k»1  )-PxV  <  I, i I 
B1 *  P  Y  V ( I,  K )-PY  V  l  I  *  1 ) 

Bc«PYV{I»K*i»-PYVU,ll 

Cl»PZV(l,K)-PZV(L»l) 

C2-PZVT I,K*1I-PZV( I,i) 

P  1  *  S*k  1 (Al**2+6l**2+Cl**2) 

P2*SuRT<A2**Zvdc**2vC2**z> 

A«tAiAAZ  +  Bi*B2v-Ci*C2)/TPl*P2) 

IF  <  a  ti  j  t  A) . LI  .1 .  )  00  10  0 

K«tSti 

iFlis,EO,NK)  Ou  TO  iv 
uL  TO  7 
S>  an*aCUSIA) 

K»1i/ I Pi*P2*S 1N( AN) ) 

A XN  I  1  )  «♦ ( 61*C2-Cl*d2 )  *R 
AYNli )*-IAl*Cz-Cl*A2)*R 
AZMl)*MAl*B2-A2*tti)AR 
l 0  CONTINUE 
R  c 1  URN 
fc  NO 

SGdROUTINc  DRY 
CGMM0N/LK/Aw(-,) 

A*(0)»1H» 

R  t  A  0 ( 2  » Boo  ) 

REAo(2»6<.1)XH,XP,YP,ZP,AN,BiV»CN,Y*,ZA 
RLAD(2,8  31)Xl1i,APi,YPl»ZPl,AM,8Ni.,CNl,YAi.,ZAl 
R=AD<2,30G) 

REA0(2»8o2)B0 
Rt  AO (2* 8 GO  ) 

R c  AO  ( 2  >  802  )  U Ci. 

RtAD(2#800) 
aOO  FORMAT ( ) 

8C1  FORMAT (2X,4(FiG»4,2X) /6(2X>FiJ.A) ) 

602  i-URMAT(2X,3(F10.A»2X),I2) 

8*60 

kCl*.S* ( Xh  +  *  2 ♦ B*  *2 )/6 

CALL  COHPR<i.»2»KCC»8»XP,  Yp,ZP»AN»BN»CN,  YA,<.A) 

6  1  *  6  v  l 
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KCCi».y*U8i.**2»bi»*2l/8i 

CALL  CJMPKd#3»PCCl#ttifXPifYPlfZPi*ANAffc,U»CN*fYAi,ZAAl 
tNCu0t(75/99:>f  A»)  OfflCC#  81*  RCCi 

995  PuRNAT  l 2Xf "SIDE  M  i  80  On  ",  2  l  2  A*  F  1C  .  A  )  »  2  A  ,  '  1  UP  «  I  N0uPM>  2  (  2  A  » 
mPIO.AI  ) 

CALL  CUNPQ 

Pc  I l.pn 
b  NO 

SUBROUTINE  COMPRdJlfiJ2»RCCfB»XPffPfZPfAN,dNfiNfYAf2d 
..OHMON/C  YL/NC  Yf  NSC  (10  I  fNSP  (  lC»iOl>AC<lv>»YcddfZC(lG)*Atllv)f 
Jdf  ( lw) *CE ( au)» PC ( 10) 

00  13  i J •  I  J  1 »  1J2 
PCI  1J1-KCC 

XC  ( IJl»-(XP-(PCC-e)*AN)*( ( -1 ) ♦* I J  ) 

Yl(1J)«YP“(PCC-B)*oN+YA 
l C( 1 j I «/P-(oCC-B) *CN»2A 
10  CONI INUb 
KE I OPN 
tNO 

SuoPuuTlNt  CuMPu 

CUMMUN/CAN/NTfNDfNCfNPfNAfNUft'<V(i.GC)fPXLdGa*fL)fPYV(i....fC)fPAVd,, 
af  cl 

CGMMuN/CYL/NCYfN*C(lL)fNSP(lLfiw)fXCdt)fYL(l.)f4.CdC)fAtdL)fit( 

OwifCbdCIfKCdC) 

COMMON/LK/AB (-  ) 

OIMtNoICJN  AU ( 3  I 
A J i j ) * 1M> 

rebind  c 

NVLL»0 

CALL  CANL(NVLl) 

L  TRANSMISSION 
C  AL  l  OR  Am 
oALL  0R*U3l 
CALL  ORtaP ( If NvLL ) 

CALc  PL TSYM( . 123f A*f . . f oCu.  f S  .  I 
tNCuOc(  L3f  995f  AO) 

9 05  bQRPAT ( 2XfM TRANs MISSION  ") 

CALL  PLTSYM(.125fA0f O.ffcOO.f-S. ) 

CALL  PLTPGt 
C  OPTICAL  DISTORTION 
CALL  DRAm 
CALL  OkM(li) 

CALL  URwP ( 2f  Nv  LL  ) 

CALl  PLTSYhi  .12Sf AWf C.f 6^v.f S.  ) 
tNCuDE(15f 996f AD) 

99o  PURMAT(2Xf "OPTICAL  POWtP") 

CALL  PL TSY«( . 12at AD,0. fbOu. f-5.  ) 

CALL  PLTPSc 

RETOkN 

END 

SUBROUTiNc  CANL(NVLL) 

COMMUN/C AN/Nlf  NBf NC»NPf NAf  NUtNV ( iOO ) f b  AV ( iO»f  B  I  f  PYY ( *Cv  »o  I  » 
QPZVdOO*  o) 

CQMM0N/LINt/AS»6Sf CSf  XSf  YSfZSfACfBCfCC 
CUMMUN/PILUT/ XPf  YP»ZP 
DATA  Pi/3. 1*159/ 

YSG • Y  P  *2  .  • 

X  S  C  ■  u  $ 

2  SO* 1 CC . 
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1C  r i o ■  r ^ ij-j. 

IP l YSO.LT .0.  1  RtTURN 

u  0  30  I  A*  U  1  80.  2 

X3»  a  SO 

ys«ysg 

ZS-ZSO 

Am* i a— s 1 

AN*  AN**- 1/ ISO. 

A3«SlM(AN) 

8S*u  . 

C i*Cu3 ( AN  I 
NC  v.  *  ML  ♦  1 
013  tv  IS*NCC»ND 
I  jK  *  j 

I F ( 1 S , L  T . N A  >  C Al L  iNTtCU.ISK.lS.XR.YK.Zf) 

IPliS.GT.NAlCALL  1NTCY(1C>13IUI3»X«,YR,/R) 

IF l IjK.lT.OIGOTU  23 
Zj  CONT  lNllt 
o 0  T  u  3C 
23  C  ON  T 1  NUc 

K*soKIUx*-XP)**<.  +  iYk-YP}**2+<Zk-ZPl**<:> 

Aj*(XR-XP)/k 
B3*  < Yk-YP ) / k 
C  3  *  ( iR-ZP  1  /R 

X3  *  A K 

Y  S  *  Y  K 
Z3*  Zk 

CaLL  C  A  L  C  <  R  .  I C  *  I  S  »  N  \t  L  L  ) 

30  CJNTlNUt 
u  0  T  u  A  0 
c  NO 

3US  R  OUT 1  Nc  CALC  IRC, 1C. IS.NVLl  ) 

CQNMGN/CAN/NUNB»NC»NP.NA»N0»NV(A0uUPXV(3b..»cUrYVliuu»oU>'ZL 

4i6) 

CuMMGN/PilQT/XG»YG»ZO 

Common /iinc/as>bs»cs»xs»Y3>2s»ac>3C>cc 

CoMMON/lrlw/Il 
DATA  Pl/3.1<il59/ 

IFlCS.Gfc.l.lCS-.VJiN 
LFUS.Lt.-l.  ICS*-. 9999 
iP(A3.fcU»u.)A3*.00C’i 
XS0*AS 

Y  j  j  *  Y  S 
Z  Su*  ZS 

A  5  0  *  A  > 
t»S  )*bS 
C  S  V  =  L  S 
NVLL*NVLL  *1 

lHIS.lT.NAlCAll  *KTtC(iilUl3»XfUYIUZK.Ak»BR>Ck) 
1F(1S»GT.NA)CALl  «  .fCY(lc»lSR#iS*XK»YR*ZR*AH»di>>LK) 

AMO*Ai*AR  +  BS*BR  +  CS  *C  K 
X  F ( AMO.GT  .1.  )AM0«1. 

IF  ( AM* .IT  .-1.  ) A  Mu  *  —  1  • 
xii.  »ACCS(  AHC)*idC.  /PI 
rri  1 1(<.,  iooo)  ri.xs.rsws 

100C  FORMAT  tRPlO.C) 

RKi  l  L  <  3#  10GC  1TUXS.YS.Z3 

MX  *  u 

RO«C. 
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CK*.sa*<T<i.-cs**i ) 

AO*-Bi/CRR 
Bu* AS/CRR 

co*c . 

C  A-DtGRtE  CUNt  DISTORTION  FItLD 
JO  AO  1 R  ■  1 »  5 
AiG» IR 

ASG-ASG*. A*PI/iBG. 

Ri*RD*TAN< Abo) 

A*  J  . 

ACG-12C/1R 
AK*  AO 
BK*  BO 
CK-CO 
AC  K  *  AC  G“1  < 

OQ  AO  K  G  *  A  »  1  b 
A*  A ♦ 2 A • 

IFCA.LT.aCIOGijTQ  AO 
ACK* aCK+AC  T 

ax* a*p  w ieo. 

CA-C-SCAX) 

XXH  »BO*AoO-AO  +  BSw 
XXwi«AiL*CA/XXF 
XXl*(CC*ASO-AO*CSuWxXr 
YYO*-BS0*CA/XXF 

Y  Y  i  *  (  (CO*Aiw-Aw*CSJ)*BSL/XXF~Cbrf)/AS.. 

AY«XX1**2*YY1**2+1. 

t}Y»2»*(XX0*XXi*YY0*YYl) 

C Y»  XxG**c*YYO**c-i. 

D  f  *  SGR  f  (BY**2-A.*AY*CY) 

Cl- (-BY+JY  )  / ( 2 • *  A  Y  ) 

6.«XX0+XAi*Cl 
ai*  rrotYn*ci 

JC-ASO*  (BK^Cl-Cl^Bil-Biv* Un*Cl-A**CKWC J-*(AX*Bi-Ai*BA) 
iFCQC.LT.O.  )GJTU  3S 
Cl* (-BY-OY  )/ 12.*AY  ) 

01-XXO+XX1*C1 
Al*YYOtYYl*Cl 
3 5  CONTINOt 
Ak  •  A  1 
Bk  «  u  A 
Ck*Ul 

XS*  xS.<  +  A1*kS 
Ya*YSO*Ba*RS 
ZSsZSG+C l*Ki 

K«SGKT(CXS-XG)**2MYS-r.»***MZW-»***) 

A j ■ <  XS~XO  I  / R 
bj»(YS-YOI/R 
CS« US-2G)/k 
xS*xo 

Y  S  *  Y  J 
ZS»  l  j 
lBIL*0 

1F(1j.LT.NA)CALL  lNTtC(l»IiR#lB»XRi,YKi,tKl) 

IF ( lG .GT  <  N A  I  CALL  lNTCY(IC»liK#lS»XKi»YRi#ZKA) 

IF l ISk.tU.OGOTO  AO 
N  X  •  N  A  ♦  1 

1FCIS.LT  .NATCALL  RRT£C(lSK»lS»Xki»YRi»ZK.i»AKi»DKi»CK») 

1F(  1S.GT  »  N  A ) C  AL  L  wKTCY(IC»ljR#IS»XKA/Ykl#ZR..,AFi»BRl*CRl» 


Art»Aj*AKi+bb*bRi*CS*CRi 

IFCAM.GT.i.)AM-l. 

IF(A,1.LT.-l.)AN«-i. 

XM»ACaS (AM)*180.  /PI 
R0-RO+ ( XMG-XMI **Z 
AO  C  3N  T  I  No  t 
XN»Nx 

IFtNA.tJ.OIGQTU  A I 
i<j*  jgKl  tkO/XNI 
ii  Rk l T c ( A. aOCG  )kO 
RRIfc(3»10CG)AN»RG 
Rt  l  GkN 
cnJ 

aUbROUllNL  lNTtGllG.lik.iS.XK.fk.Zr.) 

C  DfcURMINcS  IF  RAY  ilRIKcb  CJNYtX  SURFACc 

CGMMJN/C  An/NT.NB.NC.  NP.NA,  Ng.NY  <  i  J*  )  >  PX\,  (  ±  wl  .  b  >  #  P  YY  (  .CC  ,  B  )  ,  Pa  V  »  ..  • 
Jf  o) 

CUMMGN/NGk  M/ A  XN  (  10  G  l  »  A  Y  N  (  xG  .  >  »  A  ZM  *•*  -  1 
CGhnijN/LlNfc/AS»3S.C3»XS»Yj.Zs»AC»bt»v.C 
CK»AXNCIS)*AS*AYi«US)*BS*AZN(I:>)«,CS 
XG*  IG 

ck»uv*ag 

IFtv.K.bt.0.)  RcT  GKN 

C  kAY  3  T  k  I K 1 3  SGRFACe  IN  J  G  I  *  A  R  J  GlRtCilCN 
1*1 

iF(AS.tG,FXY(ij»i).ANO.Yi,Lv.PYV(lJ»i).nNL.i.i.tw.YZV(lj>I))i*Iti 
3*(AAN(Ib)*(PXY(13.1)  “X  S  I  ♦  A  YN  (  1  5  )  *  (  P  Y  V  (  lb#  1  )“Y-)  +  »*ZN(  1  b  J  *  l  P  2  V  1  i  b  >  i 
C  1-13  >  I  /Ck 
AR*AS*S+Xb 
YR«bi*b+  Yb 
ZR»Ci*3tZS 
Ai«PAY( IS.U-Ai 
b  i  *  P  Y  V  t Ib.ll-YS 
C1»PZV( lb. 1 )-ZS 
Pi*  XK-Xb 
PZ  «  YR-Y  S 
P  3-ZK-ZS 
I N»  N  V ( I S  ) 

U'J  1G  I*  1>  IN 
1C*  I  ♦ 1 

IF  (  i . t G . IN )  1C  *  i 
A  2  •  P  A  V  (  I  b  »  1 C  )  -  X  b 
bZ-PYY  US.  IC  >-YS 
C2-PZV<  IS#  iC  )-Zb 

J«Pl*(bl*C2-b2*Cl  >  -  P  2  ♦  (  A1  *12-1  i*A2  )  +  rJ*lAi*B2-Gi<-A£;) 

IMG.GUG.  JktTGRN 

C  kAY  b  Tt>  lKt  b  SGRFaCl  UN  lNCLJScg  jibe  G  F  S  uk  F  Ag  t  tulL 
A 1*  A  2 
bi*o2 
C 1  »C  2 

1-  C JN 1  INGE 

C  kAY  SUIkaS  l NC L ub cU  bURFACC 
1  jA  *  1 

AC ” A AN  Hi  ) 
bC«AYNUS) 

CC-aZN(IS) 
k  t TURN 
tNO 

buBRGuI  INi  InICY(  I.ISK.Ig.XR.Yr./r) 
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C  CuMPuUS  INTtRstCTiuN  POINT  UP  LINE  mI  TH  CYLIMJlR 

CUMMUN/taN/M,Nb,NC,NP,NA,NG,Nv(10tl»PAY(lt/t,B),RYV(10Of''),PZvl.^ 

vita) 

COMMON /LINE/ AS, B3,CS,  XS,  tit  ZS,aC,  Bt,CC 

CONMON/CYL/NCY,NSC<1C),N$P(1C,1GI ,XC  <lC>,rC(iu)»2C(*C)»AtC'.  > »  d  c  l 
■J  J  )  »  CE  ( 10  )  »  RC  1 1C  ) 

CUMMQN/PlLOT/XP»YP»ZP 

C  DfcltKMINt  CYLInOERICaL  StRFACt  *HICH  CuNVtX  iORPACt  IS  A  PART  UP 
Ou  2  I»1,NCY 
KP»NiC ( I  ) 

Ou  2  K"1»KP 

U  <NSP<l»R).fcO.IS)bJlU  3 

2  wGMINUc 
R  t  T  Uk  N 

3  CONTINUE 

C  OcTcRMINt  INTERStC  TION  PUIM 
XU» xt  (  1  ) 

Yu ■ YC ( 1  ) 

Z-»ZC  <  I  ) 

Rv  *  RC (  1 ) 

Aj< Ac (  1  ) 
flO«bc ( 1 ) 

Cu«CE(l) 

RGS«SuRT l <  XS-Au) **2M  YS-Y j I **2+ ( a j-ZC )**2 ) 

aCO* ( XS-x  u  )  / kOS 

60CMYS-Y«)/RwS 

COU- (ZS-20)/R0i 

A»A0*AS+oG*b5*CC*CS 

t)*RCS*(Aw*AUOaBJ*BJj  +  Ct*CC>.  ) 

Al- AS-A*AU 
A2»bS-A*bG 
A  C  S -A*C  w 

bi»aco*ros_b*au 

d2«tsG0*RGS-B*B0 
B3»C0U*k»tS-B*C|j 
A»a1^*2-pa2*a2  +  A3**2 
omAl*Bl^A2*B2*A3*B3 
C«B1**2*Ba**2+B3+*2 
AtJ«-B/  A 

AX»B**2-A*(C-RG**2) 

1  P  (  A  X  »  L  T  •  U  •  )  A  A  •  0  • 

BB-AQRT (AX)/A 
KS-AB-BB 

IFIRS.LT.  0.)RS»At»  +  BB 
XR«XS+AS*RS 
YR«YS*BS*RS 
ZR-ZS+CS*RS 

R»RS*(A0*AS  +  B0*flA  +  C0*CS)+RtS*(AC*At3  +  BC'tlww  +  CvAC»C) 

A1"X0+A0*K 
Y 1*  Y  0*  Bu*R 
Z1-ZG*C0*R 
aC«(X1-XR)/K0 
BC-(  Yl-YM/RO 
CC-(Zi-ZR)/RO 

C  OETERMINtS  taHETHtR  RAT  jlRlKti  tNCLGitO  SURPACt 
XSS-XP 
YSS«YP 
ZSS-ZP 

CALL  rkSuY<i,XjS,Y>S,Z$S) 
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XRR»XR 
YR  R-YR 
ZRk«ZR 

CALL  TRSCY(I;AkR«YKR>ZkK) 

P 1  *  XRK-X  AS 
P2«YkR-YSS 
P  3*  ZRk-Z  SS 
Xi*PXVtIS,l» 

Y 1  *  P  Y  Y  (  I  S  »  1  I 

Z1-PZy< IS»ll 

CALL  TkSCY(I#Xi,Yl»Z*l 

Al*Xl-xSS 

bl» Yl-YSS 

CWi-ZSS 

1N*NV(  IS) 

U  J  10  i  l  *  1  #  1  N 
1C  -  1 1  ♦  1 

IM  11. tG.  INI  IC  *  i 
A2-PXV  l  1a*  IC  I 
Y2«PYV (  IS»  1C  I 

zz*pzv<  is,  ic  > 

Call  TkSi. y  1  i> X2,  i  2,  lZ > 

A2*Xt-XSS 
o  2 *  Y2-Y S S 
C  t  *  Z2-Z  $A 

g«Pl*{ttl*C2-d2*Ci)-PZ*(Ai*CZ-Cl*A2)+P3*<Ai*oZ-r.l»A2) 

IF  (  G  «  G  I  •  0  «  I R  c  T  u  k  N 

C  RAY  STkIKcS  SUPFACt  ON  cNCLuSlO  SIGc  OF  jUkU.;  LOCl 
A 1*  A  ? 
a  i.  *  b2 
Cl  *C2 

i<)  CuMlNOt 

l  KAY  STRIKtS  cNCLOAtLl  jUKFACt 
I  SR*  1 
K  £  I  L  K  N 
t  NO 

AUB  Ruu  1  1 N ti  Ik  jC  Y  (  I»  XP>  Tk  ,  Zk  I 

C  CuNVCkTA  COukOInATcS  In  ClINOKILaL  CUOkOINAI^j  inIL  *tC)ANCLLAK  aPaCl 
CUMNUN/CYL/NCY,  NSC  (xv  I  »NSP  (*«.»lC)»Ak.(10l»Ys.(i.)#AL(*«)AAC(lM»Bt(i. 
Jv),CctlC),KC(iO) 

CN»S0RT(1.-Ct(l)**2) 

AN*-A£ ( I  I  *C£ ( 1  l/CN 
BN«-b£(l)*C£(il/CN 
AP*bc< I)*CN-BN*Ct (  II 
BP  *- l At  (  I  I *C  N—  AN *C  t  (  i  I  I 
C  P  *  A  t  (  II*o  N-AiM  ♦(!  t  (  1  I 

KO*  A  fc  C  I  I  *  t  XR-XC  (II»PtJt(il*(YR-YCIl))+Cc(l)*(2R-AC(II» 

X  X  *  XC ( 1 »*A£( II *R0 
YY*YC(I)+b£(I)*RQ 
ZZ*2C<  1 1  +Ct ( 1 1 ♦  RG 

k  *  Swk  T ( ( XR-XX I •  ♦2MYR-YY>**2+(zR-Z2)**2) 

AA« (XR-XX  I  /R 
fab* ( YR-YY I  /K 
CC  •  ( ZR- Z  Z  I  /  R 
A.AN*AA*bN*B8»CN*CC 
I F ( A ■ GT .1 .  I  A*i. 
i f  <  a . l  r .-1 .  I  A  »-l  . 

ANG*  ACuS ( A  I 
0* AP*Aa*BP»BB+CP*CC 
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LFlfl.LI  .0.  1  A  NG  *-  ANvi 
XR»ANG*kC  (  I  ) 

YR-RO 

ZR-R-RC(I> 

RETURN 
cN  J 

SUBROUT INt  NKTEC lISK*IS»XR*Yk»ZR»Ak»BR*CR) 

C  JHNON/ CAN/NT#  NB,NC»NP»NAtNO>NV(103»*Pxv<H,0#tsl,t'YtM  iOC.fc  lii'hllu, 

C,o) 

COMMON /LlNc/AS.8:>,CS,XS»Yi, ZS#AC,8C,CC 
COMMON/ NORM/ AXN (  i&o)»AYN(  1CU.AZNI  IOC) 

C  uMMON/  1 MQ/ T  i 

DAI  A  TFi*  XN  /  <  1*  1  •  54  / 

l) A  1  A  XT/. 5/ 

C  WALK  KAY  THRUU&H  SURFACE 
T  l«l. 

C  t N Tk  I  SUKFACc 

C«»-(AS*aC*BS*Bu  +  Cj*cu  ) 

IF(CA.GT.1.)CA«1, 

IF (CA.LT.-l.  )CA«-i. 
jA-S-Rl ( 1  •  -C  A*  +  2  I 
SAP-iA/XN 

CAP»SQkT( i.-SAP**2) 

ISr»1 

A*  A  COS  1C  A  ) 

AP»ACOS(CAP) 

-All  r*Afo<AN>A/AP»Ti) 

C  PAS>  TO  OTHER  SURFACE 
B-CA/XN-CAP 
AR» AC*B+AS/XN 
dR»ac*b*dS/XN 
CR*CC*3*CS /XN 

R«(ACMxS-PXV(lS#in  +  6C*lYS-PYVlISjin*CC*U5-PZy,(li,ll)*TH)/LAP 

XR ■ X  5+  Ak  *R 

YR»YS+BK*k 

Zft»xS*Ck*R 

T I • T1*E  XP  < -R  *X  T ) 

C  tXlT  SURFACE 

C  A • - ( AR*AC*bR*bC*CR*CC ) 

IFICA.GT.i. ICA-1, 
iF (CA.LT.-l.  )CA--i. 

SA»S0RT(1.-CAP*2J 

SAP»SA*XN 

1F1SAP.GT.1. JGGTC  1j 
C  RAY  EXITS  FRUM  SuRFACt 
CAP»S0RTU.-SAP**2» 

A* ACjS (Ca  I 
AP-ACOS(CAP) 

CALL  TkANS  1XN*A»AP» T 1  ) 
o»CA*XN-C AP 
AR» AC*bfXN*Ak 
0K«dC*t+XN*Bk 
CK«CC*tUXN*CK 
RETURN 
10  CONTINUE 

C  RAY  LUNlAlNtO  MlTHlN  SURFACE 
T  1*0. 

AK»XN*(AC*CA*AR) 
bft»XN*( BC*CA+BR> 
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uk-XN*(CC*CA+CK) 
t  T  uRN 
c.  NO 

SUBkuUUNt  *KTLY(iC»iSK»i:>>Xk,YR»Z*>*i<»fcK,C*J 

CuMhuN/C aN/NT,  Nd, NC*NP»NA,  NJ»Ny<10CF»PXV(iG-»dF>FYV(is.v,o),rZV< 
C  »  a  F 

CUMMJN/L lNfc/Ab,dS»CS»  XS»  r J#  Zif  aC»  m.»  cc 
vOMMuN/NOKM/ AXNJ  loo  )  t  A  YN  (  100  )  t  A  ZN  (  *C u  ) 

CuMhUN/L  YL/NCY#nS>^  UC)»NSHC^»iO)»XZ(lu(*ri;(*.  luiliUi«cln  )  *  c 
ioF»CL(10FfkC(iC) 
common/ lnj/T i 
DATA  IH»  XN/.i#i.5W 
DATA  XI / . 5/ 

C  *ALk  IhKOUoH  CYLlNDkiCAL  iOkFAk.E 
r  i  - 1 . 

15R-1 

C  cN  r  S  Ukc  aC  c 

v.A*~(Ad*AC+Bj*UC'*Ci*CC  ) 
i F ( C  A . GT . 1 . F  C  A  >  x  • 

I  F  ( v.  A  .  L  1  .-1.  FCA--1. 

SA-bGnMi.-CA***  ) 

A  A  P  *  j  A / XN 

CAP- jwkT l 1.-JAP**Z ) 

a-acosicaj 

AP*ACJS(LAPF 

(.ALL  ri<ANS(XN»A»AP»  I  1  ) 

C  PASS  TO  UTHlr  SLKPACt 
B-CA/XN-C AP 
AR-AC*B* Aj/XN 
BK«Bv.*B*oj/Xf't 
tK-CC*b+Ci/XN 

AA«Ac(iCF*AR+dt«lCF*BR+CL(iCI*Lk 

BB-AEUCFMXi-XClICFF+BEdLFMYS-YCULFJKeFiOFMZWCdCFJ 

A 1 1 A  l ( 1 C ) ♦ AA-Ak 

A2-BE ( iC ) * AA-BR 

A J-Cc ( 1C ) ♦AA-CK 

BL-XC ( it F-XS+BB* Ac ( 1C ) 

B2-YC  l  1C  F-Y5*BB*BC  (  1C  ) 

B3-ZL(IC)-ZStBB*Ct(lC) 

A*  AJ.**Z*AZ**Z*A3**Z 
o»A1*B1+A2*62*A3*b3 
C«Bl**2fB2**2+b3**2 

AB- ( SOk  T( B**2  + ( <RC  <  1C  »  +TH1 **2-C ) *A) ) /A 
bA«-a/ A 
R  -  b  A  ♦  A  B 

lF(k,LL.j.lK«BA-AD 
XK- XK*  Ak*K 
YR«YR*bk*K 
Zk-Zk+CR*R 

Ku-K,*AAtdB 

XO-XC ( IC ) *Ac I  1C ) *RC 

YC-YC(ICF»Bt(lCF*Ro 

ZO-ZC(iCF*CE(IC)*RJ 

ml- ( XO-XR F /kC (  iC  ) 

BC*(YG-YRF/KC(ICF 
CC ■ ( ZO-Zk  F  /RC (  IC  F 
T  l  ■  1 1  ♦  t  X P  (  —X  r*R  » 

C  cXir  iUkF  ALL 

LA--(At'*aCxbR*BC«’CR*Lc  F 
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iFICA.GT.  A.  )L A»i. 

1FICA.LT.-1.  IU»1. 

iAP-XN^iA 

ifuap.gt.a.igotu  iJ 

C  KAY  fcXITS  F K Oh  SURFACE 
CAP»GJKT l l.-SAP*»t) 

A-ACUS(CA) 

AP«ACl)S(CAP  ) 

CALL  TkANS(XN#A»AP»T1  ) 

O*  XNACA-L ap 
ar»ac*b*x^aak 
br»  bc*b  +  xnabr 

CR«CC*b+XN*CK 
Kc I OkN 
AC  C  UN  I  I NOE 

C  KAY  CONTAINtO  WiTHiN  i  ok  F  AC  c 
T  I  *  -  . 

Ak*  XN* ( AC*CA+AR ) 
i R  »  AN* ( Be  *  C A  +  6K ) 

Ck»aN*(CC*CA+Ck» 

HE  1  U  K  N 

cNo 

jUBkoJT  INL  TRAnS(hN»A»AP»  I  ) 

C  FRcjNtL'A  L  A  n  FUK  ON  P  UL  AK  I  i  1 0  LIGHT  AT  0  I  c.  L  L  o  1  k  l  o  AokFACi. 
F  «  K  N** i 

F1«SkKT<F-C>IN(A) )**2) 

F  ^  *  C  JS ( A  ) 

RPa((F2-FA)/lF2*Fl))**2 
R$»  l(F*F2-FA)/(F*F2*F a)»**2 
k*.3*IRP+PS) 

T«T*(l.-R»*Fx/CCB(AP) 

R  t  T  URN 
tNO 

iUBKOoT INE  ORwPI AW/NVtL ) 

COhhON/FAC  T/$X» 4j< bJ#Al#bi 
C  OPTICAL  OliTORTIONi#  ZERO  To  INFINITY 
C  T »  TRANSMISSION#  UNITY  To  ZtRU 
RcWlNO  4 
DO  A-  I  *  1  #  N  V  LL 
RcAD(4»10uC)FI»XK#Yk#ZK 
RcAUUt  lOOOD 
LOCO  FURMAT<4FAC.4» 

I F ( F  L • L  T • • COCCI .  OR  .  F  I .  (»1 •  a  • )FAa.C00wA 
lF<lfl.cG.I)iC«-ALGG(FI ) 

IF  C IU.LU.2) IC«ALOG(Utll 

IFUC.GT.lsUOAS 

XXaCYR-b.5)*SXtiOO. 

VVatZR-lll.b2l*$x«iOG« 

CALL  P  L  T  0  T  S  < i#lC#XX, YY,i,C) 

1J  CONTINUE 
Rt«ilNO  4 
00  2v  I ■ I #  NW L C 
RfcA0(4#iwG0)Fi#xK#YK»ZR 
K;  AJ ( 4# 1000)0 

IF(FI#LT«<oOO'jI<uR»Fl»GTai<  tFI*aCov.Jl 
IF(  lU.tU.l  )  1C--ALUGIF1 I 
IF ( ia.c0.2)iCaAL0G(0+lJ 
I F  C  IC.GT.i5J  iC-ii. 
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Xx«  i  Tk-5  .  s  /♦Sa*1oO. 

YY*( XR+104. )*SX*433.44 
CALL  PLT0TSI3,  IC,XX, YY,1,C) 

20  -liNT  INUt 
REklNO  4 
OU  30  I • 1 #  NVLL 
REA0(4,lG0G)Fl»XR,YR,ZR 
Rt  AO  ( 4»  1000)0 

If (Fl.LT  ..jOOOI.OR.  FI.  GT.l  .  )H -.GuCCi 
If ( IG.EQ.niC— ALOUlFI  ) 

IF(  IQ.L0.2  )IC«ALG(,(D*1  > 

IF(iC.GT.15)lC*l- 
XX*(ZR-111.52>*3X»690.42 
YY* UR*1G4.  )*SX*433.44 
CALL  PLTDTSl3,IC,XA,YY,l,C) 

30  Continue 
r  t  r  jr  4 

t  NO 

lUdRuuTINt  QRA» 

C  OkARS  GRAPHIC  PICToRt  OF  CANOPY  IN  THRtt  FOLD  LATuCT 

COMMON /ORA /MX l 10U),N3 (10Q)*XS(100u),YS(l3:C),NSY(lC:),NtxfA 
COMMON  /  CAN/NT»Nb»NC#NP,NA,Nu»NV(lot  I,  PXV  I  iu^c  )>  rl  v  i  K  j.d  I  >  i-Zv  I 

0*3) 

COMMON/ NORM/ AXN( »0o) , a YN( it. ), A2N  <  -V.  0  ) 

CuMMUN/PILUT/XC,YC, ZC 

tOMMGN/FdCT/SX>A^#tJd»Ai»di 

NA*0 

3X«2.122 
NS  S  *  NC ♦ 1 

rtG*-&, i 

du* 100 . 

Ai*-iix.52 

ol«i00. 

XX* ( Y0*At  >*3X*Bt 
YY*  Ud  +  Al  ) *S  X+  B 1 
CALL  d£LT« 1, XX, Yr»l,U) 

CALL  t  MARK (  X  X,  Y Y  ) 

N  E  *  1 

Ou  2  1=N5S»Nl) 

U  5  *  ~  A  X  N  ( 1  ) 

1 F ( (I.LL.nC.ANo.MS.LT.G.) •  0  R  •  t l.GT.NC.AND.US.tl . v .) )  G  o  T  0  c 
C  SURFACE  fACtS  VUtatR  FROM  Slot  VIE* 

XX«(PYW(I»1)+AG)*SX+B0 
YY«IPZVU*ll*AU*SXm 
NE»Nt*l 
MO*  1 

if ( 1 .LE.Nt.LR.l •  t  G . N A ) MO  *  4 
Call  b  EL  T ( Nt  »  X  X, Y Y»  MD»  o ) 

CALL  LINS(I) 

2  CONTINUE 
A  1  *  1 04 . 

31*433,44 
XX*( Y0*A0)*SX+B0 
YY*IX0+Ai)*5X+Bi 
Nt*Nt>  1 

CALL  BELT (Nt,XX,YY,i,U) 

CALL  tMARKI XX, YY  ) 

OU  3  1»NSS,NL) 
oT*+AZN(I) 
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if  (  (  1  a  L  t  a  NC  a  AND  •  0  I  a  L  (  aCa  )  aUK  a  (  i  a  la  1  aNCaANJawl  a  T  a  w  a  I  )  G  u  T  0  3 
C  SURFACc  FACES  VIEWER  FROM  l OF  vie* 

XX-  (PTV  ( If  1)*aCi)  •SXfttO 

YY-(PXV(1#1»*A1>*SX+B1 

Nc-NE»1 

MO-1 

IF  (  I.LEaNCaCRaI.EQ.NA»HO-<* 

CALL  BELT(NE»XX,YY#MD,C> 

CALL  LINT (  I ) 
i  CONTINUE 
AD--111.S2 
B0-690. AZ 
xx» ( *0* AO ) *3A»8c 
YY- ( XO*Ai  )*SX*Bi 
NE ■ N£ ♦  1 

CALL  BeLT(N£*XX*YY;1*G) 

CALL  E  TIC ( XX#  Y Y  J 
Ob  1C  I-N$S#NO 
0F--AYM1) 

IF  (  (  1  aLEaNCaANOaUF  a  L  T  aCa  )  a  0  R  a  (  1  a  G  T  a  N  C  a  A  N  0  a  0  F  a  U  T  al'a  I  )  C  U  T  U  1  > 
C  iURFACt  FACE}  VIEWER  FROM  FRONT  Vic* 

XX-(PZv(1«1)#aO)*SX*bC 

YY-(PXV(1#1)+A1)*Sa»b1 

NE»NE*1 

M0-* 

1F( l.Lt.NC.0K.I.cU.NAIND«4 
CALL  BELT (Nc#XX#YY#Mu»C) 

CALL  L INF ( 1 ) 

10  CONTINUE 

NS( NE*i) -NX+1 
Kc  T  UK N 
cN  0 

SUBROUTINE  BElTI  1 e# x , Y , MuD# NS  3 ) 

COMMON /DR  A /HD  l  iCt  )  #NS(  10.  )  #  Xi(  100C  »»YS(10CC)»\iYUCC»#Nt»NP 

NP-NP+i 

NS ( I c • ■ NP 

NSYdEJ-NSS 

M J ( IE ) -MOD 

XS(NP)-X 

Y } (NP  )  ■ Y 

RETURN 

END 

SUBROUTINE  PLEM(XjY) 

C  OMMUN/ OR  A/MDI  10c  )#N}  (100  )  jXSUGDC  ),y  SUC  DC  I  #  NS  Y  ( lOc  >  »  NE  *  NP 
NP- NP ♦ 1 
X  S  (  NP  )  -  X 
YS(NP J-Y 

return 

cNO 

SUBROUTINE  EMAKk(X#Y) 

C  MARK  EYE  POSITION  IN  CANOPY 

x-x-ia . 

Y-Y-S. 

CALL  PLEM(X#Y) 

X-X+10. 

Y-Y+Sa 

CALL  PLEI1(A#Y) 

X-X-l j. 

Y-  Y  ♦  5  a 
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call  PLLnu*n 

X  *  A  ♦  i  , 

CALL  PltNlX,Y) 

Y* Y *1. 

call  PLtnu*n 

i 

CALL  PLtM(A,Y) 

Kt I UFN 

c  N  J 

juoKUUllNt  tTIC(X»Y) 

Y*  r-i . 

CALL  PLfcM(X,Y) 

Ys  Y  ♦  u • 

CALL  PLEM(X,Y) 

Y  *  1  -  >  , 

ull  p l t n  t  x »  y  ) 

X  *  A  ™  ‘J  • 

lall  P l  t  m i x»  y  ) 

X«X*aU. 

CALL  PLtttU,Y) 

K  c  I  UKN 

NU 

juoKU  UTiNc  L  I  N  S  (  A  i 

COMMCN/C  AN/M  *  Na»NC»NP»NA,Nu»NV  (  LOC  J  »PXV  ( xu.  >  C  I  .  P  YV  (  LC0»  t  I  ,  PZV  (  luC 

J  t  3  ) 

CuMMuN/PACT/SX, AC»BJ» A1»B1 
N  A - N  V (  I  ) 

Ou  ac  A*a»NA 
A  1  *  A  ♦  1 

IP  l  A.cU.NA)Al*i. 

XA=»PY\i(I»Kl)  +  AO)*iX  +  BCJ 

Y  Y* l PZn ( 1« A1 JtAl )  ♦  S  X  ♦  B  i 
LALL  PLtM(XX.YY) 

10  L JNTINUt 

At  i  UK  II 
t  NC 

jLiokuo  flNt  LlNF(i) 

UUMHUN/C AN/NT, NB»NC»NP»NA, NU»NV(*0C  »«PX\.(*OC»d)#PYW(*CL.bJ#PZV(lu- 

0  t  >  I 

COFMwN/F  AC  T/Sx*  AGfBCs  A  1  /  B1 
NA*NV (  1  ) 

Ou  l'j  A»i,NA 
A  l  *  A  ♦  i 

IF ( A.fcw.NK  )K1»1 
XX*  <PZV<  I.KU  +  au  l*$X*BL 
YY*(PXV(I#a1)kAi)«SX+Bi 
CALL  P  L  t  M  ( XX >  Y  Y ) 

12  u  Uis  T  1  NU  t 
«L  TUAN 
cNJ 

judAOui iNt  l ini ( i ) 

CUMIuN/C AN/NTt NB*NC>NP>Na,Nu,NV  ( 100  )  *  P  XV  (IOC»c)>pYVl  1  OC  »  e  )>PZV  (  l-’v 

u»  o  ) 

CunhuN/FACT/jX»AC»B0»Al»81 
NA  *  N  V (  I  ) 

DO  ac  A  * i »  N A 
A  A  *  A  K  1 

IF (A.cU.NA)A1«a 
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XA«<PYv(I*Kl)*Au)ASX4BC 

yy- «  pxv  < i*ri )tAi )*sx*bi 

CALL  PLEM(XX,YY) 
lu  CONTINUE 
RETURN 
END 

SUBROUTINE  ORm(IUNLT) 

C0MHQN/DRA/HU(*00»*NS1 100)*XSI  10U0)  .  rS  (1CCO)*nsY  (IOC  l,m»NK 
COMMON/ L 8/ STKG ( A ) 

COMMON/ DG/XPIl 00 )*YP(  100) 

DIMENSION  LABEH4J 

DATA  LABEL/"SMYTH","b520"*,,A36S4,,*"Ht  379"/ 

C  PROGRAM  COMPUTES  LAYOUT  ON  1924  UNITS  FITTEL  lu  ij  iNuHtS 
C  BY  25  INCHES  DISPLAY  SHEET 

CALL  PLTBtGT25.*25.,l.,i UN*T, LABEL) 

CALL  PL T SC  A <5. *3. *20-.*  * . * 4U . 98* 40 . 9b ) 

CALL  PLTSYMt .25*STkG*u.*201.**.  ) 

DO  10  1  £  ■  1  *  N  E 
MDD-MDt  It » 

NSS-NSY (  It  I 
Ni-NSCIE) 

N2-NST l  E ♦ 1 1 -N 1 
00  5  I-1»N2 
XPl  1)-XS(N1+1-1) 

YP(I )-YS(Nl*I-i) 
b  CONTINUE 

CALL  PLTuTS(M00,NSS*aP,YP,N2,J) 

10  CONTINUE 
RETURN 
ENO 

C  A AH  CANOPY  DATA 

c  numb t r  of  vertices*  surfaces— constraint*  flai*  anu  cylindrical 


1 6o 

lb 

5o 

58 

39 

64 

NO. 

VERTICES  1 

PER  SURFACE 

•* 

4 

4 

4 

4 

4 

4 

4 

4 

4 

5 

3 

4 

4 

4 

4 

4 

4 

4 

A 

4 

4 

4 

4 

7 

7 

7 

7 

4 

A 

6 

4 

4 

4 

4 

4 

4 

4 

5 

4 

4 

5 

7 

7 

/ 

7 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

0 

6 

b 

6 

b 

4 

VlRTICE 

ASSIGNED  TO 

EACH 

SURFACE 

J  89 

090 

C  66 

121 

125 

129 

133 

135 

*39 

*43 

145 

149 

153 

157 

137 

145 

DC  * 

002 

C„3 

0  v  4 

0u7 

0C1 

000 

012 

Oil 

021 

015 

028 

022 

032 

029 

033 

lie 

035 

036 

111 

038 

039 

038 

109 

U4C 

041 

112 

053 

047 

06C 

054 

061 

C61 

069 

094 

093 

100 

C  99 

098 

087 

C  77 

C  8  1 

101 

066 

071 

094 

ICO 

u85 

117 

11 A 

113 

122 

126 

130 

i  34 

136 

1 V* 

*  44 

14b 

15C 

154 

138 

136 

148 

C02 

CJ3 

Gu4 

003 

Ul4 

007 

004 

C  14 

012 

0  2- 

C  16 

027 

0  23 

v  3 1 

030 

034 

109 

G3o 

037 

112 

039 

03o 

111 

111 

041 

11C 

110 

C52 

048 

959 

0  55 

u64 

062 

u  9c 

089 

094 

095 

100 

099 

098 

C  76 

082 

102 

..69 

072 

093 

099 

C  86 

114 

115 

119 

123 

127 

131 

i  30 

137 

*41 

*40 

147 

131 

166 

159 

162 

163 

009 

010 

0*1 

012 

012 

005 

003 

Cud 

U'09 

u  19 

0*7 

026 

v  2  4 

..30 

130 


031 

109 

033 

043 

044  044 

046 

034 

1C  9 

045 

112 

042 

051 

049  056 

056 

C63 

C63 

0fe6 

065 

0  70 

071 

076 

U  73  074 

079 

083 

103 

*7: 

g  7  3 

092 

096 

ub  7 

w  4  G 

w  G  6 

065 

124 

126  132 

129 

138 

142 

±39 

l«»d 

132 

155 

160 

161  164 

006 

C09 

OxC 

DU 

GG  5 

Ol  2 

002 

009 

010  018 

016 

025 

023 

032 

110 

042 

042 

043  037 

043 

033 

034 

*36 

04b 

C  4  5 

C43 

05u 

050  057 

057 

062 

u  64 

Oo^ 

u  70 

k64 

0  76 

0  73 

074  C  86 

Obk 

104 

v66 

074 

j  9  X 

097 

ceb 

j 

L 

0 

0 

C  \J 

0 

0 

u 

i 

0 

lo  6 

136 

0 

0  0 

c 

0 

c 

r 

j 

w 

0 

0 

j  J  1  7 

0  a  9 

0  Zb 

J 

v 

04 1 

L 

C  J 

x 

j 

* 

L  23 

v 

u 

044 

044 

051  g56 

o5c 

0 

w 

k 

J 

C 

C 

0  C 

G 

o 

l.G 

t  66 

0  75 

O  *y  w 

u9o 

0 

g 

w 

L 

G 

0  i 

Z 

3 

J 

C 

V 

G 

0 

W  * 

0 

w 

«.* 

. 

w 

G 

0 

0  01b 

020 

C23 

g  2  7 

0 

V 

to 

t 

0 

C  V 

c 

C 

u 

046 

052  055 

C  3  9 

t 

0 

J 

G 

V 

j  : 

V 

/ 

iv.b 

*67 

C  76 

j  fc  0 

w95 

r 

* 

* 

0 

i* 

0  0 

t 

c 

C 

V 

u 

>* 

V 

L  w 

w> 

w 

V 

w 

0 

C 

C  013 

l21 

022 

JLC 

J 

0 

c 

v* 

G  0 

C 

rj 

VP 

C 

0 

047 

053  054 

060 

V 

J 

w 

•- 

c 

w 

0  0 

W 

'J 

iw  ( 

0 

O 

J 

0 

<J 

j 

J 

L 

0 

w  U 

0 

c 

J 

Sj 

O 

0 

c 

0 

C  0 

0 

L 

V* 

u 

V 

0 

0 

k  C 

»J 

V 

•„ 

■* 

0 

U 

c 

0 

0  G 

J 

0 

w 

* 

0 

C 

0 

0 

0  0 

0 

0 

3 

Sp 

c 

V 

0 

0 

0  o 

si 

0 

1j8 

0 

0 

0 

0 

0 

C  VcKTiOt  X 

-POSITION 

2 

•  C  <j 

2. 

oc 

2 . 0  C 

2  .  C  G 

2.00 

2  •  g  . 

6.0k 

"  w  #  -  «p 

-2 

•  00 

-2. 

Oo 

-2.  Do 

-2.0g 

-2.0. 

-2.x: 

19.06 

1  e .  7  0 

.23 

12. 

90 

li.20 

14.76 

il.fi 

-  X  C  •  1  G 

—  A  7 . OG 

-1  3 . 6  ■» 

-13 

.30 

-13. 

24 

-15.90 

17.96 

0 . 3k 

—  1  A  .  3  G 

-12.04 

9  .  g  C 

19 

.6  3 

19. 

63 

11.46 

1  i  .4b 

11.46 

19. g3 

19.  bi 

-x9.c * 

-19 

.63 

-11. 

4b 

-11.4b 

11. 4b 

- 

19.63 

-a  *  •  t  3 

X4.lt 

Ib.b. 

14 

.22 

13. 

iC 

11.24 

x  4 , 7  o 

1  7 . 3  U 

—  1  7 . 9  g 

-17. bO 

—  1  4  .  G  k 

-13 

.30 

-13. 

26 

-13.96 

x  7 . 9b 

3.5b 

—  a  *  «  4  6 

-12.14 

4  .  A  C 

23 

.06 

2  0. 

9c 

18.06 

x  3  •  7  2 

22.72 

2  3.2c 

23.20 

1  3.6g 

13 

.00 

13. 

16 

1  b .  3b 

26.64 

x  6 . 7  2 

-lb.  72 

-14.2c 

14.2c 

13 

•  oO 

-13. 

60 

-9.2g 

9.2g 

10.56 

—  aI . 5  b 

-1a  .  56 

1 .  .  3  6 

-23 

.Jb 

-20. 

92 

-18.08 

13.  f  2 

- 

22. 72 

-2  J . 2  o 

-23.20 

-  1  J  .  C  k 

-13 

.00 

-15 . 

16 

-lb.3o 

■22 .04 

a9. 63 

13.2b 

1a. 46 

-  1  A  .  4  C 

-13 

.  26 

-19. 

b3 

-12.00 

12.00 

13.26 

-11.46 

a  J  .  26 

-1  1 . 4t 

-23 

.75 

-4 

• 

4. 

23.75 

-4  . 

-23.73 

•. . 

c3. 7. 

if 

• 

3o . 

3b. 

2<i  . 

- 

5b. 

-24  . 

-  c.  4  , 

-3c  . 

Zb 

• 

104  . 

104. 

2t . 

26. 

1  G  4  . 

6  7 . 6 

t  A  . 

ol 

• 

67. 

6  -it4. 

26. 

- 

26. 

-  a04  • 

-  X  v4  , 

~7t. 

-bl 

• 

-67. 

6 

-b7.6 

•61  . 

- 

4o. 

-24  . 

-2  2. 

—  5  G  . 
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29  . 

96. 

56. 

22. 

61. 

61. 

-61. 

-61 . 

c 

VtrtTlCt  Y- 

POSITION 

37.50 

57. 5C 

08.37 

71.57 

57.50 

68.37 

71.57 

56.98 

8  3.96 

66.96 

96.18 

98.66 

86.98 

98.16 

98.00 

96.29 

98.95 

103.99 

100.92 

110.92 

103.99 

100.92 

110.92 

115.61 

192.59 

196.  OC 

157.52 

158.93 

19o . 0 J 

137.52 

136.90 

156.32 

39 . 2D 

59.20 

70.19 

108.92 

139.09 

156.97 

.50.97 

19C  .09 

69.57 

69.57 

lu6 . Oi 

iOo.Ci 

59.20 

5  9.23 

7d.19 

160.92 

1  39 . o9 

156.97 

156.97 

19*  .07 

115.34 

115.52 

113.01 

113. Cl 

37.5 

3.5 

v  » 

>  (  3 

57.3 

39 . 2 

65.2 

.32.6. 

1 5  o  .  6 

iab.Oi 

191. ox 

226.6. 

230.61 

1 8  0 . 6 1 

*  6  8 . 6 1 

19. .61 

180.6. 

i  09  .  o  1 

109.61 

160.61 

1  8  0  .  t  . 

219.01 

219.61 

219.61 

219.61 

198.61 

229.61 

.96.6 

229 .61 

c 

VtKlICt  L- 

POSHION 

131.81 

138.66 

196.27 

191.77 

1.38.66 

190.27 

191.77 

133.9/ 

195.26 

198.59 

198.06 

136.39 

198.59 

196.08 

136,39 

12  7.9o 

129. 20 

198. CC 

157.38 

1 60  •  51 

1 9  o • OC 

157.56 

1  b\j  .  3 1 

1 60  .  t  o 

lo9 . 26 

107.69 

167. 6t 

136.26 

167,69 

167.08 

1 3  6 . 2u 

197,39 

129.05 

191.91 

139.30 

172.90 

179.26 

179.29 

176. 9« 

.  3  v  .tic 

156.72 

15b. 72 

.73.13 

173. *3 

129.05 

191 . 91 

.59. 8v 

172.96 

179.26 

179.29 

176.00 

1  3  C  ,8  <j 

132.59 

152.59 

175.79 

173.79 

199.2 

129.1 

129.1 

199.,; 

129.12 

129.12 

129.12 

129. *2 

192.52 

190.52 

190.52 

192.52 

111.52 

111.52 

190.52 

192.52 

132.52 

132.52 

111.52 

.11.52 

197.92 

197.92 

176.52 

176.32 

111.52 

132.52 

111.52 

1 o2  •  52 

c 

CYLINDRICAL  DATA 

1 

3 

2  2 

1 

60  61 

62  63 

69 

-31.9959 

.35.89 

197.9713 

•  C 

31.9959 

135.89 

197.9713 

.0 

C. 

133.6172 

139.3729 

.0 

c 

PHJT  tYL 

POSITION 

-1.5 

192.33 

17. .2 

C  CANOPY  EDGES  FlirtO  TO  CYLINDRICAL 
20.783  17.875  130.91 


-10. 

1  0  • 

10. 

-K  . 

-36. 

6  • 

38.98 

5  8.98 

Ov .  30 

5  7.56 

3  o  .  96 

6  *  ,  3  6 

69 . 06 

82.92 

96.29 

3  9.06 

Ot  .  92 

8  3 . 9  o 

8  7.98 

6  7  .  9  0 

9  7.50 

97  ,,6 

113.61 

1 2 1 .90 

121.50 

98.9  3 

121.90 

121.30 

199.20 

19  6.32 

150.32 

199 .26 

x94. 52 

192.59 

196.0; 

l  9  c  .  w  y 

x  3  6  •  6  0 

i  5  o  .  fcc 

110.96 

C9.ro 

110.99 

1  13.43 

3  8.16 

5  0 , 4  0 

0  7  .  *  6 

t/.j" 

112.13 

112.13 

195.99 

x  9  5 . 9  i 

ilO. 90 

tv .  t  c 

116.97 

113.23 

1x5.52 

1  1  5  .  3  2 

119.82 

i . 9 . 32 

163.99 

1  »  3 . 9  9 

121.96 

1 2  x  .  9  c 

3  .  5 

3  7.3 

✓  •  ^ 

3  7.3 

65.2 

5  9.2 

13  8.6 

1  32  .  t  i 

I86.0I 

191.61 

1  6  9 . 0  1 

1  8  9  .  C  l 

2  3  o  .  o  x 

2  2  6 . 0  X 

1  9  x  .  1 1 

1  o  6  •  0  1 

219.01 

219,0* 

2.9.61 

2*9 .01 

198.01 

198.01 

X9fa . cl 

198.0* 

2  *  9  .  o  1 

2  i  9 . 6  1 

x  33 . 7  7 

135.97 

152.83 

121.8,. 

133.97 

*32.03 

X  4  6 , 9  r 

186 . 39 

127.96 

1 2  r  .  •,  0 

1  30 . 5  9 

x  9  5 . 4  0 

1x6. ol 

118.01 

136.00 

1  3 1 .  *  C 

io*. . 6d 

x  2  9 . 2 1. 

129.21 

1  L.  9  1  4  * 

139. 

129,*. 

197.36 

199.96 

19  7. 59 

197.30 

x  9  9 . 9  w 

It  9 ,2o 

x37.8o 

13  7 .06 

173.08 

1  7  3  .  C  8 

150.68 

1 3C . c  5 

196.36 

179.39 

193.32 

193.-2 

139.90 

139.9 1 

177.93 

177.9.1 

182.12 

1  6  c  ,  x  2 

153.08 

130.83 

1 9c . 36 

179.59 

132.59 

x  52 • 59 

*39.03 

1  5  9  .  *  J 

198. OC 

1 98  . *  * 

189.2* 

139.2* 

122.1 

1 2  7  •  C  - 

122.x 

127.13 

149.12 

1  4  9  ,  X  2 

149.12 

24  9,  *  4 

132. 12 

132.92 

132.52 

132.32 

192.52 

x9u  .5  ’ 

132.72 

132.12 

197.92 

197.92 

17b.32 

x  7  d  ,  3  2 

161.32 

lc 1  •  52 

221.52 

221.52 

It  7.92 

137.92 

.960  5 

•  196c 

53.9916 

.98)5 

•  *  96o 

55 . 99ic 

.983  7 

-.179o 

92 . 9 i7i 

SIDE  l 
iti.175 
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.9592 

.0 

.2820 

<>•90  .  9991 

11.25 

0. 

126.175 

100.131 

w  • 

-.1796 

.9837 

j . 

CANuPY  OiiPLACthtNT  OF 
u  . 

CTL  INORICAU 

SlJtb  f K L rl  i-L*.I  kuNOuw 

C  luP  talNDGto 
1  •  2 
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